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\  ABSTRACT 

\l 

Two-dimensional  steady-state  solutions  (suitable  for 
computer  evaluation)  are  given  for  moving  normal  contact 
loads  (sliding  or  rolling)  on  the  surface  of  a  semi-in¬ 
finite  viscoelastic  base4  wltl^otttnre^t'r^tion 
profiles  (e.g.  circular  cylinders)  or  to  the  simplest 


vlscoelastlc-'me'dla  (e.g.  a  single  relaxation  time,  same 
behavior  in  shear  and  dilatation) . 

-^Examples  worked  out  Include  (a)  the  rigid  cylinder 
rolling  on  a  material  haveng  5  relaxation  times;  (b)  the 
flat  punch  with  corners  rounded  to  eliminate  infinite  pres¬ 
sures.  The  results  for  (a)  show  that  the  variation  of 
rolling  resistance  over  a  range  of  velocities  becomes  smooth¬ 
er  with  the  greater  number  of  relaxation  times.  The  results 
for  (b)  show  the  dependence  of  the  pressure  peaks,  near  the 
ends  of  the  contact,  on  the  nature  of  rounding  of  the 
corners  and  the  tilt  of  the  punch. 

Three  dimensional  solutions  are  considered  for  the 
purely  viscous  material.  It  is  shown  that  there  is  a 
steady  state  for  a  base  of  finite  thickness,  but  not  for 
the  semi -infinite  base.  ' 
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INTRODUCTION  AND  SUMMARY 


Distributed  loads  moving  over  the  surface  of  a  de¬ 
formable  material  are  of  frequent  occurrence  in  engineering 
Many  technical  operations  involve  a  moving  roller  or 
slider.  Impact  of  a  high-speed  Jet  on  a  moving  surface, 

and  certain  cases  of  explosive  loading  and  oblique  impact 

1 

equivalent  to  such  Jets  ,  provide  soir  of  the  less  obvious 
examples.  When  the  velocity  of  the  load  is  small  compared 
with  stress-propagation  speeds,  inertia  may  be  neglected 
and  the  problem  treated  as  quasi-static.  In  this  case,  if 
the  material  is  elastic  the  solution  is  the  same  as  for  the 
stationary  load  problem,  but  now  carried  along  unchanged 
with  the  load.  If,  however,  the  stress  depends  not  only 
on  strain  Itself,  but  also  on  time  derivatives  of  strain, 
the  solution  to  the  moving  load  problem  has  no  such  immed¬ 
iate  relation  to  the  corresponding  statiot  ary  load 
solution. 

In  this  investigation,  a  steady  localised  load  moves 

at  constant  velocity  in  a  straight  line  on  the  initially 

flat  surface  of  a  semi -infinite  base.  The  base  material 

is  linearly  viscoelastic.  Isotropic  and  homogeneous,  with 

]_  ~  ”  2 
See,  for  example,  Abrahamson  [1]  . 

2 

Numbers  in  brackets  refer  to  the  Bibliography. 
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stress-strain  relations  expressed  in  hereditary  integral 
form  (or  as  differential  equations  in  special  cases).  It 
is  free  from  body  force,  at  constant  temperature,  and 
initially  at  rest  and  undisturbed.  Inertia  effects  are 
neglected.  The  problems  are  considered  with  the  usual 
assumptions  of  linear  theory:  infinitesimal  strains, 
boundary  conditions  applied  to  the  original  undeformed 
surface. 

The  principal  features  of  Interest  are: 

(i)  the  possibility  of  a  steady  state  when  the 
load  has  been  moving  uniformly  for  a  sufficiently 
long  time; 

(li)  the  solution  for  such  a  steady  state; 

(ill)  the  rate  of  energy  dissipation  in  the 
viscoelastic  material. 

A  consequence  of  energy  dissipation  is  that  work  is  re¬ 
quired  to  maintain  the  uniform  velocity  of  a  moving  body 
applying  the  load,  and  the  moving  body  encounters  a 
resistance  (e.g.,  rolling  resistance). 

The  general  problem  of  a  moving  load  on  a  semi- 
infinite  base  is  treated  in  Chapter  I.  It  is  shown  that 
a  steady  state  is  always  reached  if  the  base  is  a  solid, 
but  there  is  no  steady  state  if  the  base  is  a  viscoelastic 
fluid.  General  expressions  for  surface  displacement  are 
developed,  and  some  examples  given.  Moving  loads  on  a 
purely  viscous  material  are  considered  in  Chapter  II. 
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Specific  results  are  given  for  a  uniform  circular  load, 
showing  the  change  in  surface  displacement  with  time. 
Although  there  is  no  steady  state  for  a  semi-infinite  base, 
it  is  shown  that  for  a  base  of  finite  thickness  thei*e  is 
a  steady  state,  its  character  depending  on  the  way  in 
which  the  lower  boundary  of  the  base  is  attached  to  its 
support . 

The  rest  of  the  investigation  is  concerned  with  the 
steady  state  moving  load  problem  resulting  from  the  moving 
contact  of  a  rigid  body  of  given  shape  with  the  surface  of 
a  viscoelastic  solid.  The  general  theory  of  such  contacts 
is  considered  in  Chapter  III.  Further  detailed  analysis 
is  limited  to  two-dimensional  (plane  strain)  problems.  The 
particular  problem  of  a  rolling  rigid  cylinder  is  considered 
ill  Chapter  IV.  Some  results  available  in  the  literature, 
for  materials  of  very  restricted  behavior,  are  examined 
and  their  limitations  discussed.  A  numerical  method  for 
solving  two  dimensional  moving  contact  problems  is  devel¬ 
oped  in  Chapter  V.  '.’he  method,  suitable  for  computer 
evaluation,  can  be  used  for  very  general  contact  profiles 
and  any  viscoelastic  solid.  Solutions  are  given  for  sev¬ 
eral  rolling  cylinder  examples,  illustrating  the  results 
for  various  materials  over  a  range  of  contact  times. 

Results  for  a  material  with  five  retardation  times  show 
that  viscoelastic  creep  effects  are  much  less  pronounced, 
and  the  rolling  resistance  vs.  velocity  curve  much  smoother, 
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than  for  a  simple  material  with  a  single  retardation  time. 
These  results  indicate  that  a  single  retardation  time  does 
not  adequately  characterize  the  behavior  of  actual  materials 
in  moving  contact  problems.  Illustrative  solutions  are 
given  also  for  a  flat  punch  with  corners  slightly  rounded. 
The  pressure  is  everywhere  finite,  but  there  are  sharp 
peaks  at  the  ends  of  the  contact.  Results  indicate  how 
the  pressure  peaks  depend  on  the  nature  of  rounding  of  the 
corners  and  the  tilt  of  the  punch. 
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CHAPTER  1 


GENERAL  FEATURES  OF  MOVING  LOAD  PROBLElvIS 

The  quasi-static  proble.,i  of  an  arbitrary  load  moving 
uniformly  on  the  initially  flat  surface  of  a  semi-infinite 
linear  viscoelastic  base  is  considered  in  this  chapter.  In 
rectangular  coordinates  x*,  y*,  z*  fixed  in  the  viscoelastic 
material,  the  base  is  taken  as  z*  >  0  ,  with  the  surface 
z*  =  0  .  In  this  discussion,  "load”  refers  to  a  pressure 
y*)  applied  normal  to  the  surface  over  an  area  A. 

There  are  no  tangential  tractions  and  no  other  normal  trac¬ 
tions  applied  to  the  surface.  The  load  is  moving  at  constant 
velocity  V  in  a  straight  line  parallel  to  the  x"*  axis 
(Fig.  1.1).  For  later  convenience  are  introduced  coordi¬ 
nates  X,  y,  z  moving  with  the  load: 

X  =  X*  -  Vt  ,  y  =  y*  ,  z  =  z*  (l.l) 

where  t  is  the  time  elapsed  from  some  initial  point. 


PIG.  1.1 
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1.1  Elastic  Solution 

The  problem  of  a  load  on  the  surface  of  a  semi-infinite 
elastic  solid  is  the  well-known  Eoussinesq  problem^  and  the 
complete  solution  is  readily  available  (see,  for  example, 

Love  [22],  p.  192.).  For  the  quasi-static  moving  load  prob¬ 
lem,  the  results  are  the  same  as  for  the  load  fixed,  but 
carried  along  with  the  load.  The  moving  load  solution  is 
thus  the  solution  for  the  fixed  load  with  x*  replaced  by 
X*  -  Vt,  or  expressed  directly  in  x,  y,  z  coordinates. 

The  elastic  solution  makes  use  of  the  potential  of  the 
pressure  dl^'tribution; 

P(x,y,z)  =  fj  dx-dy'  (l.S) 

where  r*  =  )^4'(y-y '  )^+  distance  from  a 

element  at  (x',  y',  0)  to  the  field  point  (x,  y,  z).  In 
particular,  the  vertical  displacement  at  the  surface  is 


u®  (x,  y,  0)  = 


1-v 


P(x,y,0) 


(1.3) 


where  u^(x,y,z}  is  the  displacement  in  the  z  direction 
(the  superscript  e  indicates  the  elastic  solution),  v  is 
Poisson’s  ratio  and  m-  the  shear  modulus  of  the  elastic 
material.  All  displacements  can  be  expressed  in  terms  of 
P(x,y, z)  and  yp(x,y, z)  dz  ,  and  from  these  the  strains  and 
stresses  can  be  found. 
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In  the  two-dimensional  problem  of  plane  strain,  the 
load  does  not  vary  with  y,  and  extends  infinitely  far  in 
the  y  direction.  The  problem  can  then  be  considered  in  the 
(x,z)  plane,  with  a  load  q(x)  (for  a  unit  lenj;th  in  the  y 
direction)  applied  over  a  segment  of  the  x-axis  of  length 
£.  For  convenience,  this  segment  is  taken  as  0  <  x  < 

The  displacement  u  from  the  initially  flat  surface'^  is 

Z 

infinite  (because  as  a  three-dimensional  problem  the  load 
is  infinite),  so  displacement  must  be  expressed  relative 
to  some  arbitrary  point  on  the  surface.  Prom  (1.2), 


a 

=  ^  y'q(x')log|x-x' 
o 


dx’  +  constant 


(1.4) 


Then, 


(.-)  = 


A/ 

^  J  q(x')log|x-x’ |dx' 


(1.5a) 


where  A  Ls  an  arbitrary  constant.  The  displacement  rela¬ 
tive  to  a  point  x^  on  the  surface  is 

Z 

q(x' )  [log|x^-x‘ I  -  log|x-x' I  ]dx’ 

°  (1.5b) 

^2  alway/  refer  to  the  vertical  displacement  at 

the  surface  z  =  0. 
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1.2  Linear  Viscoelastic  Materials 

The  stress-strain  relations  of  viscoelastic  materials 
are  time  dependent.  At  any  instant  they  depend  on  the  past 
history  as  well  as  the  present  state  of  the  material.  For 
any  particular  material^  this  time  dependence  is  conveniently 
determined  by  a  creep  test.  A  constant  stress  is  suddenly 
applied  to  an  initially  undisturbed  specimen^^  and  the  re¬ 
sulting  strain  is  measured.  The  function  giving  the  strain 
variation  with  time  for  a  unit  applied  stress  is  called  the 
creep  function,  J*(t),  of  the  material^  i.e., 

€(t)  =  J*(t)  for^  a(t)  =  H(t) 

J*(t)  is  necessarily  zero  for  t  <  0  ,  since  the  material 
is  initially  undisturbed;  J*(t)  increases  monotonically  for 
t  >  0. 

A  linear  viscoelastic  material  has  a  ratio  of  stress 
to  strain  that  is  independent  of  the  magnitude  of  the  stress; 


Viscoelastic  materials  will  always  be  considered  ini¬ 
tially  undisturbed,  with  the  load  applied  at  t  =  0;  thus, 
€(t)  =  a(t)  =  0  for  t  <  0. 


H(t)  is  the  unit  step  function: 
Its  derivative  is  the  delta  function. 


=  {?:  t>o  • 
6(t)  =  dK(t)/dt  , 


which  has  the  properties 
6(t)  -  0  for  t  ^  Oj 
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e(t)  =  Gq  J*(t)  for  G(t)  =  Gq  H(t) 


Many  materials  of  practical  interest  can  be  considered 
linearly  viscoelastic,  at  least  for  small  enough  applied 
stresses.  This  has  the  important  consequence  that  the  cumu¬ 
lative  effect  of  several  loadings  is  just  the  sum  of  the 
individual  effects.  The  response  to  arbitrary  time-varying 
stress  is'  then  obtained  by  superposition,  leading  to  the 
hereditary  integral  stress-strain  relation  for  linear  visco¬ 
elastic  materials: 

t 

e(t)  =  J  J*(t-T)  (1.6) 

o 

If  there  are  suddenly  applied  loads,  the  stress  and 
strain  histories  will  not  be  continuous  functions  of  time, 
but  will  have  Jump  discontinuities  when  these  loads  are 
applied.  Treated  in  the  classical  manner,  integrals  such 
as  (1.6)  would  need  to  be  considered  in  several  parts,  with 
the  effects  of  the  jumps  added  as  extra  terms.  It  is  more 
convenient  to  work  with  generalized  functions--the  step 
function  and  its  derivati ves--which  take  care  of  the  effects 
of  jump  discontinuities  without  special  consideration.  This 
approach  can  be  interpreted  as  the  limiting  case  of  rapidly- 
changing  continuous  functions.  The  use  of  generalized 
functions  can,  however,  be  established  on  a  rigorous 
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mathematical  basis^.  If  the  discontinuity  occurs  at  t  =  0, 
the  lower  limit  in  integrals  such  as  (1.6)  is  taken  as  0“, 
to  include  the  effects  of  the  jump.  In  this  way,  initial 
values  of  stress  and  strain  functions  and  derivatives,  taken 
at  t  =  0“  ,  are  all  zero  for  initially  undisturbed  materials 

a. 

Initial  values  at  t  =  0  due  to  sudden  loads  at  t  =  0  are 
then  accounted  for  by  the  generalized  functions.  This  often 
eliminates  the  need  for  explicit  evaluation  of  initial  values 
of  unknown  functions,  thus  simplifying  the  solution  of  many 
viscoelastic  problems.  The  use  of  generalized  functions  in 
viscoelastic  problems  is  discussed  by  Corneliussen  and  Lee 
[5]. 

Viscoelastic  materials  fall  into  two  general  categories- 
sollds  and  fluids .  A  viscoelastic  solid  has  a  finite  upper 
limit  to  its  creep  after  a  long  time.  This  means 
a  finite  value,  as  t  -♦  «  .  A  viscoelastic  fluid  continues 
to  creep  without  limit  after  a  long  time;  thus  J*(t)  » 

as  t  -*■  09  .  The  characterization  of  a  viscoelastic  material 
as  solid  or  fluid  depends  on  the  extent  of  time  appropriate 
to  the  particular  problem.  After  a  long  enough  time,  a 
viscoelastic  solid  acts  essentially  as  an  elastic  solid. 

The  material  exhibits  its  final  elasticity,  with  final  elas¬ 
tic  modulus  E^.  Most  viscoelastic  solids  also  have  initial 
1 

See  Friedman  [15l  for  the  use  of  generalized  functions 
in  mathematical  analysis,  and  for  further  references  on  the 
theoretical  foundations. 
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elasticity.  A  suddenly  applied  stress  produces  an  Instan¬ 
taneous  response,  then  additional  creep  as  t^me  passes.  In 
this  case,  7^  Oj  Is  the  Initial  elastic 

moduxus.  If  there  Is  no  Initial  elasticity,  J*(0'^)  -  0. 
Representative  creep  curves  are  sketched  In  Fig.  1.2. 


It  Is  convenient  to  express  creep  functions  In  non- 
dlmenslonal  form.  An  appropriate  time  parameter,  T  ,  Is 
chosen.  Viscoelastic  materials  may  exhibit  one  or  more 
characteristic  retardation  times,  one  of  which  can  be  used. 
Then,  a  non-dlmenslonal  time  variable  Is  taken  as  t  =  t/T, 

For  a  viscoelastic  solid  with  Initial  elasticity,  a  dimension¬ 
less  creep  function  Is  conveniently  defined  by  J(i:)  =  E^J^tT), 
and  the  final  elastic  response  represented  by 
F  =  Eq  =  Eq/E^.  Then,  the  stress-strain  law  (I.6) 

becomes 


e('r) 


=  k  /  ^ 


(1-7) 
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with  =  1  ,  J(oo)  =  F  ,  T  =  t/T  .  This  non-dimensional 

form  based  on  an  initial  elastic  response  will  be  used  in  all 
that  follows,  unless  otherwise  noted. 

The  discussion  so  far  has  considered  a  one-dimensional 
state  of  stress  and  strain.  In  a  general  three-dimensional 
state,  the  response  of  a  linear  isotropic  homogeneous  visco¬ 
elastic  material  can  be  described  by  two  independent  creep 
functions.  This  is  analogous  to  an  isotropic  elastic  ma¬ 
terial  where  two  elastic  constants  are  sufficient.  The  two 
functions  are  conveniently  taken  as  the  shear  creep  function, 
J*(t)  ,  and  the  dilatation  creep  function,  B*(t)  .  For  a 
material  with  initial  elasticity,  these  are  non-dimension- 
alized  in  terms  of  the  initial  elastic  response  and  a 
common  characterisitc  time  parameter  T: 

J(t)  =  2\i^  J*(tT)  ,  J(0+)  =  1  (1.8a) 

B(t)  =  3Kg  B*(tT)  ,  B(0+)  =  1  (1.8b) 

where  is  the  ir.itial  elastic  shear  modulus  and 

is  the  initial  elastic  bulk  modulus.  If  the  material  is 
a  solid,  with  a  final  elastic  response,  the  shear  and 
dilatation  creep  functions  will  each  approach  a  finite 
limit  as  t  ->■  cxj  .  In  this  case,  J*(c»)  =  l/2ix^  and 
B*(oo)  =  1/3K|.  ,  where  and  are  the  final  elastic 

moduli.  Then, 

12 


/ 


(1.8c) 


j(eo)  = 


B(=o)  =  F, 


where 


Fj  =  >  1  ,  and  F^  =  >  1, 


General  three-dimensional  stress  and  strain  is  described 
in  standard  suffix  notation^,  i,j  taking  on  the  values  1,2,3 
corresponding  to  the  coordinate  axes  x^  -  x,  x^  =  y, 

X3  =  z  . 

Stress  components:  a^j[x,y,z,t)  ;  d  = 

2  1 
Stress  deviator  components  ;  s^j(x,y,z,t)  =  ^ 

Strain  components:  e^j(x,y,z,t)  ;  dilatation  e  = 

Strain  deviator  components:  e^j(x,y,z,t)  =  j  et>^  ^ 

Corresponding  to  (I.7),  the  stress-strain  laws  are^ 


■.'Z 


^^1 j ( ^ ) 


(1.9a) 


Kt)  =  3^  / 


(1.9b) 


The  summation  convention  is  used,  so  that  repeated 
indices  mean  a  sum  of  terms  with  index  values  1,  2,  3. 

2 

is  the  Kronecker  delta  symbol:  6^^  =  0  for  i  j  , 
<5ij  =  1  for  i  =  J. 

3 

The  stress-strain  components  are  of  course  functions  of 
the  "ace  variables  (x,y,z)  as  well  as  time;  this  is  not 
explicitly  indicated. 
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For  an  elastic  material  the  creep  functions  are  just  the 
step  functions,  H(t).  In  this  case,  (1.9a,TD)  reduce  to  the 
stress-strain  laws  for  an  elastic  body: 


e 


(1.9c) 


Mechanical  models,  made  up  of  suitable  combinations  of 
springs  and  dashpots,  caii  be  used  to  represent  viscoelastic 
materials.  The  load-defomation  behavior  of  such  models 
exhibits  many  of  the  important  features  of  actual  materials. 
By  taking  a  large  number  of  elements  in  a  model,  a  partic¬ 
ular  material  can  often  be  approximated  quite  closely.  The 
springs  allow  for  storage  and  recovery  of  energy  and  the 
dashpots  allow  dissipation  of  energy  that  is  characteristic 
of  viscoelastic  materials.  In  fact,  molecular  theories  of 
viscoelasticity  have  been  based  on  this  type  of  model  (see, 
for  example.  Bland  [3]).  For  three-dimensional  states,  two 
such  models  would  be  necessary  in  general,  to  represent 
shear  and  dilatation.  The  advantage  of  using  models  is 
that  the  creep  function  can  then  be  expressed  in  a  simple 
analytic  form,  as  a  sum  of  exponential  terms.  This  permits 
mathematical  analysis  of  viscoelastic  problems  to  be  carried 
much  further  than  would  otherwise  be  possible,  bringing  to 
light  significant  features  of  the  solution  that  are  charac¬ 
teristic  of  actual  materials. 
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Model  representation  permits  a  stress-strain  law  to  be 
expressed  as  a  linear  uirrerentlal  cQUatlon^  rather  than  whs 
integral  equation  (1.7).  This  can  be  expressed 

P(a)  =  Q(6)  (1.10a) 


where  P  and  Q  are  linear  differential  operators  of  the 
form 


n 

P  =  2 

k=0 


n 

Pk  n? 


m 

Q  =  2 

k=0 


^k 


(1.10b) 


Pj^  and  being  constants  for  the  particular  model.  In 

three-dimensional  problems,  two  sets  of  operators  are  needed. 


For  shear  (l.lla) 

and  for  dilatation  P’(®)  =  Q'*(e)  (1.11b) 


More  details  and  further  references  on  viscoelastic  models 
and  their  differential  equations  are  given  in  Lee  [20], 

Bland  [3],  and  Kelly  [l8]. 

A  convenient  model  to  use  is  the  general  Voigt  material. 

This- is  made  up  of  a  series  of  Kelvin  elements  (spring  and 

dashpot  in  parallel),  with  a  separate  spring  to  represent 

initial  elasticity  if  it  is  present  (see  Fig.  1.3).  A 

model  representing  a  viscoelastic  fluid  must  also  have  a 

separate  dashpot,  giving  unlimited  creep  under  load.  For 

each  Kelvin  element,  the  stress-strain  law  is  of  the  form^ 

}  A  dot  indicates  the  derivative  with  respect  to  time; 
f(t)=  df/^t  . 
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a  =  E(e  +  T^),  where  T  =  t\/E  Is  the  retardation  time  of 

the  element.  A  Voi^t  material  thus  has  n  characteristic 

retardation  times,  one  for  each  Kelvin  element,  with  an 

initial  elastic  modulus  and  a  final  elastic  modulus 

n 

E-  =  E^  +  2  E,,  .  The  simplest  model  that  exhibits  the 

°  k=i 

initial  response,  delayed  creep,  and  final  elasticity  char¬ 
acteristic  of  a  viscoelastic  solid  is  a  Voigt  solid  with  a 
single  retardation  time,  called  the  standard  linear  (or 
three  parameter)  solid  (see  Pig.  1.3)* 

General  Voigt  Material  (fluid) 

(A  Voigt  solid  has  no  free  dashpot,  i.e . ,  -  »  . ) 


PIG-  1.3 


16 


More  Kelvin  elements  with  different  retardation  times 
can  be  added  to  give  a  more  complex  model  that  may  in  certain 
cases  satisfactorily  approximate  a  real  material.  However, 
real  materials  appear  to  have  no  discrete  retardation  times, 
so  a  model  with  a  large  number  of  elements  is  usually  neces¬ 
sary  to  accurately  represent  an  actual  material  over  a  long 
time  span.  The  detemination  of  constants  to  fit  such  a 
model  to  a  test  curve  can  be  very  tedious,  and  the  result 
is  still  only  an  approximation.  It  is  frequently  more  de¬ 
sirable  to  use  the  test  curve  directly  in  an  integral  law 
such  as  (1.9),  which  can  be  evaluated  numerically.  The 
approximation  is  then  numerical  rather  than  physical,  and 
the  resulting  error  can  be  estimated  and  controlled  more 
readily. 

The  creep  function  for  a  model  material  can  be  determined 

directly  from  the  differential  law  (1.10)  by  using  the 

Laplace  transform^.  The  linear  differential  operators  P 

and  Q  become  polynomials  in  s  ,  so  that  (1.10)  transforms- 

to  _  _ 

P(s)  a  =  q(s)  €  (1.12) 

where 

n  m  , 

P(3)  =  2  p,  s^  ,  q(s)  =  2  q,  s^' 

k=0  k=0  ^ 

T - 

The  Laplace  transform  of  function  f(t)  is: 

go  A 

f(s)  =^{f(t))  =  /  e  ^  f(t)dt  .  It  is  taken  in  the  general- 

0“ 

ized  sense,  starting  at  t  =  0“  ,  to  include  the  effects  of 
discontinviities  at  t  =  0. 
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Prom  the  definition  of  a  creep  function,  €(t)  =  J*(t)  when 
a(t)  =  H(t)  ,  thus  i(s)  =  J*(s)  when  a(s)  =  1/s  .  Then 
from  (1.12), 

J*(s)  =  i|  (1-13) 

For  a  general  "Voigt  material  (Fig.  1.3) >  the  result  is 


n 

2 

k=l 


1 


n 

2 

k=l 


1  +  sT* 


(1.14) 


th 

where  Tj  =  is  the  retardation  time  of  the  k  K.^lvin 

element,  f^^’  =  ,  and  T*  =  .  This  model  is  a 

fluid  material,,  because  of  the  presence  of  the  free  dashpot 
with  viscosity  .  In  a  solid,  there  is  no  free  dashpot, 
and  this  terra  in  (l.l4)  is  absent  (ii^  =  »,  T^  =  ») .  From 
(1.14),  the  creep  function  is 

J*(t)='^fl  +  —  +  2  fjl  -  (1.15a) 

%  '■  T*  k=l  ■' 

This  has  the  initial  response  J*(0‘^*)  =  ,  and  the  long- 

time  unlimited  creep  J*(t)  ~  ss  t  -►  «.  With  the 

Y]  term  absent,  the  solid  has  a  final  limiting  value  of 
°  n 

creep  J*(«')  =  1/^^  =(1/Eq)(1  +  2  f^^)  • 
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In  dimensionless  form. 


j(t)  =  [1  +  ^  +  S  fjl  -  e"'/'^k)]H(T)  (1.15b) 

o  1 

where  t  -  t/T  ,  T.  =  Tj/T  and  T  is  a  characteristic  time 

n 

parameter.  Then,  J(0+)  =  1  ,  and  J(oo)  =  F  =  1  +  2  fi^  . 

k=l  ^ 

For  a  standard  linear  solid,  n  =  1,  T  =  T*  .  Then, 

J(t)  =  1  +  f(l  -  e"*^)  ;  J(oo}  =  F  =  1  +  f 

Most  experimental  data  on  actual  viscoelastic  materials, 
such  as  high  polymers,  provide  information  on  the  behavior 
in  shear  or  simple  tension.  Very  few  tests  have  measured 
dilatation  alone,  so  data  for  the  creep  curve  B(t)  are 
frequently  not  available.  In  the  absence  of  specific  in¬ 
formation,  several  assumptions  of  dilatational  behavior  have 
been  made.  Each  may,  in  certain  circumstances,  be  a  good 
approximation  to  actual  behavior^.  The  common  assumptions 
are: 

1)  elastic  dilatation,  B(t)  -■  H(t) 

2)  incompressible,  =  «,  B(t)  =  0 

3)  similar  behavior  in  shear  and  dilatation,  B(t)  =  J(t) 

In  terms  of  Poisson’s  ratio,  v  ,  which  is  discussed  in  the 

next  section,  2)  corresponds  to  v  =  1/2,  and  3)  corresponds 
1 

For  further  discussion,  see  Staverman  and  Schwarzl  [29], 
Section  6d;  Ferry  [9],  Chapter  x3. 
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to  V  as  =  constant.  Even  when  actual  data  are  available, 
o 

use  of  one  of  these  assumptions  may  lead  to  significant 
simplifications  in  the  analysis  of  a  problem,  without  being 
too  seriously  in  error. 

1.3  Solution  of  Viscoelastic  Moving  Load  Problems 

Many  problems  in  viscoelasticity  can  be  solved  by 
removing  the  time  dependence  with  a  Laplace  transform^, 
which  reduces  the  viscoelastic  problem  to  an  "associated 
elastic  problem"  (Lee  [19,  2o]).  Taking  a  Laplace  transform 
of  the  viscoelastic  stress-strain  laws  ^1.9a,b)  gives 

These  are  the  same  as  the  transforms  of  the  elastic  laws 
(1.9c)  if  iIq  is  replaced  by  M-q/sJ  and  by  K^/sB  . 

If  a. Laplace  transform  of  the  boundary  conditions  is  pos¬ 
sible,  the  viscoelastic  problem  becomes  an  elastic  problem 
in  the  transformed  variables,  with  the  elastic  constants 
now  functions  of  the  transform  parameter  s.  This  is  the 
associated  elastic  problem,  and  if  its  solution  can  be 
found,  the  viscoelastic  solution  is  obtained  by  inversion 
of  the  Laplace  transform. 

^  The  Laplace  transform  is  taken  with  respect  to  the 
dimensionless  time  t  =  t/T  .  The  generalized  transform, 
starting  at  t  =  0“,  wi”  always  be  used.  In  this  way, 
initial  conditions  are  ail  zero  for  an  initi'^lly  undisturbed 
material. 
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For  an  associated  elastic  problem  to  exist,  it  must  be 
possible  to  transform  the  boundary  conditions.  There  are 
two  types  of  boundary  conditions:  tractions  prescribed  on 
a  region  R-j^  of  the  boundary  surface,  and  displacements  pre¬ 
scribed  on  a  region  Rg  of  the  surface.  If  the  regions  R^ 
and  Rg  do  not  change  with  time  (the  quantities  prescribed 
on  them  may  vary  with  time,  however),  then  a  Laplace  trans¬ 
form  is  possible.  If  R^  or  Rg  change  with  time,  a  Laplace 
transform  of  the  boundary  conditions  is  not  generally  pos¬ 
sible,  and  the  above  method  of  solution  is  not  applicable. 

A  procedure  equivalent  to  solving  the  associated  elas¬ 
tic  problem  (when  it  exists)  is  to  solve  the  original  vis¬ 
coelastic  problem  as  an  elastic  problem  (with  the  same 
boundary  conditions  as  the  viscoelastic  problem),  take  the 
Laplace  transform  of  this  solution,  and  replace  the  elastic 
constants  by  their  appropriate  viscoelastic  analogues  from 
(l.l6).  This  is  often  much  easier  than  solving  the  asso¬ 
ciated  elastic  problem  directly;  use  can  be  made  of  exist¬ 
ing  solutions  from  the  theory  of  elasticity.  Frequently, 
the  elastic  solution  is  in  the  form  of  functions  of  space 
and  time  multiplied  by  rational  functions  of  the  elastic 
constants.  The  transformed  viscoelastic  solution  is  then 
a  product  of  the  transformed  elastic  solution  and  creep 
functions,  and  the  inversion  is  Just  a  convolution  integral 
In  replacing  the  elastic  constants  according  to  (l.l6),  the 
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term  sJ  (or  sB)  occurs.  This  Is  the  transform  of 
j(T)  =  dJ(T)/dT  [or  B(t)  =  dB(T)/dT].  The  viscoelastic 
solution  is  thus  a  convolution  of  appropriate  creep  function 
derivatives  with  the  elastic  solution. 

As  an  example,  suppose  an  elastic  solution  has  the  form 

g®(x,y,z.T)  =  i  f(x,y,z,T) 


Taking  Laplace  transforms  gives 

g  ®(x,y,z,s)  =  ^  f(x,y,z,s) 

Replacing  |i  by  M-^/sJ  gives  the  solution  of  the  associated 
elastic  problem,  and  inverting  gives  the  viscoelastic  solution 

T 


g(x,y,z,T) 


J(C)  f(x,y,z,T-C)dC 


Elastic  solutions  frequently  contain  Poisson's  ratio, 

V  =  (3K  -  2m-)/2(3K  +  |x)  ,  which  depends  on  both  shear  and 
dilatation  properties  of  the  material.  In  forming  an  asso¬ 
ciated  elastic  problem,  v  is  replaced  by  this  same  function 
of  the  transformed  creep  functions: 


V 


2m..  3K  ul  3K.sj  -  2m.„sB 

+  4)  =  — r - 2__ 

sJ  sB  sJ  2(3K^^sJ  +  M-q  sB) 


(1.17a) 


This  can  be  considered  the  transform  of  a  viscoelastic 
Poisson ‘3  ratio,  v(t),  which  varies  with  time,  so  that 
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(1.17b) 


S  v{s)  5 


r.T  -  Pit  rR 

■'"Q--  -'"O - 

2(3KQsj  +  ^LqSB) 


Then,  v  is  replaced  by  sv(s)  in  the  associated  elastic 
problem. 

From  the  initial  and  final  value  theorems  of  Laplace 
transforms  (see,  for  example,  Thomson  [33]), 


3K  -  2H- 

v(0’^)  =  sv(s)  =  2(3K  +  u  ; 


(l.l8a) 


where  v  is  the  Poisson’s  ratio  of  the  initial  elastic 
o 

response,  and 


v(oo)  =  lim  sv(s)  =  v„ 

S"^  a> 


(1.18b) 


where  is  the  final  value  of  Poisson's  ratio.  If  the 
material  has  a  final  elastic  response,  then  its  Poisson's 
ratio  is  =  (3K^  -  2m.|.)/2(3K^  +  |x^) .  Using  (1.8c)  this 


becomes 


1  3K„  1 
"f  -  2  3K„  + 


Thus,  if  Fj  >  ,  and  if  P^  <  F^  . 

If  the  material  is  fluid  in  shear  and  solid  in  dilatation,  so 
that  J(oo)  =  00  while  B(oo)  is  finite,  then  v^.  =  1/2, 

An  example  of  a  simple  model  that  exhibits  features 
similar  to  more  general  viscoelastic  solids  is  one  with 
standard  linear  solid  behavior  in  shear,  and  elastic  dilatation. 
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1 


Then^. 


J(t)  =  [1  +  ]H(t)  }  bCt)  =  H('r) 

J(s)  =  sj(s)  =  1  +  j  B(s)  =  sB(s)  =  1 
3Kq[1  +  f/(l+s)]  - 


sv(  s) 


=  V.  1  + 


t9iTC  u.  jiv  T|x  T^A.  i  1 

^  ( 3K^ -2|i^ )  ( 3k^-Hi^+3K^ f }  3K^+I^^  'r)l]H(T) 


s  +  tl  +  3K^f/(3VnJT 
3K^+H^+3K^f 


■q  *^0 '  '  O  o 


O  '  o 


For  example,  if  =  (8/3)  ix^,  then  ^  1/3  , 

=  C3+4f)./(9+8f)  , 

v(t)  =  [l-exp(-  T)]]H(Tr) 

For  additional  examples  and  discussion  of  Poisson's  ratio, 
see  Kelly  [18],  Freudenthal  et  al.  [14], 


PIG,  1.4 


The  retardation  time  T*  is  taken  as  the  time  parameter  T. 
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For  the  present  problem  of  a  given  load  moving  on  the 
surface  of  a  semi-infinite  base,  boundary  conditions  on  the 
surface  z  =  0  are 


rr  vt  V  =  /■■<l(2<*“Vt,y)  Under  the  load 

)  *  Q  outside  the  load 


(1.19a) 


=  ty2(x*,y*,0)  0  everywhere  (I.19b) 


Since  the  tractions  at  each  point  on  the  surface  are  known 
for  all  time,  a  Laplace  transform  of  f  boundary  conditions 
is  possible,  and  the  problem  can  be  solved  as  an  associated 
elastic  problem.  The  vertical  surface  displacement  u^ 

(on  z  =  0)  will  be  considered  in  detail.  The  same  procedure 
can  be  used  with  the  known  elastic  solution  to  find  the 
viscoelastic  displacements  and  stresses  at  all  points  in  the 
semi-infinite  material. 

The  elastic  solution,  as  given  in  (1.2),  (I.3),  is 


u^(x*-  VTT,y*,0) 
P(x*-  VTTr,y*,0) 


1-v 

M- 


P(x*-  VTT,y*,0) 


(1.20) 


q (x',y‘ )dx‘dy ' 


VTT:-x')^+(y*-y')2 


The  superscript  e  indicates  the  elastic  solution.  Taking 
a  Laplace  transform  gives 

Ug  (x*,y*,s)  =  ^  P(x*,y*,s) 
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■ 


_ 


/  -  ^  \ 


ncpxa'jxtig  i/Xit:  fixasuic;  c;Qrxstf4nt.s  accoraxrig  xo  i^x.ifa; 

gives  the  associated  elastic  solution 


a  ^  sj(s)[l  -  sv(s)]  P(x*,y*,  s) 
’  ^o 


(1.21) 


Inversion  then  gives  the  desired  viscoelastic  surface 
displacement . 

Although  (1.21)  could  be  formally  inverted  as  it  stands, 
in  the  form  of  a  convolution  integral,  the  final  result  is 
seen  more  clearly  by  proceeding  further.  Let 


J(t)  =  [1  + 

B(t)  =  [1  +  B3^(t)}H(t) 

with 

J^{0)  =  B^(0)  =  0 

Then 

sJ(s)  =  14-  3j-j^(s)  ;  sB(s}  =  1  +  sB^(s; 


(1.22a) 


(1.22b) 


Also,  let  K  =  Then  (1.18a)  becomes 


Vq  =  (3k-2)/2(3k+1) 


(1^23) 


From  (1.17b),  after  simplification. 


1  -  sv(s)  =  (1-Vq) 


1  - 


9k 


sJt  -  sBt 
i  i 


3k+4  3K^'1  +  Sksu^+sB] 


(1.24) 
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Then,  using  ( 1 . 22d )  and  ( i . 24 ; , 


—  sJ( s) [ l-sv( s) ] 

kJL 

o 


1-v^ 

- 2.  1-1-  sj 

u  1 

o 


sj,  -  sB, 

- r~~r(  i+s J^) 


Now  let  7(t)  be  a  new  function  such  that 


sJt  -  sB. 


7(s)  -  sj^  -  3k+1  +  3k6J'i+sS^‘  (1-25) 


Then  (1.21)  can  be  written 


\  P(x*,y*,s) 


Inverting  this  gives  the  desired  surface  displacement 
l-v„  l-v„  J 

UjCxfyfT)  =  P(x»-vn:,y*,0)  +  J  7(C)p(x*-VTT+VTC.yfO)dC 

^  ^  O 

T 

=  u®(x*-VrT,y*)  +  J  7(0  u®(x*-VTt  +  VTC,  y*)dC 


where  is  the  displacement  due  to  the  Initial  elastic 

response,  given  by  (1.2,3)  for  a  three-dimensional  problem, 
and  by  (l.5b)  for  a  two-dimensional  (plane  strain)  problem. 
Expressed  In  coordinates  moving  with  the  load 

(x  =  X*  -  VTt,  y  =  y*), 

T 

u^Cx.y.t)  =  u®(x,y)  +  J 7(C)  u®{x  +  TCC,y)dC  (1.26) 


The  explicit  evaluation  of  requires  the  function 

7(t).  If  the  creep  functions  J(t)  and  B(t)  can  be  ex¬ 
pressed  analytically,  so  that  their  Laplace  transforms  can 
be  found,  then  direc'-  inversion  of  (1.25)  gives  7(t).  How¬ 
ever,  for  actual  materials,  the  creep  functions  cannot 
usually  be  expressed  In  simple  analytic  form,  but  are  avail¬ 
able  only  as  curves  or  tabulated  values  from  tests.  Instead 
of  the  several  approximations  and  considera'-^le  labor  that 
may  be  needed  to  fit  a  curve  with  an  analytic  function,  take 
its  transform,  then  invert  (l.25)>  it  is  better  to  obtain 
7(t)  directly  from  the  given  data  by  a  numerical  procedure. 

For  this  purpose,  (1.25)  can  be  written: 

r  3i£.sJ-|-i-  sS-.  1 

=  "3k  "+  i  '  j  =  +  SPT  (3ksJ;^+  sB^)  sJj 


Inversion  of  this  gives 


3Kj;^(0+Bj(C) 

^55+1 


7(T-C)dC  =  jj^(T) 


3kJ,(0+'i6,(0  . 

- -3K+it  - 

(1.27) 
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This  is  a  Volterra  integral  equation  of  the  second  kind  for 
the  unknown  function  tCt).  It  can  be  solved  numerically  by 
a  finite  difference  method,  requiring  numerical  values  of 
the  derivative  of  each  creep  function  at  intervals  from  0 
to  T.  A  procedure  for  solving  such  integral  equations  has 
been  developed  and  discussed  by  Lee  and  Rogers  [21] .  Using 
a  slight  modification  of  their  procedure,  equation  (1.27) 
can  be  solved  numerically  as  follov/s:  Let 


P(t)  =  j^(T) 


M{C)  “  3^  +  %(01 

,•  -s  r  3Ki,(C)+'IBi(C) 
'  (3k+4)(3k+1)^'^1‘®1^  ,/  SkW 


j\(T-C)dC 


lyi(^)  and  PCt)  are  known  or  can  be  found  from  the  given  creep 
functions  J('r)  and  B(t).  Then,  (1.27)  becomes 

T-C)7(C)dC  =  F(t) 
o 

Divide  the  time  into  intervals  ,  i  =  1,2, n+lj 
■'l  =  ■'n+l  =  ^n+1  = 

=  7(0)  =  F(0)  »  i,{0)  -  [j,(0)-B,{0)l 


■^n+1  1  +  1/2  r  (■'n+l)  " 

j.n  =  1,2, 

29 

In  this  way,  7(t)  can  be  evaluated  in  successive  steps, 
starting  from  t=0.  The  time  intervals  can  be  varied  as 
desired j  equal  intervals  in  log  i  may  be  convenient. 

The  function  7(t)  will  be  qualitatively  similar  to 
J^(t),  as  can  be  seen  from  (1.27).  It  is  a  monotonic  de¬ 
creasing  function,  concave  upward.  If  the  material  is  a 
solid,  the  creep  functions  approach  asymptotically  a  limit¬ 
ing  value]  thus,  j^(oo)  =  =  0  ,  so  7(00)  =  0  .  If 

the  material  is  a  fluid,  the  creep  functions  approach  a 


linear  function  of  time  after  a  long  time,  so  y{<y>)  =  >  0. 

For  a  general  Voigt  model,  from  (1.15b), 


(t)  = 


2  e 
1  -^k 


If  the  material  is  a  solid,  T^  •=  »  ,  and  J2^(t)  decreases 
exponentially  to  zero.  .An  actual  material  will  have  the 
same  qualitative  behavior  for  J-j^(t)  ,  and  thus  for  7(t). 
A  representative  sketch  is  shown  in  Pig.  1.5* 
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If  the  material  is  assumed  to  have  constant  Poisson's  ratio 
[B(t)  =  or  incompressible  (v  =  1/2)],  then  7('r)  = 

j^(T)  .  If  the  dilatation  is  assumed  elastic ,  =  0  in 

(1.27).  As  an  example,  for  a  material  with  standard  linear 
solid  behavior  in  shear  and  elastic  dilatation. 


J(t)  =  [1  +  f(l-e‘'')]  H(t)  ;  B(t)  =  H(t)  (1.28a) 


j^(T)  =  f  e'"'  ,  83^(1)  =  0 


Inversion  of  (1.25)  gives 


=  211^  IsIfTF 

°  (1.28b) 


where  = 


_  3k-2 


o  "  2(3k+1}  " 


or  K  = 


1+v 

2  o 


3  l-2v. 


For  any  material  whose  shear  and  dilatation  behavior  are 
represented  by  mechanical  models,  such  as  the  general  Voigt 
solid,  7(t)  will  be  a  sum  of  exponentials,  in  the  form 

7(t)  =2  . 

1  ^ 

tor  further  applications,  it  is  desirable  to  express 
the  surface  displacement  (1.26)  in  dimensionless  fom.  .This 
is  most  conveniently  done  by  using  as  a  characteristic  length 
parameter  the  quantity  VT  .  Then,  the  following  dimension¬ 
less  variables  are  introduced: 
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Coordinates : 

(moving  with  the  load) 

fi  .•=  2L.  n  =  ^ 

-  yip  >  '  -  VT 

(1.29a) 

u_(x,y,0) 

(1.29b) 

Displacement : 

v(4>n)  = 

1-v  , 

Pressure : 

Q(4>n)  ^  Tj.  q(x,y) 

(1.290) 

Using  these,  the  elastic  displacement  is,  from  (1.20), 

=  i  ff  (1.30) 

y  ’/(e-l')®  +  (n-n')® 


The  viscoelastic  result  (1.26)  then  becom.es: 

T 


Hi.n.-c)  =  v®(5,t,)  +  J  -y(c)  v®(|+c,n)dc  (1.31) 

O 


For  a  general  three-dimensional  problem,  using  (1.30), 


-A-  {n-r,') 

T 


I 


+ 


I 


7(C)c3C 


de’dii'  (1.32) 


o  ^(e-|Vc)^+(Ti-Ti')^  ' 


For  a  plane  strain  problem,  from  (l.5a)> 

A 


ieo-5'l 


'(5<t)  -  v{5g,T)  =  J  Q(c')  Jlog 


/7(C) 


log 


dC 


(1^33) 
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T 


0 

where  \  ^  is  the  dimensionless  length  of  the  loaded 

region. 


1,4  Steady  State  Solution 


The  stress  and  deformation  patterns  in  the  base  change 
with  time  due  to  the  delayed  creep  and  recovery  of  visco¬ 
elastic  materials.  After  a  sufficiently  long  time,  however, 
the  transient  effects  may  disappear,  so  that  these  patterns 
are  unchanged  with  further  passage  of  time.  In  this  case, 
a  steady  state  exists  for  an  observer  moving  with  the  load. 
With  respect  to  coordinates  fixed  in  the  base,  functional 
dependence  on  time  then  occurs  only  in  the  combination 
X*  “  Vt  .  In  coordinates  x,y,z  moving  with  the  load,  there 
is  no  explicit  time  dependence.  The  question  of  whether  and 
under  what  conditions  a  steady  state  occurs,  and  some  of  its 
consequences,  will  be  considered  in  this  section. 

The  viscoelastic  solution  (I.3I)  contains  time  explioity 
only  as  the  upper  limit  of  the  integral.  If  this  integral 
exists  as  t  -►  «  ,  a  steady  state  is  reached,  and  is  given 
by  (1.31)  wit>'  T  =  00  .  If  the  value  of  the  integral  in¬ 
creases  without  limit  as  t  ,  no  steady  state  exists. 
Consider  the  integral  in  equation  (I.31): 

r  » 

I(t)  =  J  7(C)  .’®(|+C.'))dC 


o  A 


(1.34) 
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where 


If  l(«»)  exists,  there  is  a  steady  state.  Let  be  a 

large  but  finite  time.  Then, 


l(oo)  = 


00 

r 


J 


7(0  JJ  Ti' 


(1.35) 


I(Tq)  is  finite,  so  only  the  second  term  must  be  considered. 
Since  ^  is  always  large  in  this  integral,  the  integrand 
can  be  expanded  in  terms  of 


f  -fyU)[ 

*1'^  A 


Jdc 

(1.36a) 


where  N^,  ....  are  finite  terms  depending  on  (I,!))  but 

not  C.  If  the  material  is  a  viscoelastic  solid,  7(C)  de¬ 
creases  monotonically  to  zero  as  ^  .  In  this  case, 

the  Integral  (1.36a)  is  finite.  Thus,  !(«>)  exists;  a 
steady  state  is  always  reached  when  the  semi-infinite  mate¬ 
rial  is  a  viscoelastic  solid.  If  the  material  is  fluid, 
7(0  decreases  monotonically  to  a  constant  value,  7^  , 
as  ^  00.  Then,  7(0  =  7^^  +  7]^(C)  ^  where  7^(0  decreases 

monotonically  to  zero  as  C  “*■  «  (see  Pig.  1.6).  The  inte¬ 
gral  (1.36a)  becomes 


+  ...]dC 
(1.36b) 
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The  second  Integral  is  finite,  but  the  first  term  integrates 
to  7^Q[log  ,  which  is  unbounded,  so  l(oo)  does  not 

exist.  Thus,  a  steady  state  is  not  reached  when  the  material 
is  a  viscoelastic  fluid. 

If  the  problem  is  two-dimensional, 

^  ^  I  ?  t‘C  *1 

i(t)  =  JyU)\J'i{i')  log  iirrTT 

o  o 


Following  the  same  procedure  as  before,  an  integral  of  the 
same  form  as  (I.36)  results.  Thus,  the  conclusions  are  the 
same  for  the  two-  and  three-dimensional  cases;  There  is  a 
steady  state  when  the  material  is  a  solid;  there  is  no 
steady  state  for  a  fluid. 

For  a  viscoelastic  fluid,  w’here  ^(C)  =  +  7^(C)> 

the  displacement  (I.3I)  becomes 

T  1 

v(e.n,T)  =  v®(e,ii)  +  7f  Jv®{4+c,Ti)dc  +  j  7i(c)v®(i+c,ii)dc 

o  o 

^  The  symbol  ~  means  "asymptotically  approaches";  l.e., 
f(x)'^g(x)  as  X a  means  lf(x)-g(x)|  =  0 

X”^oo  ' 
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The  term  with  is  due  to  the  unlimited  creep  of  the 

material,  and  this  is  what  fails  to  reach  a  steady  state. 
The  fluid  behaves  after  a  long  time  as  if  it  were  puxely 
viscous,  as  given  by  the  term  with  7^.  Thus,  to  decide 
if  a  steady  state  is  reached  for  a  general  viscoelastic 
fluid,  it  is  necessary  to  consider  only  an  ideal  viscous 
fluid.  A  viscoelastic  fluid  will  have  a  steady  state  if 
and  only  if  the  ideal  viscous  fluid  does. 

The  failure  to  reach  a  steady  state  when  the  base  is 
a  fluid  has  been  shown  here  for  a  semi -infinite  base.  In 
this  case  there  are  no  constraints  or  supports  to  prevent 
the  displacement  from  increasing  without  limit  as  the  fluid 
material  undergoes  unlimited  creep.  If  the  base  was  not 
semi-infinite,  but  a  layer  of  finite  thickness  supported 
from  below,  the  creep  would  be  constrained  and  a  steady 
state  might  result.  This  possibility  v;ill  be  considered 
in  more  detail  in  Chapter  II. 

When  the  material  is  a  viscoelastic  solid,  the  steady 
state  solution  is  given  by  (I.3I)  with  t  =  c»  : 


(1.37) 


The  steady  state  is  evident  only  to  an  observer  moving  with 
the  load  at  velocity  V;  the  coordinates  (^,11)  are  also 
moving  with  the  load.  It  is  convenient  to  think  instead 
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of  a  fixed  load  and  coordinate  system,  with  the  semi-infinite 
base  moving  uniformly  under  the  load,  in  the  negative  | 
direction  (Fig.  1.6).  A  particle  of  the  base  comes  from  the 
right  (I  =  oo)  and  has  been  traveling  horizontally  for  a  long 
time.  It  continues  to  move  with  constant  horizontal  velocity 
V  ,  but  it  also  acquires  a  vertical  velocity  as  it  comes  under 
the  influence  of  the  load.  The  path  of  the  particle  is  de¬ 
formed  in  the  vicinity  of  the  load,  but  the  particle  is 
eventually  moving  horizontally  again  far  downstream.  Parti¬ 
cles  initially  on  the  surface  remain  on  the  surface,  and 
undergo  the  vertical  displacement  u  (x*-Vt,y*)  =  u  (x,y). 

The  vertical  velocity  of  surface  particles  is 


U*  ^ 


3u 

z 

SF" 


du 


=  -  V 


z 


Sx 


=  -  V 


av(i,T)) 

— 


(1.38) 


Thus,  the  vertical  velocity  of  a  particle  on  the  surface  is 

prooortional  to  the  slope  of  the  deformed  surface. 

% 


PIG.  1.6 
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when  a  steady  state  does  exist,  the  resulting  expressions 
jr  surface  displacement  (or  any  other  dlj.placement  or 
stresses)  can  be  found  directly,  without  the  necessity  of 
considering  time  explicitly.  Consider  the  load  and  coordinate 
system  (x,y,z)  fixed,  and  the  viscoelastic  base  irioving  at 
constant  velocity  V  in  the  negative  x  direction.  For  s.1mplic- 
ity,  a  one -dimensional  state  cf  stress  is  considered.  The 
argument  is  easily  extended  to  shear  components  and  dilatation 
in  a  general  three-dimensional  state. 

An  element  A  in  the  base,  now  at  coordinate  x,  moves^ 
on  the  pathline  A' A  (Fife.  1.6).  Suppose  while  at  x‘  it  re¬ 
ceived  a  stress  increment  da  =  [^a(x')/^x’]  dx’.  This  has 
influenced  the  element  for  a  time  (x'-x)/V,  and  the  resulting 
strain  increment  is 


d€(x)  =  J('"^-.^)  da 


Adding  up  all  such  increments  from  the  initial  position  of 
the  element  at  x*  =  oo  to  its  present  position  a  x’  =  x 
gives  the  total  strain 


€(x) 


x’-x,  6a(x') 
"VT  ' 


03 


dx 


(1.39) 


1 

Deformations  due  tc  the  applied  load  are  infinitesimal 
and  can  be  neglected,  so  the  orily  motion  is  in  t.ie  negative 
x-direction  at  velocity  V. 
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The  elements  are  coming  from  an  undisturbed  state,  so  that 
0(00)  =  0  .  Thus,  Integrating  (1.39)  by  parts,  and  letting 
C  =  (x'-x)/VT, 


e(x)  = 


.,.(x+VTO  j(C)dC 


(1.40) 


This  Is  the  steady  state  equivalent  of  the  stress-strain  law 
(1.7).  Similar  laws,  analogous  to  (1.9a,b),  for  shear  and 
dilatation  can  be  obtained  in  the  same  waj  ,,  The  stress-strain 
law  (l.40)  can  be  considered  an  operational  extension  of  the 
elastic  law  e(x)  =  l/E^  a(x)  ,  with  l/E^  replaced  by 
J(C)/Eq  and  the  integration  carried  out  as  indicated  in 
(1.40).  In  the  same  way,  (1-v  Vfi  [6(^)  +  7(C)]  replaces 
(1-Vq)/m,^  in  the  elastic  suri'Ce  displacement  (l.3)>  which 
leads  to  the  viscoelastic  steady  state  displacement^ 


U2(x,y) 


w 

U2(x,y)  =  u®(x,y)  ^  J  7(C)u®(x+VTC,y)df; 


In  dimensionless  terms; 


v(^yTj)  =  V^(4,T^)  +  J  7(C)v®(^+Cy1l)dC 


which  is  identical  to  the  previous  result  (I.37). 

This  operation  can  be  justified  more  rigorously  using 
Fourier  transfoms  in  a  way  similar  to  the  use  of  Laplace 
transforms  n  the  preceding  section. 
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The  steady  state  solution  for  the  surface  displacement 
of  a  semi-infinite  base  is  given  by  (1.32)  or  (1.33),  with 
T  =  00.  The  result  for  the  three-dimensional  problem  can  be 
written 

v(l,Ti)  =  /TQ(?Vn')K*(5-?'  =  n-’i')ae'<an'  (i.ina) 


where 


00 

1  _L_  .  1  r _ 'y(c)dc 


(i.4lb) 


For  a  two-dimensional  problem, 

A 

v{4)  -v{?o)  =  J  «{5')[K(eo-i')-K(e-e')]de'  (1.42a) 

o 

where 

K(0  =  login 

o 

These  results  are  for  a  viscoelastic  solid  with  initial 

elasticity.  The  dimensionless  terms  are  defined 

in  (1.29),  while  7(C)  is  obtained  from  the  creep  functions 

of  the  material  by  (1.25)  or  (1.27). 

It  is  of  interest  to  examine  the  surface  displacement 

far  away  from  the  load.  The  loaded  area  will  be  taken  to 

include  the  origin  of  the  coordinate  system.  In  the  three- 

2  2  2 

dimensional  case,  let  p  =  ^  '  be  very  large;  p  is 

the  (dimensionless)  distance  of  a  point  on  the  surface  from 
the  origin.  Then,  as  p  -*  <»  , 


00 


J  7(0  log|^+C|  dC 


(1.42b) 


40 


~  [1  +  [  -^COdC] 

o 


Thus,  when  +  r\^ 


— 2  i  J CifOdC  +  ••. 

2p  ^ 


v(^.Tl) 


00 


ri')d|'dr,'  [1  + 


Using  (1.23)  and  the  final  value  theorem  of  Laplace  trans¬ 
forms,  it  can  be  shown  that 


1  +  /\(C)d4  =  [1  +  y{s)]  = 


1-V. 


J 

o 


^f 


( 1. ^3) 


Then, 


v(^.n) 


1-v 


The  total  load,  N*,  applied  to  the  surface  is 


N*  =  ^q(x*,y’)dx‘dy’  =  ( VT)^  JJq(  ^  ‘ H  ’  )d| ‘dr]  ‘ 
A  X 


Thus,  far  away  from  the  load, 

^  N* 


v(^,Ti) 


1-v, 


2ir{VJ})'^  ^f 


,  or 


U  (x,y) - 1 _  in 

Sir  IXf 


(1.44) 


This  is  the  same  displacement  that  would  result  from  a 
concen-crated  load  of  magnitude  N*  applied  at  the  origin. 
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The  result  is  that  far  from  a  localized  pressure  distribution 
mo’  on  the  surface  of  a  viscoelastic  solid,  the  surface 
di  .placement  is  the  same  as  that  for  a  concentrated  load  of 
the  same  total  magnitude,  acting  at  the  origin  on  the  surface 
of  an  elastic  solid  with  the  final  elastic  moduli  of  the 
viscoelastic  material.  This  is  true  in  any  direction,  whether 
ahead  of  or  behind  the  moving  load.  The  initially  flat  sur¬ 
face  is  deformed  in  the  vicinity  of  the  load,  but  gradually 
returns  to  its  initial  level  after  the  load  passes.  Far 
enough  away  from  the  load,  the  actual  pressure  distribution 
has  no  effect;  only  the  total  load  is  significant. 

In  a  two-dimensional  problem,  as  |^|  -►  «  , 

v(C)  -  ~  11^  Uogl^^l  -  log|e|]  (1.^5) 

Here  the  result  is  the  same  as  for  a  line  load  of  the  same 

total  magnitude,  N*  =  /  q(x‘)dx’  ,  per  unit  length  in  the 

o 

y-dlrectlon.  Again  the  base  exhibits  the  final  elastic 
behavior  of  the  viscoelastic  material.  The  surface  levels 
far  ahead  of  and  behind  the  load  are  the  same,  and  the  slopes 
tend  to  zero  as  }^j  -►  «  ,  For  any  two-dimensional  (plane 
strain)  problem,  elastic  or  viscoelastic,  the  load  is  dis¬ 
tributed  on  an  area  unbounded  in  the  y-direction,  and  thus 
the  total  load  is  infinite.  There  is  no  lower  boundary  to 
the  base  to  restrict:  the  resulting  deformation,  so  the  dis- 
placeme^^t  far  from  the  load  is  infinite  relative  to  points 
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i 

i 


I 


near  the  load.  The  displacement  rn'ii-it  thus  be  expressed 
relative  to  an  arbitrary  point  near  the  load. 

For  the  extremes  of  very  large  or  very  small  velocity, 
the  viscoelastic  solution  approaches  as  a  limiting  case  the 
corresponding  elastic  solution.  For  V  =  0  ^  the  viscoelastic 
material  exhibits  its  final  elastic  response  when  a  steady 
state  is  reached,  and  the  result  is  the  same  as  if  it  were 
an  elastic  solid  with  the  final  elastic  behavior.  For  very 
large  V  ,  the  load  passes  so  fast  that  there  is  time  only 
for  the  initial  elastic  response.  Thus,  as  V  ,  the 
result  is  the  same  as  if  it  were  an  elastic  solid  with  the 
initial  elastic  behavior^.  These  results  are  confirmed  from 
the  limiting  values  of  the  surface  displacement.  Written 
out  completely,  the  steady  state  solution  is,  from  (l.4la,b) 


1-v. 


1_ 

2'w 


A  Mx-x«)  +(y-y‘) 


+  ■  . . -  j  dx  ’  dy  ’  (1,46) 

o  -X ’)^+  ( y -y* 


For  7=0,  (1.46)  becomes 


1-v  ”  1 

U2(x,y)  =  [1  +  J"  y(OdC]  P(x,y,0)  =  P(x,y,0) 

*^0  •  ^ 


The  solution  in  this  case  is  no  longer  physically  mean¬ 
ingful,  because  the  problem  cannot  be  treated  as  quasi-static. 
The  limiting  case  may  be  approached,  however,  for  relatively 
large  velocities  still  within  the  restrictions  of  the  quasi¬ 
static  assumption. 
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This  Is  the  elastic  solution  (1.3)^  with  elastic  constants 

final  elastic  response.  For  V ^  the 
inner  integral  in  (1.46)  vanishes,  so 

1-v 

U2(x,y)  =  P(x,y,0) 

v^hich  is  the  elastic  solution  for  the  initial  elastic 
response. 


1.5  Energy  Dissipation 

When  a  perfectly  elastic  material  is  loaded,  work  done 
on  the  body  as  it  deforms  is  stored  internally  as  strain 
energy.  If  the  body  is  subsequently  unloaded,  recovery  takes 
place  instantaneously.  The  strain  energy  is  released  and 
there  is  no  net  work  done  in  a  complete  cycle.  In  a  visco¬ 
elastic  body,  however,  recovery  takes  place  gradually  after 
unloading,  and  not  all  of  the  stored  energy  is  released. 

Energy  is  dissipated  within  the  material  during  the  loading 
and  unloading  process.  Thus,  net  work  must  be  done  on  a 
viscoelastic  body,  even  though  it  may  eventually  return  to 
its  Initial  shape.  The  amount  of  work  depends  not  only  on  the 
initial  and  final  states,  but  on  the  history  and  rate  of  load¬ 
ing  during  the  process.  In  a  mechanical  model  of  a  visco¬ 
elastic  material,  the  dashpot  elements  account  for  internal 
energy  dissipation,  v/hile  the  springs  provide  for  storage  of 
free  energy.  The  dissipation  in  actual  materials  is  often 

attributed  to  "internal  friction"  or  "elastic  hysteresis  losses." 
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Consider  a  moving  load  applied  to  some  body  in  contact 
with  the  surface  of  base  material.  The  contact  is  perfectly 
lubricated  (or  "frictionless")  so  there  are  no  tangential 
surface  tractions.  The  body  would  move  at  constant  velocity 
with  no  resistance  if  the  base  were  perfectly  elastic.  How¬ 
ever,  resistance  is  encountered  with  actual  base  materials, 
and  a  force  is  needed  to  maintain  constant  velocity.  A 
common  example  is  the  resistance  to  rolling  of  a  cylinder 
or  sphere.  Experimental  evidence^  indicates  that  rolling 
resistance  is  nearly  independent  of  the  roller  material  and 
of  the  lubrication  between  roller  and  base.  Also,  when  the 
surface  is  well  lubricated,  resistance  is  nearly  the  same 
for  rolling  and  sliding.  This  indicates  the  resistance  is 
due  to  bulk  properties  of  the  base  material  and  not  to  sur¬ 
face  phenomena.  Energy  dissipation  in  the  base  material  as 
it  is  deformed  is  the  primary  factor.  Work  must  be  done  by 
the  moving  body  to  provide  the  energy  dissipated  in  the  base, 
and  this  accounts  for  the  resistance  to  motion.  In  the  al¬ 
ternate  view  adopted  here,  the  base  moves  under  a  fixed  load, 
and  external  work  must  be  done  on  the  base. 

Since  a  viscoelastic  material  dissipates  energy,  the 
problem  of  a  moving  load  on  a  viscoelastic  base  being  con¬ 
sidered  in  this  chapter  is  one  way  to  account  for  and  pre¬ 
dict  resistance  to  moving  loads.  The  initially  flat  surface 

^  Bowden  and  Tabor  [4],  Tabor  [30,31 ]>  Tabor  and  Atack 
[32],  Flom  [10,11]. 
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of  the  baL^e  is  deformed  by  the  load  q(x,y).  Because  there 
are  no  tangential  stresses  on  the  surface,  the  load  must 
remain  normal  to  the  deformed  surface.  Consistent  with 
other  assumptions  of  the  linear  theory,  the  effect  of  the 
infinitesimal  horizontal  components  of  the  load  can  be  neg¬ 
lected  compared  to  the  vertical  components.  Thus,  in 
solving  the  problem  for  stress  and  displacements,  the  load 
is  considered  to  remain  vertical  and  applied  to  the  unde¬ 
formed  surface.  To  evaluate  the  energy  dissipation,  however, 
the  horizontal  component  of  the  load  and  deformation  of  the 
surface  must  be  considered. 

At  each  point  on  the  surface  under  the  load  there  is  a 
vertical  component  =  q  ccs0  and  a  horizontal  component 
qj^  =  q  sin0,  where  tan0  =  -  du^/^x  is  the  slope  of  the 
deformed  surface  (Fig.  1.7) •  With  the  restriction  to  small 


PIG.  1.7 
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slopes,  the  total  vertical  load  Is 


N' 


=  JJ  q(x,y)  dx  dy 


The  resultant  horizontal  force  at  the  contact  is 


F  =  -  Jj  q(3c,y)  - -  dx  dy 

A 


(1.^7) 


It  is  convenient  to  introduce  dimensionless  forces  N  and  F, 
given  in  terms  of  the  quantities  of  (1.29)  hy 

(1.48) 

^  -  -  Jf  ^U,n)  or,  (1.49) 

In  a  two-dimensional  problem,  the  body  is  an  infinite  cylin¬ 
der  (not  necessarily  circular),  with  axis  in  the  y  direction. 
Then,  N  and  F  are  forces  per  unit  length,  given  by 

1— V  ^ 

N  =  ?  W  =  / (4.50) 

o 

^  ?  W=  (4.51) 

0 
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The  force  F"*  is  entirely  a  consequence  of  the  visco¬ 
elastic  behavior  of  the  base.  The  deformation  will  not  be 
symmetric,  even  for  a  symmetric  load  distribution,  because 
the  viscoelastic  material  experiences  delayed  recovery  as  it 
moves  under  the  load.  There  will  always  be  a  horizontal  re¬ 
sultant  of  the  contact  pressure,  acting  on  the  moving  base 
(or  body)  opposite  to  the  direction  of  motion.  This  will  be 
shown  in  the  two-dimensional  case  for  simplicity  (the  three- 
dimensional  case  is  similar).  Prom  (l.42a,b), 

o  o 

Substituting  this  in  (1.51)  gives; 

70^  “Kh-  00 

^  =fj «(?)«(?')  +  ff  <i(i)«(e')  f 

00  'bo  o 


The  first  term  is  the  contribution  from  the  initial  elastic 
response ; 

XX 

^0=  Jf  «(«)«(!) 

00 

XX  XX 

Fo  -Jf  dUMV)  |4r|^  =  -JJ'iiiMi')  |4|4^=  -  Fo 

00  oo 


interchanging  ^  and  gives 
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Thus,  F  =  0  :  there  is  no  contribution  from  the  elastic 
o 

response.  Then, 

'h'h  00 

F  =  )«(«')  J  deal'  (1-52) 

oo  o 


The  horizontal  force  depends  directly  on  the  function  7(C) 
which  characterizes  the  viscoelastic  behavior  of  the  material. 
For  an  elastic  material,  7(0  =  0  ,  so  F  =  0. 

At  the  contact,  the  force  F*  acts  to  resist  the  motion, 
and  can  be  referred  to  as  a  ’’resistance"  or  ’’friction’’  force. 

A  convenient  dimensionless  measure  of  F*  is  the  coefficient 
of  friction  Xi  givsn  by 


X 


_  F* 
= 


F 


"pOUTnF 


(1.53) 


This  coefficient  varies  with  the  velocity,  total  load,  and 
load  distribution.  With  the  restriction  to  small  slopes, 
^u^/3x  =  dv/3^  «  1  ,  ana  thus  F*  «  N*  (x  «  1).  This 
agrees  with  observed  evidencei  for  example,  resistance  to 
rolling  on  hard  surfaces  is  usually  found  to  be  very  small. 

For  equilibrium,  there  must  be  resultant  horizot.tal 
forces  equal  to  F*  acting  on  the  body  and  on  the  remote 
boundaries  of  the  base  (Fig.  1.8).  These  forces  are  needed 
to  move  the  base  and  hold  the  body  .fixed  (or,  alternately, 
to  move  the  body  and  hold  the  base  fixed).  Work  is  done 
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on  the  base  at  Its  remote  boundaries,  at  a  rate 

W  =  VP*  =  -  V  JJ"  ci{x,y){hu^/hx)  dxdy  (l‘5^) 

This  equals  the  rate  of  energy  dissipation  within  the  visco¬ 
elastic  base.  No  work  is  done  at  the  contact,  since  particles 
are  moving  under  the  load  at  right  angles  to  the  direction 
of  the  load. 

In  the  limiting  cases  of  very  fast  or  very  slow  velocity 
(V  -*>  0  or  V  -*•  «)  it  was  shown  in  Section  1.4  that  a  viscoelastic 
solid  behaves  like  an  elastic  solid,  and  the  viscoelastic  solu¬ 
tion  is  the  same  as  the  elastic  solution.  As  these  limits 
are  approached,  it  is  a  plausible  expectation  that  energy 
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dissipation  and  resistance  to  motion  will  decrease  to  zero, 
since  they  are  zero  in  the  elastic  problem.  This  is  now 
shown  to  be  the  case.  Prom  (1.52)> 

l_v  - 

=  ?  JJ  'l(x)q(x')  J  dxdx'  (1.55) 

°  oo  o 

For  V  =  0  ,  this  becomes 

-t  OLi  ££i 

=  i  /  .(Oac  .  JJ  ,(x)a(x') 

o  oo 

The  double  Integral  has  been  shown  previously  to  be  zero,  so 
P*  =  0  .  For  V  “  00  ,  the  Integrand  vanishes,  so  again 
F*  =  0  .  The  resistance  increases  from  zero  as  V  increases 
from  zero,  and  decreases  to  zero  as  V  gets  very  large.  It 
will  have  at  least  one  maximum  at  some  intermediate  velocity. 
For  a  given  load,  N*  is  fixed,  so  x  =  F*/^*  the  same 

behavior  as  F*.  The  resistance  to  a  given  load,  due  to 
energy  dissipation  as  the  velocity  varies  can  be  seen  best 
by  plotting  x  vs  Pj  Fig.  1.9  shows  how  this  might  appear. 
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There  is  an  alternate  view  of  the  way  in  which  work  Is 
done  on  the  moving  viscoelastic  base.  The  load  q;(x,y)  Is 
taken  as  strictly  vertical  on  the  deformed  surface.  Particles 
under  the  load  have  a  vertical  velocity  u*  =  -  V  du  /3x 
(1.38)  In  the  direction  of  the  load.  Thus,  a  load  element 
qdA  does  work  at  a  rsite- 

dW  =  U*  q(x,y)dA  =  -  V  q(x,y)dA 

The  total  rate  of  work  done  is  then 

3u 

W  =  -  V  //  (3^)  q(x,y)  dxdy 
A 

Ihls  is  (in  the  linear  theory)  the  same  result  as  before 
(1.5^).  But  here  the  work  is  done  at  the  contact  of  load 
and  base;  no  work  is  done  at  the  remote  boundaries.  Thus, 
the  source  of  energy  is  quite  different  from  before.  This 
alternate  view  gives  rise  to  difficulties  in  physical  inter¬ 
pretation.  The  vertical  pressure  is  not  normal  to  the  de¬ 
formed  surface,  so  there  are  tangential  stresses.  This  con¬ 
tradicts  the  assumption  of  a  perfectly  frictionless  contact, 
since  the  base  and  load  are  in  relative  motion. 

2.6  Examples  of  Moving  Loads 

The  results  given  so  far  in  this  chapter  apply  to  any 
load  distribution  and  any  linear  viscoelastic  material.  Some 
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examples  of  particular  loads  are  now  considered.  The  material 
used  will  he  a  general  Voigt  solid.  This  adequately  repre¬ 
sents  characteristics  of  actual  materials  such  as  polymers^, 
while  allowing  results  to  be  expressed  in  analytical  form. 

For  such  a  material, 

7(^)  =  1  f  e-V  (1.56) 

i=l 

The  steady  state  surface  displacement  (1.37)  is  then 

00 

v(l,n)  =  V®(|,n)  +2  f  v®(e+C,T,)dC  (1.57) 

i=l  ^  J 
o 

Integrating  by  parts  gives 

m  f.  ^  m  f .  „ 

vce,i)  =  (1  +  2  gi  )  v®(e,t))  +  2  gi  v''(e,Ti}b.) 

i=l  1=1  '^l  ^ 

(1.58a) 

where 

v^(^,T);b)  =  J  dC  (1.58b) 

o 

The  elastic  part  of  the  displacement  gives  no  contribution 
to  the  resisting  force.  Thus, 

F  =  -  2  gi  F''(b.)  (1.59a) 

i=l  ^1  ^ 

I - 

Ferry  [9]  and  Tobclsky  [34]  give  creep  or  relaxation 
data  for  many  polymers,  with  discussions  of  their  properties 
and  references  to  original  literature. 


where 

F''(b)  =  JJ  <4(4, T))  (1.59b) 

A 

The  results  (1.58)^  (1«59)  sums  of  terms  which  ure 
Identical  except  for  the  value  of  the  constants  f^  and  b^; 
there  is  one  term  for  each  exponential  in  (1.56).  In  vh*; 
examples  to  follow,  only  a  single  term  is  indicated,  and  th'^ 
functions  v'^(4,T);b^)  and  are  eval'jated  for  an  arbi¬ 

trary  value  of  b.  The  results  thus  presented  correspond  to 
a  single  term  in  (1.56),  i..e.,  7(C)  =  f  e  .  A  more 

general  model  is  then  easily  treated  as  a  sum  of  terms,  as 
indicated  in  (1.58)>  (1.59)* 

Even  when  the  material  behavior  is  expressed  in  the 
simple  analytic  form  (I.56),  the  surface  displacement  in  a 
.general  three-dimensional  problem  cannot  usually  be  given 
in  simple  terras.  This  is  evident  even  in  the  most  elemen¬ 
tary  case,  that  of  a  concentrated  load  at  the  origin  (Example 
a) .  Therefore,  the  other  examples  will  be  two-dimensional 
(plane  strain)  problems,  for  which  results  can  be  given  in 
relatively  simple  expressions. 

A)  Three  Dimensions  Concentrated  Load 


<5(4. n)  =  26(4)6(n) 


vu.Ti;  = 


^  x.Du; 


The  integral  in  (l„60)  has  no  simple  evaluation  In 
i;he  general  case;  some  special  cases  are  given. 


For  r,  =  0: 


v(4>0)  =  ^  #  i  > 


as  ^  ^  +  00  f 


v(4>0)  ~  i  ( 1  +  £)  - 


,  i  <  0 

f  1  /I  2 

FT  ^ 


For  i  =  0: 


v(0,T))  »  +  f  I 


is  Weber’s  Bessel  function  of  the  second 
kind,  order  zero. 

Sq  is  Struve’s  function  of  order  zero. 


as  In]  —0,  v(0,ii)  ~  -  f  I  log  |b7)| 

as  in!--,  v(0.n)~(i^|)-i^-f^ 


?  fi - ^ 

3  (bn)' 


The  surface  displ? cement  is  infinite  at  the  origin  and 
all  along  the  negative  ^  axis;  at  all  other  points  it 
is  finite.  Some  typical  surface  profiles  are  sketched 
in  Fig.  1.10. 


Surface  \|/ 

Profiles  yj 

V 

i' 

!v 

-5  =  0 


Parallel  to  4-axis 


Parallel  to  “n-axls 
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Two  Dimensions 


Concentrated  Load 


Q(e)  =  ^U) 

v®(4)  -  v®(4q)s  log  (|2.|  ,  arbitrary 
v(C)  -  vdp)  =  (1  +  f)  log||^|  +  f  [e‘’^(-be)-e*’?OEj^(-bep)] 

»  -t 

E^(-x)  =  ~  J  the  Exponential  Integral. 

X 

as  III  —O  ,  v(4)  -  v(4^)  ~  log  ||^| 

as  |4|  -^00  ,  v(4)  -  v(4^)  ~  (1  +  |)  log  ||^| 


A  surface  profile  is  sketched  in  Pig.  1.11. 


Example  B  Kelvin  solid 

PIG.  1.11 


A  concentrated  load  is  an  Idealized  case  which  is  use¬ 
ful  for  some  considerations,  and  the  results  can  serve  as  a 
fundamental  solution  to  be  integrated  for  more  general  prob¬ 
lems.  At  distances  far  from  the  load,  the  results  for  any 
distributed  load  are  nearly  the  same  as  those  for  a  concen¬ 
trated  load.  Near  the  load,  however,  the  actual  distribution 
is  significant.  A  concentrated  load  is  a  singularity  which 
has  an  infinite  displacement  at  the  point  of  application. 

The  elastic  displacement  under  a  concentrated  load  is  also 
infinite.  Thus  all  points  of  the  viscoelastic  surface  that 
-have  passed  under  the  load  retain  infinite  displacement.  In 
the  two-dimensional  problem,  the  concentrated  load  is  actually 
a  line  load;  the  singularity  is  not  quite  as  pronounced.  The 
displacement  immediately  under  the  load  is  infinite,  due  en¬ 
tirely  to  the  initial  elastic  response.  However,  points  that 
have  passed  under  the  load  have  recovered  to  a  finite  dis¬ 
placement.  If  there  is  no  initial  elasticity,  the  displace¬ 
ment  in  the  two-dimensional  problem  is  finite  even  under  the 
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load.  For  example,  a  model  with  a  single  Kelvin  element 
gives: 

V(5)  -  v(l^)  =  log  1^!  +  Ej^(-t)5)-  e’^^oEj^(bej,)j 

for  1=0  ^  ^  _gb|o  E^(-bl^) 

A  sketch  of  this  displacement  profile  is  shown  in  Fig.  1.11. 

A  distributed  load  can  be  considered  an  integral  of  a 
concentrated  load.  The  singularity  is  no  longer  present, 
and  displacement  under  the  load  is  finite.  Two  simple  cases 
of  distributed  loads  are  given  in  examples  C)  and  D)  as  two- 
dimensional  problems.  The  resisting  force  and  coefficient 
of  friction  are  evaluated  in  ,.  3se  examples,  along  with  the 
surface  displacement. 

C)  Two  Dimensions  Uniform  Load 

<}(e)  =  Qo  -  0  <  i  <  -K 

N  =  Qo 

1  /  / 

v®{4)  H®(o)  =  [(e-^)iog|4f?v!  -  eiogiei  +  ?>  log  m 

c 

( 

\ 

i 

i-  [v(e)  -  V  0)]  =  (1  +  |)[{4->>)log|4-?i|  -  4log|4|  +  log  X] 

(1.62) 

+  ^  [log(b|4->.|)  -  log  b|^| 

+  e^^Ej;(-be)  -  log  bA  +  e"^^E^(b>,)  -  7  ] 
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The  low  point  (where  =  0)  always  occurs  for 

0  <  ^  <  ^  .  The  slope  ov/h^  is  infinite  at  ^  -  0 
and  4  =  A. 

X 

^  ^ 

O 

x  =  |  =  «oi^  (1 

where  h^(x)  =  [27  +  2  log  x-e^E^( -x) --e“^E^(x)  ] 

as  X  -*•  0  ,  h^(x)  ~  I  <2  -  7  -  log  X  +  I  +...)-*•  0 
as  X  —  «  ,  h^(x)  ~  i  (log  X  +  7  ~  -*-0 

X  X 

See  Fig.  1.12. 

7  =  .577216..  ,  Euler's  constant 


63) 
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D)  Two  Dlmenalons  Elliptical  Load 

ft;o  =  I  ^  vtt?'-?)  .  0  <  e  <  A 

2 

[v®(5)  -  v®(o)]  =  +  0  <  i  <  -h 

=  -  +  A4  +  i?  -  ||  vITF^ 

2 

“  ^  cosh  ^  |^-1|  ,  ^<0,  ^>1 

^  [v(0  -  v(0)]=  (1  +  (-4^  +  7v|) 

+  ^l-2i  +  (e^^-l)K^(|2l)]^  0<  ^  <>v 

(1.64) 

=  (!  +  §)(  +  U  ■  |l 


-  ^  cosh-^lfi  -  l| 

+  ^  [-2|  + 

2eA-i 

- 1  bx2  e'=(e-V2) 


-  l)Ki(|^) 


du] 


lower  limit  of  Integral  is 


r  +  1  for  I  >  A 
^  -  1  for  I  <  0 
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Slopes  are  finite  at 

e  -  0,  X 


Low  point  occurs  for 

f 

0  <  4  <  I 


where  h^(x)  =  ^  [1  -  2K3^(x)  I^(x)] 


as  X  -  0  ,  h  (x)  ~  §  ( .6159  -  log  x)  0 


Some  results  of  examples  C)  and  D)  indicate  general 
conclusions  that  apply  to  other  distributed  loads.  The  dis¬ 
placement  relative  to  4  =  0,  v(|)  -  v(0)  ,  Is  everywhere 
finite  In  the  vicinity  of  the  load,  and  increases  like 
-log|4|  as  |4|  -♦  00  .  Example  C)  has  a  discontinuity  In  the 
pressure,  a  sudden  Jump  from  zero  to  ,  at  each  end  of 
the  loaded  area.  This  results  in  an  Infinite  slope  3v/34 
at  each  end  where  the  Jump  occurs.  This  is  a  general  result; 
if  there  is  a  finite  Jump  in  the  load  at  4  =  >  -he  slope 

will  contain  a  term  log 1 4-^0 1  »  will  thus  be  infinite 
at  4q  •  example  D),  the  pressure  is  continuous  and  goes 
to  zero  at  the  ends;  the  slope  is  finite  everywhere.  In  each 
example,  the  pressure  distribution  is  symmetric  about  the 
center  of  the  loaded  area  (4  =  ?^/2)  .  However,  the  displace¬ 
ment  is  not  symmetric  because  of  delayed  recovery  in  the  mov¬ 
ing  viscoelastic  base.  The  low  point  in  the  surface  profile 
will  always  be  displaced  downstream  from  the  center  of  the 
load.  The  velocity  is  conveniently  representea  in  dimension¬ 
less  form  by  VT/i  =  1/A  .  The  variation  of  resistance  with 
velocity  for  a  given  load  is  then  shown  by  a  plot  of  x  vs 
1/A  .  For  the  two  examples,  this  variation  is  Indicated  by 
the  plots  of  h^(x)  and  hgCx)  in  Fig.  1.12,  Each  example 
shows  that  x  increases  from  zero,  reaches  a  single  maxi¬ 
mum,  and  then  decreases  to  zero  again,  as  the  velocity  goes 
from  0  to  00  . 
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CHAPTER  II 


MOVING  LOADS  ON  A  PURELY  VISCOUS  MATERIAL 

The  problem  of  a  load  moving  on  the  surface  of  a  purely 
viscous  material  Is  considered  in  this  chapter.  By  “viscous 
material”  is  meant  an  ideal  incompressible  viscous  fluid,  in 
which  stress  is  linearly  related  to  strain  rate.  This  is  a 
special  case  of  the  general  viscoelastic  fluid  considered  in 
Chapter  I.  Such  materials  as  dough  or  road  asphalt  may  some¬ 
times  be  treated  as  purely  viscous.  Even  metals  may  be  con¬ 
sidered  viscous  fluids  when  subject  to  high  speed  impact 
(Abrahamson  [1]).  The  problem  is  treated  as  one  of  slow 
steady  motion  of  a  very  viscous  fluid,  so  that  inertia  forces 
are  negligible  and  linear  theory  can  be  used.  As  in  Chapter 
I,  the  load  q(x,y)  will  be  normal  pressure  on  the  initially 
flat  surface  z  =  0  of  a  base  of  infinite  extent  (see  Fig. 
1.1). 

The  principal  interest  will  again  be  the  surface  defor¬ 
mation  as  the  load  moves  with  constant  velocity  V  in  the  x 
direction.  The  base  is  undisturbed  for  t  <  0  ,  and  the 
load  is  applied  at  t  =  0.  The  resulting  deformation  will 
vary  with  time  due  to  both  the  motion  of  he  load  and  the 
viscous  nature  of  the  base.  The  relative  simplicity  of  the 
purely  viscous  material  allows  solution  of  some  three- 
dimensional  problems.  These  give  some  idea  of  the  much 
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greater  difficulties  encountered  in  the  same  problems  with 
a  general  viscoelastic  material.  In  Section  1.4  it  was 
shown  that  the  question  of  a  steady  state  for  a  general  vis¬ 
coelastic  fluid  can  be  an'jwered  by  considering  the  same  prob¬ 
lem  for  a  purely  viscous  material.  In  particular,  since 
there  is  no  steady  state  for  a  semi-infinite  base,  the  pos¬ 
sibility  of  a  steady,  when  the  base  is  a  viscous  layer  of 
finite  thickness  is  considered. 

2.1  General  Solution  for  Viscous  Base 

For  an  incompressj^le  viscous  fluid  with  viscosity 

A 

the  relations  between  st:?ess  and  strain  are 

®ij  “  ®lj  ^  ®  ° 

For  an  Incompressible  elastic  solid,  with  v  =  ^  , 

^ij  =  ^  e  =  0  (2.2) 

In  a  quasi-static  problem  with  prescribed  tractions  on  the 

boundaries,  equilibrium  and  boundary  conditions  are  the  same 

for  elastic  and  viscous  materials  at  any  instant  of  time. 

Also,  for  an  elastic  material  e. .  =  l/2(u.  u.  .),  while 

i  j  J  j  >  1 

for  a  viscous  material  1^*  Comparison 

of  these  relations  makes  evident  the  "viscous  analogy"  for 
slow  motion  of  a  viscous  fluid:  the  viscous  solution  is  the 
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same  as  the  incompressible  elastic  solution  with  |i  replaced 
by  Tj*  ,  V  =  1/2  ,  and  strains  replaced  by  strain  rates. 

In  particular,  elastic  displacements  become  viscous  veloci¬ 
ties.  This  viscous  analogy  corresponds  to  the  use  of  the 
associated  elastic  problem  for  general  viscoelastic  materials. 

In  moving  load  problems,  the  surface  displacement  for 
an  Incompressible  elastic  material  is  u®(x*-Vt,y*)  in 
coordinates  (x*,y*)  fixed  with  the  base,  where  u®(x*,y*)  = 
(l/2|i)  P(x*,y*)  is  the  displacement  for  the  same  load  at 
rest.  Thus,  from  the  viscous  analogy 

Ig-  U2(x*,y*,t)  =  P(x*-Vt,y*) 

Since  the  surface  is  undeformed  at  t  --  0,  u_(x*,y*,0)  *  o. 

The  viscous  displacement  is  then 

t 

U2(x*,y*,t)  =  ^  J  P(x*-Vt’ ,y*)dt  ‘ 

o 

Changing  to  coordinates  (l.l)  moving  with  the  load 

t 

U2(x,y,t)  =  ^  J  p(x+  Vt  -  Vt'  )  at’ 

o 

x+Vt 

UgCx.y.t)  =  J  P(x',y)  dx' 


(2.3) 


In  theiic  coor-ulnates,  it  will  be  more  convenient  to  tninic  of 
the  load  at  rest  (with  the  origin  included  inside  the  loaded 
area),  and  the  case  moving  under  it  with  velocity  V  in  the 
negative  x-dlrection.  The  base  may  be  either  semi-infinite 
or  a  layer  of  finite  thickness. 

The  general  solution  for  surface  displacement  is  given 
by  equation  (2.3).  This  can  be  expressed  in  terms  of  a 
function  G(x,y),  defined  by 

X 

o 

Since  P(x‘,y)  is  (except  for  a  constant  factor)  the  elas¬ 
tic  displacement  for  the  given  load  distribution,  G(x,y) 
represents  the  area  between  the  initial  and  deformed  sur¬ 
faces  of  an  elastic  base,  from  x‘=0  to  x'«x,ina 
slice  parallel  to  the  x-axis  at  the  given  value  of  y  (see 
Pig.  2.1).  Then 

u  (x,y,t)  =  G(x+Vt,y)  -  G(x,y)  (?.5) 

The  term  G(x+Vt,y)  represents  a  "wave"  moving  with  the 
base.  The  other  term  G(x,y)  is  a  displacement  of  the  same 
shape,  fixed  relative  to  the  load.  The  difference  of  these 
terms  gives  the  net  displacement,  which  varies  with  time  as 
the  "wave"  is  carried  downstream  (see  Fig.  2.2).  If  the 
area  under  the  complete  elastic  curve  is  finite,  so  that 
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G 
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G„  »  0  / 

G  >  0 

00 

\y 

u 

“z 

O  O  vCt  vw 

profiles 

z 
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2.4a 

PIG.  2.4b 

nf _  Vie»0  o  ««  /v  «♦  _ _ — '14«m4^ 

'4\w^jr/  v9.^wi^^  Oppl'UOUiiCtJ  ®  xxmxt/”^ 

ing  value  as  t  -»  «  ,  and  there  is  a  steady  state.  If 
G(oo,y)  is  Infinite,  there  is  no  steady  state. 

When  the  elastic  displacement  is  always  positive  (this 
is  the  case  for  a  semi-infinite  hase),  G(x,y)  is  a  monotonic 
Increasing  function  as  shown  in  Pig.  2.2.  When  the  base  is 
a  layer  of  finite  thickness,  the  elastic  displacement  may  be 
negative  (rise  above  the  initial  surface)  at  some  distance 
from  the  load.  Then,  G(x, y)  will  have  the  form  shown  in 
Pig.  2.3.  Let  G  (y)  =  G(3c,y)j  there  is  a  steady  state 

00  X  “•*00 

if  this  limit  exists.  If  G  >  0  ,  the  steady  state  dis- 

00 

placement  is  u  (x,y)  =  G  (y)  -  G(x,y).  The  result  is  a 
depression  of  finite  depth  extending  infinitely  far  down¬ 
stream  from  the  load,  and  zero  displacement  far  upstream 
from  the  load  (Plg.  2.4a).  In  the  case  shown  in  Pig.  2.3, 

it  is  possible  to  have  G  =0.  Then,  the  steady  state  is 

00 

u_.  -G(x,y).  The  displacement  is  zero  far  from  the  load 

in  any  direction,  and  a  definite  localized  hump  is  left  near 
the  load  (Pig.  2.4b). 

When  the  load  is  symmetric  with  respect  to  the  y-axis, 
so  that  q(x,y)  =  q(-x,y),  the  elastic  displacement  will  be 
symmetric  (thus  even  in  x) :  P(x,y)  =  P(-x,y).  Tien,  the 
function  G  will  be  antisymmetric  (odd  in  x) :  G(x,y)  = 
-G(-x,y).  In  this  case. 
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U2(-x-Vt,y)  =  G(-x,y)  -G(-x-Vt,y) 

=  -G(x,y)  +  G(x+Vt,y) 

=  u^Cx^y) 

Thus,  the  viscous  displacement  -will  be  syiTimetric  with  respect 
to  a  line  parallel  to  the  y-axis,  halfway  between  the  fixed 
origin  (x=0)  and  the  point  on  the  base  originally  at  the 
origin  (x  =  -Vt).  Thi:  line  of  symmetry,  to  be  called -the 
trough  line,  is  then  x  =  -  Vt/2  .  It  moves  downstream  at 
half  the  velocity  of  the  base  material.  The  displacement 
along  the  trough  line  is 

I  =  2G(|  Vt,y)  (2.6) 

If  G(x,y)  is  a  monotonic  function,  this  is  the  maximum 
displacement  for  given  y  ?nd  t  (see  Fig.  2.5), 
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For  a  pressure- q(x,y)  on  a  semi -infinite  base,  the  dis¬ 
placement  function  is 


P(x,y) 


1_ 

27r 


q(x!,y»)  dx'dy* 

(x-x‘)^+(y-y') 


2 


The  displacement  is  always  positive  [for  q(x,y)  >  0],  with 
a  maximum  near  the  load,  and  approaches  zero  far  from  the 
load  (see  Fig.  2.1  for  typical  profile).  For  large 
r  =  v^- 2 ,  2  ,  P(x,y)  is  nearly  the  same  as  for  a  concen- 

X  +y 

/  ,  N 

trated  load,  i.e.,  P(r)  ~  2??  r -►  »  ,  where 

N  =  //  q{x,y)dxdy.  Thus,  the  function  G(x,y)  is  nearly 
A 

the  same  far  from  any  distributed  load  as  for  a  concentrated 
load  N  at  the  origin.  It  is  useful,  therefore,  to  consider 
the  case  of  a  concentrated  load  in  some  detail. 

For  a  load  N  at  the  origin. 


P(x  v)  -  ^  _ ^ 

nx,yj  -  2Tr  r  ”  27r  /  x- ■  p 


Then,  from  (2.4) 


G(x,y)  = 


N 


4TrTi*V 


I 


dx 


N 


o  "^  x  ’  ^+y^ 


4Trr)*V  * 


[x  -tyx^-tv^ 


y 
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Prom  (2.5), 


47rn*V 
N  ^ 


U2(x,y,t)  =  log  ^ 


f  x+Vt  +  v/f  x-hVt  r  +  7 


*  +  '/x2+y2 


(2.7) 


For  a  given  y,  G.(oo,y)  is  infinite,  confirming  the  conclusion 
in  Chapter  I  that  there  is  no  steady  state.  At  any  point  on 
the  surface,  as  t  -♦  «  , 


-  U2(x,y,t)  ~  log  Vt  -  log(x  +v^x^y^  +  log  2  +  ^  +  .. 

(2.8) 


Thus,  the  displacement  increases  without  limit  as  t  -►  ». 

It  can  be  shown  that  the  vertical  velocity  ^u  /^t  of  the 
surface  approaches  zero,  and  the  surface  slopes  (Su^/^x, 
du  /^y)  approach  finite  values  t  -*•  «  . 

For  a  given  finite  time,  the  displacement  (2.7)  is 
positive  everywhere;  it  approaches  zero  far  from  the  load, 
as  |x|  or  |y|  -*>  «,  Thus,  appreciable  Sv.rface  deformation 
occurs  only  in  the  vicinity  of  those  points  on  the  surface 
that  have  passed  under  or  near  the  load.  The  displacement 
is  symmetric  about  a  trough  line  at  x  =  -  Vt/2  ,  as  dis¬ 
cussed  in  Ccction  2.2.  The  maximum  displacement  at  any  given 
y  and  t  is,  from  (2.6) 


u. 


max 


N 

27rTi*V 


N 

2i7-n*v 


tanh~^[-- - . 

^+{2y/Vt)^ 


(2.9) 
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Because  the  load  is  concentrated  at  a  point,  the  elastic  dis¬ 
placement  is  infinite  at  that  point.  The  viscous  displace¬ 
ment  is  thus  infinite  at  all  points  that  have  passed  directly 
under  the  load,  i.e.,  for  y  =  0  ,  -Vt  <  x  <  0.  The  contour 
lines  of  constant  displacement  at  a  given  time  are  ellipses, 
with  common  foci  at  the  initial  and  present  load  positions: 
y  =  0;  X  =  0,  -Vt.  The  contour  ellipse  for  (^irr) *V/N)u  =  P 
has  a  constant  eccentricity  e  =  (e^-iy(e^+l) ,  and  an  expanding 
major  axis  Vt/2G .  Figure  2.6  shows  typical  profiles  and  a 
plan  view  of  the  surface  for  a  concentraced  load. 

Except  near  points  that  have  passed  under  the  load,  the 
results  for  a  distributed  load  (symmetric  about  the  y-axis) 
are  qualitatively  similar  to  those  for  the  concentrated  load. 

The  general  shape  of  the  surface  is  like  a  shallow  bowl,  with 
displacements  decreasing  monotonies lly  in  any  direction  from 
the  low  point  at  y  =  0,  x  =  -  Vt/2.  The  surface  approaches 
its  initial  flat  position  far  away  from  the  loaded  regions. 

As  time  goes  on,  the  displacement  increases  without  limit, 
but  at  a  continue] ly  slower  rate.  For  the  semi-infinite  base, 
there  is  no  localized  "hump"  formed  near  the  load,  but  rather 
a  downstream  depression  that  is  continually  deepening  and 
moving  away  from  the  load.  The  contours  of  constant  dis¬ 
placement  will  be  simple  closed  curves,  symmetric  about  the 
trough  line.  To  demonstrate  these  results  for  a  specific  dis¬ 
tributed  load,  an  example  will  be  considered  In  the  next  section. 
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PIG.  2.6 

Surface  displacement,  fixed  t  (eqn.  2.7) 


2.3  Uniform  Pressure  on  a  Circular  Area 

As  a  three-dimensional  example  of  a  distributed  load, 
a  uniform  pressure  acts  on  a  circle  of  radius  a. 

Then, 

q(x,y)  =  q(r)  =  q^  ,  0  <  r  <  a 


where  r  -J  2.  c  .  The  viscous  base  is  semi-infinite. 
^  X  +y 

The  following  dimensionless  quantities  are  used; 


(2.10a) 


«(p)  =  ^ 

g(e.n)  =  ^  G(x,y) 


v(?,n.T)  =  u  (x,y,t) 


(2.10b) 

(2.10c) 

(2.10d) 


The  elastic  displacement  can  be  expressed  in  several 
ways.  ^  The  most  convenient  form  for  use  here  is  derived  using 
the  Hankel  Transform  (Sneddon  [28] y  p.  ^69) y  and  gives 

00 

w(p)  =  J  J;^(u)  J^(pu)  2“  (2.11) 

o 

This  integral  is  a  special  case  of  the  Weber-Schaftheltlih 
discontinuous  integral  (Watson  [35],  p.  398),  and  has  the 
series  form 
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(2.12) 


00 


J  J^(u)  Jq(pu)  ^  =  2^1  ^  f  ^  f  P  1  ^ 


o 


:^1  (|  ^  I  i  2  i  p  >  1 


2‘1  ^  2 
2 

T  * 


P  =  1 


where 


(a,b;  c;  x)  = 


1  +  I  ri^  r,Lc), 

k=l  ^Tar  r(b)  r(c+k)  IT. 


is  the  hypergeometrlc  function.  Thus, 


00 


.2k 


w(p)  =  1  -  Sc,  p^^^  ,  p  <  1 
k=l  ^ 


=  ^  ^  ’’k  ,  p  >  1 


2p 


k=l 


TT 


P  =  1 


(2.13) 


where 


°k  = 


1  [r(k4)]g 

^  (k')^(2k-l) 


1m2 


1  [r(k-f|)] 

^  ”  2Tr  /^,i2 


2k- 1 


(k')‘^(k+l)  2XktlT  k 

(2.14) 


(Numerical 'values  of  C  and 

k 

are  given  in  Table  2.1  at 
the  end  of  the  chapter.)  A  pro¬ 
file  of  the  elastic  displace¬ 
ment  is  sketched  in  Fig.  2.7. 
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Since  p(x,y)  =  2ix,  u®(r)  =  q^a  w(p)  ,  (2.4)  and  (2.10c)  give 

Z  C/ 

g(«,n)  =  J  w(  )  dC  (2.15) 

o 

and  (2.5)  and  (2.10d;  give 


v(4>il^'r)  =  6(^'H^,ti)  -  g(CiTi)  (2.16) 


The  function  g(4i‘n)  is  independent  of  any  physical  param¬ 
eters,  so  the  numerical  values  need  be  determined  only  once 
for  all  problems  of  uniform  circular  load.  Then,  the  viscous 
displacement  at  any  point  on  the  surface  for  any  time  is 


U2(x,y,t) 


^0^^  f„/X+Vt 
2q*V  I®'  a 


|)  - 


Because  there  are  two  different  expressions  (2.13)  for 
w(p),  one  for  p  <  1  and  one  for  p  >  1  ^  there  is  more 
than  one  form  of  g(^,Ti)  evaluated  from  (2.15).  Since 

=  -  6(-^^Ti)  and  g(4>n)  =  g^^-T]),  it  is  sufficient 
to  evaluate  g(^,ri)  for  (^j'H)  in  the  first  quadrant  only. 
There  are  three  forms  for  ,  one  for  each  of  the 


three  regions  shown  in  Fig.  2.8. 
Numerical  values  of  g(^,T))  can 
be  determined  from  the  expres¬ 
sions  which  follow.  Some 
values  are  given  in  Pig.  2.9, 


and  contour  lines  are  sketched 


PIG.  2.8 


in. 
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(  P  ^  ^  <  <s/  l-ll^ 


=  gi(?,n)=  /[I  -  z  cJc^+Ti2)^]dc 
i  k=i  ^ 

=  ?  -  f  ‘^k/  (CW)*'  d? 


Let 


ii(«.n)  .  /(cW)"  dc  =  (i)"' 


2k+l 


Then 


“  5raT  ? 


Si(?.n)  =  «  -  2  c^ijce.n) 

K=1 


= «  -  ?  sSr 


(2.17a) 

(2.17b) 

(2.18a) 

(2.18b) 


(2)  P  :^i  ,  n  >  1 


■id.n)  =  gp(4.n)  = 


=  /  [„7^+  2  D  (c2+Ti2)‘'‘‘4MdC 

d  '•2v^‘^+r|2  k=l  J 


? 


=  5  log  ^  +  J^.  y  {cW)"^.‘'4> 


dC 


Let 


4  (5.1)  =  /(cW)‘^‘''*i^  dC  ,  (n  >  1)  (2.19a) 
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With  the  c-uhstitution  C  =  cot  0  ,  this  becomes 

■w/2 

(t.n)  =  J  <P  a<t>  (2.19b) 


where 


9  -  tan”^  ^  (see  Pig.-  2.8). 


For  numerical  evaluation,  this  can  be  expressed 
In  terms  of  the  incomplete  Beta  function  B  (a,b). 
With  X  =  sin^e  (0  <  ®  <  5)^ 

9 

sin^^"^  0  d0  =  i  B^  (k,|) 


I 


2  X 


Then,  (2.19b)  becomes 

J-l  (e.n)  =  I  CB3^(k,|)  -  Bx(k,|)]  (2.19c) 

Values  of  tabulated  in  Pearson  [28]  for 

X  =  0( .01)1.00,  k  =  1(1)50. 

_  -r2/^  1  ^-2k„  lx  1  ^-2k  .  2k+l  ^2^-{2k+2) 

Ase-co,  ij^(^,T])  ~  -  Ti  -Bi(k,^)-^e  +7r(¥FiT^^  ..: 

As  e  -*•  0,  i^(e,Ti)  ~  T]"(2k+3)|3  _ 

Using  (2.19a,c), 

gi^(e.'l)  =  I  log  ^  \  4  (2.20a) 
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=  -  5  log  tan  I  +5  2 

(2.20b) 

2 

where  x  =  sin  0. 

For  a  given  0,  as  p-^oo,g  ~-i  log  tan  S  +  t 

^  ^  ^  ID 


(3)  p>l,0<ri<l,e>  v'^ 

^  I 

=  J  w(  +  y'  w( 

°  y — 2 

Vl-T] 

=  gi(v^  -  n)  +/  t  s  D^(cV)-'''4)]dc 

vi-n^ 


Let  ilu.n)  =  y  dC  .  (n  >  1. 

5 

V  1-T] 


=  I  [B^2(>'4)-V<*'’|)^ 


(2.21a) 


where  s  =  ii/p 


ifeU-o)  =  iic  (1  - 


(2.21b) 


As  C  -♦  00  ,  I^(C,Tl)  ~  ^  ^  +  .  . . 


I 


,ri  /  .  \  /  n — ~/y  _  \  ,  1  ,  .  /  ^  +  O  V 

•men,  goU??V  =  gi  W-L“'l^»^U-<'  P  \-^‘-—z:=j 

^  ~  i+vT:;^ 


4  2  B.  (C,T)) 

k=l  ^  ^ 


(2,22a) 


g3(4.0) 


»  Cv  D, 

1-2  (  ^ 

,  -^2k+l  2k-' 

k=i 


°®  -2k  1 

'  2  ^  I  +  i  log  ^ 

k=l 


(2.22b) 


As  ^  -*•  00  for  a  given  t),  both  gg  and  g,;^  behave  like 
1/2  log  ^  and  increase  without  limit.  This  again  confirms 
that  there  is  no  steady  state. 

There  are  several  regions  of  che  surface,  each  having 
a  different  expression  for  the  displacement  v(^,T),'r)  = 
g( ^4-1,11)  -  g(4,t))  ,  because  of  the  different  forms  of  g. 
These  regions  are  shown  in  Fig.  2.10a  for  0  <  t  <  2  , 
and  Fig.  2.10b  for  t  >  2.  Displacements  are  symmetric 
about  the  ^  axis  and  the  trough  line  ^=-t/2,  so  only  the 
quadrant  t]  >  0,  i  >  -  t/2  is  considered. 


) 

II 

n 

/" 

LA 

'v  ^ 

j 

* 

PIG.  2.10a 


PIG.  2.10b 
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Points  always  under  the  load  (beginning  at  t  --  0) . 


-  g^(4,Ti) 


II.  Points  never  under  the  load. 


=  g.;)(^-K,Ti)  -  g^(^,r,)  for  T]  >  1 


=  g3(^+'r>il)  -  for  r|  <  1 

(Both  cases  give  the  same  analytic  expression. ) 

III.  Points  originally  outside  the  load  which  have  moved 
under  it. 


v(^,ri,T)  =  g3(?+T,Ti)  - 

IV.  Point  originally  outside  the  load  which  then  moved 
under  it  and  are  now  outside  again. 

v(^,t),t)  =  +  g^viehn) 

Displacements  and  slopes  are  continuous  across  the 
boundaries  of  these  regions.  The  maximum  displacement 
occurs  at  i  t]  =  0: 

T  ^Ir  r  2k  1 

%ax^  2g(^,0)  =  T  1  -  2.  (2)  I  ^  ^ 


2^1  -  A  21^+1  ~  2k'  ,  ■  2k  ^2^  2  " 


T  >  2 


TnS  dcfoi'mcd  SUrfSCc  ics  Smooth^  Wioil  riuitc 

displacement  everywhere  for  finite  t.  As  t  -►  «, 

^  log  at  all  points  on  the  surface.  The 

general  description  in  the  last  paragraph  of  Section  2.2 

applies  in  particular  to  this  problem,  and  conclusions  given 

there  can  be  verified  directly  from  this  section.  At  any 

finite  time,  the  displacement  far  from  the  load  is  nearly 

2 

the  same  as  for  a  concentrated  load  N  =  q  ira  . 


2.4  Base  of  Finite  Thickness 

■  *  II  .1  . . .  ...........  R  I  . 

When  the  base  is  semi-infinite,  the  viscous  material 
flows  without  limit  due  to  the  vertical  velocity  resulting 
from  the  load  on  the  surface.  There  is  no  lower  boundary 
to  the  base,  so  the  vertical  displacements  continue  to  in¬ 
crease  with  time  and  no  steady  state  is  reached.  If,  however, 
the  base  had  a  lower  boundary,  so  that  flow  is  constrained  at 
some  depth  below  the  surface,  vertical  displacements  would 
not  increase  indefinitely.  It  seems  quite  likely  on  physi¬ 
cal  grounds  that  this  situation  would  lead  to  a  steady  state 
when  the  base  is  moving  horizontally  under  the  load.  The 
base  of  finite  thickness  does  of  course  represent  a  real 
situation  more  closely  than  does  the  semi-infinite  base. 

The  thickness  of  the  base  is  thus  a  critical  factor  for 
estabiis'nment  of  a  steady  state  in  problems  of  this  type, 
when  the  base  is  a  fluid  material 


To  investigate  this  situation  more  closel.y,  the  base 
will  be  considered  as  a  layer  of  constant  finite  thickness 
h  ,  but  extending  infinitely  far  in  the  x  and  y  directions. 

As  before,  the  load  q(x,y) 
is  applied  to  the  otherwise 
free  surface  z  =  0  ,  and  the 
viscous  layer  moves  with  velo¬ 
city  V  in  the  negative  x  direc¬ 
tion  (Fig.  2.11).  The  bottom 
of  the  layer,  z  =  h  ,  rests 
on  a  rigid  support,  so  that  vertical  velocity  and  displace¬ 
ment  are  prevented.  One  condition  is  thus 

U2(x,y,b)  «  0  on  z  =  h  (2*23) 

The  way  in  which  the  layer  is  supported  at  the  lower  boundary 
is  quite  important.  The  two  extremes  are  considered  herej 
other  possibilities  will  fall  in  between.  One  extreme  is 
rigid  attachment  so  that  tangential  displacements  are  pre¬ 
vented.  This  gives  the  conditions  (for  a  moving  base) 

u^(x,y,t)  =  -V,  u„(x,y,h)  =0  on  z  =  h  (2.24a) 

-  y 

The  other  extreiae  is  no  attachjnent  (i.e.,  frictionless 
support),  so  there  will  be  no  tangential  stresses.  The.n, 

■r.^2  =  (-x,y,h)  =  ^0  on  z  =  h  (2.24b) 


q(x.y) 

_ 


//  /  y  ~~77V777'7-/^7-ry'7~^/  .  jrnif 


PIG.  2.11 
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The  problem  of  an  axlsymmetric  load  will  be  considered 
in  more  detail.  The  results  for  other  distributed  loads 
would  bo  qualitatively  similar.  The  load  is  q(r)  on  a 
circle  of  radius  a.  The  dimensionless  quantities  (2.10a-d) 
are  used.  Also,  let  q(r)  =  q^  Q(p)  ,  where  Q(p)  is  di¬ 
mensionless  and  q^  is  some  measure  of  the  load  magnitude. 
Following  the  analysis  of  Sneddon  [28],  the  incompressible 
elastic  solution  is  determined  using  Hankel  Transforms.  The 
resulting  surface  displacement  is 


For  a  uniform  load  q^,  Q(p)  =  1  and  Q(u) 


(2.25) 


(2.26) 


The  function  F(>xu)  depends  on  the  method  of  support.  For 
rigid  attachment  (2.24a) 


F(x) 


cosh  X  sinh  x  -  x 
2 

X  +  cosh‘'x 


(2.27) 


~  ^  x^(  l-^x^+...)  asx-^0 
~  1  -  2(  l+2x+2x^)e“^^+. . .  as  X -*•  00 


85 


~  1  -  2(  1  +  2x)e~^^"+. . .  as  x  -►  » 

In  both  cases,  F(x)  increases  monotonically  from  zero  and 
approaches  1  asjmiptotically  as  x  goes  from  0  to  <». 

For  a  semi-infinite  base,  A  =  »  and  F(Au)  =1  in  each 
case.  Then,  for  a  uniform  load,  (2,25)  reduces  to  (2.11). 

From  the  elastic  displacement  [2 ,2j) ,  equation  (2.15) 
gives  4 

gCen)  ^-  [  \  f  «(u)  Jo(  (2-29) 

J  I  J 
o  o 

The  possible  steady  state  depends  on  SgpC'n)  =  b(  Cj  t)  * 
Letting  ^-*■00  in  (2.29),  and  using 

Ic —  1 

00 

/j^(u?)de  =  i 

o 

the  (  integration  in  (2.29)  can  be  carried  out.  This  gives 

00 

=  J  5('>)  F(’'u)[J^(ut,)  +2  2^ 
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The  integrand  of  (2.30)  is  finite,  continuous,  and  behaves 
like  (m  >  1)  as  u  “*■  00  .  Although  (2.30)  cannot 

be  evaluated  analytically,  it  can  be  verified  that  goo('n) 
has  a  finite  value.  Thus,  the  base  of  finite  thickness  does 
reach  a  steady  state. 

The  nature  of  the  steady  state  depends  on  the  value  of 
discussed  in  Section  2.1.  This  would  require  a 
numerical  evaluation  of  the  Integral  (2.30),  or  of  the  sur¬ 
face  displacement  and  the  area  it  encloses.  For  frictionless 
support,  F(Xu)  of  equation  (2.28)  can  be  approximated  by 
exponential  terms  (Sneddon  [28]),  and  the  resulting  integrals 
can  be  evaluated.  This  gives  goo(Tl)  >  0  for  any  q  ^  thus 
Indicating  a  steady  state  with  a  finite  depression  downstream 
(Fig.  2.4a). 

The  two-dimensional  problem  can  be  treated  in  the  same 
way.  The  elastic  displacement  can  be  determined  using  Fourier 
transforms.  The  result  is  similar  to  the  axisymmetric  prob- 
lejTij  for  a  layer  of  finite  thickness,  the  expression  for  dis¬ 
placement  (corresponding  to  equation  2.25)  is  the  same  as 
for  a  semi-infinite  base  (corresponding  to  equation  2.11)  with 
the  additional  function  F(Au)  .  For  the  same  support  con¬ 
ditions  (2.24a,b),  F(au)  is  the  same  as  for  the  axisymmetric 
problem  (2.27,  2.28).  Thus  it  is  reasonable  to  expect  that 
the  effect  of  the  kind  nf  support  at  the  lower  boundary  (z  =  h) 
on  the  nature  of  the  steady  state  will  be  similar  for  two- 
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dimensional  and  three-dimensional  problems.  The  two-dimen¬ 
sional  problem  of  a  uniform  load  on  a  base  of  finite  thick¬ 
ness,  rigidly  attached,  has  been  solved  by  Abrahamson  and 
Goodler  [2] .  This  solution  gives  g  =  0  ,  and  the 

.  C4 

resulting  steady  state  is  a  localized  hump  near  the  load 
(similar  to  Figi  2.4b). 

These  examples  suggest  general  conclusions  about  the 
steady  state  displacement  for  any  load  distribution  on  a 
base  of  finite  thickness.  When  the  base  is  rigidly  attached 
at  the  lower  boundary,  S^(u)  =  The  deformation  would 
be  IdcaTlzed,  forming  a  hump  at  the  load  (Fig*  2.4a).  When 
the  lower  boundary  is  not  attached  (i.e.>  frictionless 
support),  g  (fj)  >  0.  There  would  be  general  deformation 
downstream  from  the  load,  forming  a  depression  of  finite 
-depth  (Fig>  2.4b).  The  overall  form  of  the  steady  state 
diSpia; cement  is  thus  determined  by  the  way  in  which  the 
base  is  supported  at  the  lower  boundary.  Regardless  of  the 
kind  of  support,  a  purely  viscous  (and  thus  any  viscoelas¬ 
tic  fluid)  material  of  finite  thickness  will  always  reach 
a  steady  state. 


COEFFICIENTS  FOR  UNIFORM  CIRCULAR  LOAD  DISPIACEMENT 

(Section  2.3) 
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CHAPTER  III 


MO^NG  CONTACT  PROBLEMS  | 

I 

When  a  moving  body  of  given  shape  is  pressed  into  another 
(stationary)  body,  the  problem  is  one  of  "moving  contact." 

There  is  some  load  distribution  within  an  area  of  contact 
where  the  deformed  surfaces  of  the  two  bodies  must  match. 

If  one  or  both  bodies  is  viscoelastic,  dissipation  of  energy 

I 

■will  produce  res'istance  to  the  moving  body,  ^nd  in  general  I 

1 

i 

the  contact  conditions  will  change  with  time.  | 

\ 

1 

In  this  and  following  chapters,  a  particular  viscoelas-  | 

tic  moving  contact  problem  will  be  considered.  A  rigid  body 
moves  with  constant  velocity  V  in  the  x-direction  over  the 
initially  flat  surface  (z=0)  Of  a  semi-infinite  viscoelastic  I 

I 

solid i  It  is  held  in  contact  by  a  constant  normal  fores  N*,  | 

j 

which  produces  a  contact  area  A  of  unknown  shape,  extent,  and  | 

placement.  The  surfaces  are  "frictionless"  (i.e.,  no  tangen-  I 

tial  stresses)  and  non-adhesive,  so  the  load  distribution  is 
normal  pressure  q(x,y)  >  0  within  A.  The  moving  contact  is 
thus  a  moving  load  problem  as  discussed  in  Chapter  I,  although 
in  the  contact  problem  the  load  q  is  initially  unknown.  It 
was  shown  in  Section  1.4  that  a  steady  state  is  always  possi¬ 
ble;  only  these  will  be  considered.  Then,  the  contact  area 
and  pressure  do  not  change  with  time.  Once  they  are  deter¬ 
mined,  the  "friction"  force  F*  (which  resists  the  motion) 

can  be  found  as  in  Section  1.5- 

90 


It  is  again  convenient  to  consider  the  rigid  body  sta^ 
tionary  (in  Xjy,z  coordinates),  with  the  base  moving  in  th^; 
negative  x^direction  at  velocity  V  (see  Fig.  3.1) •  For  some 
fixed  orientation  of  the  body,  its  lower  surface  has  a  known 
shape  given  by  z  =  constant  +  w*(x-b,y-c)  ,  where  (bjc)  are 
the  coordinates  of  a  fixed  (but  arbitraiy)  reference  point-  E 
in  the  body.  The  body  is  held  in  contact  with  the  moving 
base  by  a  load  N*  normal  to  the  initially  flat  Surface  of 
the  base  (i.e.,  vertical)  and  any  horizontal  load  F*  (the 
resisting  or  ’'frictibn”  force)  that  may  be  hscessa:^^  For 
a  given  N*,  the  size,  shape,  and  plaeemMt  of  the  contact 
area  A  is  initially  unknown>  except  that  it  must  of  course," 
coincide  with  some  part  of  the  surface  of  the.  rigid- body. 

In  the  actual  ana lysis >  it  will  be  convenient  to  allow  the'- 
rlgid  body  to  rotate  through  a  small  angle  about  the  hori^ 
zontal  axis,  from  the  orientation  implied  by  the  function 


■ 

iH| 

in 

w*(x,y).  The  aiiouaj;  and  axis  of  this  slight  tilt  Is  also 
Initially  unknown. 

To  formulate  a  definite  problem,  the  conditions  of  con¬ 
tact  must  be  further  specified.  There  is  no  adhesion  between 
the  contact  surfaces  so  q(x,y)  >  0  inside  A.  If  the  sur¬ 
face  were  sufficiently  uneven,  breaks  in  contact  would  occur 
within  the  overall  ax^ea  whenever  complete  contact  required 
negative  pressure  (adhesion).  It  is  postulated  here  that 
the  surface  is  such  as  to  maintain  contact  over  the  full  area. 
This  means  the  area  A  will  be  a  simply  connected  region, 
bounded  by  a  simple  closed  curve  called  the  "edge"  of  the 
contact . 

To  insure  a  unique  .solution  to  a  given  contact  problem, 
only  smooth  contacts  will  be  considered.  For  a  "smooth"  con¬ 
tact,  by  definition,  the  displacement  u  and  its  derivatives 
^u  /^x  and  3u  /by  are  continuous  across  the  edge  of  the 

2j  Z 

contact.  This  requires  the  pressure  q(x,y)  to  be  zero  all 
along  the  edge^.  The  alternative  to  smooth  contact  is  a 
sharp  corner  at  the  edge,  with  infinite  pressure  there.  If 
the  rigid  body  had  such  a  sharp  corner,  any  vertical  load  N* 
(above  a  certain  minimum)  would  give  proper  contact .  This 
possibility  is  excluded  by  requiring  smooth  contact. 

^  A  finite  non-zei’o  pressure  at  the  edge  wou'’i  produce  a 
continuous  but  vertical  slope  (see  example  C,  Section  1.6), 
which  would  not  in  general  fit  the  lower  surface  of  the 
rigid  body. 
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The  moving  contact  is  a  mixed  boundary  value  problem, 
since  normal  displacement  is  prescribed  over  part  of  the  sur¬ 
face,  and  normal  stress  is  prescribed  elsewhere.  For  the 
quasi-static,  steady  state  problem,  the  equilibrium  and 
strain-displacement  relations  to  be  satisfied  are 


“ij,j  =  ° 


(3.1a) 


and  the  stress-strain  relations  can  be  written  (from  Section 
1 

00 

=  /  ^3^  Bj^_j(x+VTC,y<z)ac  (3.1b) 


3K^e(x,y,z) 


f  e(x+-vTC.y,z)dC 


(3.1o) 


The  boundary  conditions  on  z  =  0  are 


""  Ty2(x,y,0)  =  0  for  all  x,y 


(3.2a) 


=  0 


outside  A 


(3.2b) 


u  (x,y, 0)  =  6  +  a  X  +  py  +  w*(x-b,y-o)  inside  A 


where  a,p, 6  are  constants  initially  unknown  (a  and  p  are 
components  of  the  possible  small  angle  of  tilt),  and  A  is 
the  initially  unknown  contact  area.  The  unknown  contact 
pressure  inside  A  is  q(x,y)  =  -  a^(x,y,0)  >  0, 


For  smooth  contsct>  u  and  its  derivatives  are  continuous, 

z 

and  q(x,y)  =0  on  the  edge  of  the  contact.  The  load  com- 
ponents,  from  Section  1.5?  are 

N*  =  Jj  q(x,y)dA  ;  F*  =  -  ^  dA 

A  A 

It  is  evident  from  the  mixed  boundary  conditions  and 
the  other  conditions  to  he  met  that  solution  of  a  moving  prob¬ 
lem  will  in  general  be  very  difficult.  It  seems  likely  on 
physical  grounds  that  a  solution  will  exist,  provided  the 
lower  surface  given  by  w*(x,y}  is  sufficiently  smooth. 

This  means,  for  example,  that  for  a  given  load  N'*,  there  will 
be  a  definite  area,  q,  a,  p,  and  6  which  allow  all  the  con¬ 
ditions  to  be  satisfied.  However,  to  attempt  to  find  a  solu¬ 
tion  there  must  be  some  knowledge  or  assumption  about  the 
contact  area.  In  a  general  three-dimensional  problem,  the 
area  is  a  two-dimensional  region,  thus  encompassing  an  infi¬ 
nite  number  of  possible  shapes,  as  well  as  the  other  unknown 
parameters  of  size  and  placement.  The  motion  of  the  base  in¬ 
troduces  an  asymmetry  not  present  in  stationary  contact  pi?ob- 
lems,  so  that  the  methods  of  the  Hertz  theory  for  elastic 
contact  problems  are  not  useful  here.  In  view  of  these 
difficulties,  a  solution  for  even  the  simplest  three- 
dimensional  problems  is  out  of  reach  at  this  time. 
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A  considerable  simplification  results  if  the  problem 
is  two-dimensional,  i.e.,  the  rigid  body  is  an  infinite 
cylinder  (with  arbitrary  cross-section).  Taking  the  cylinder 
axis  to  be  in  the  y-direction,  the  problem  can  be  considered 
in  the  x, z  plane.  The  area  of  contact  is  then  simply  a  line 
segment,  specified  by  two  parameters:  its  length  i,  and  a 
distance  b  fixing  its  placement  on  the  rigid  body  surface. 

The  surface  is  nov/  given  by  z  =  6  +  w*(x-b)  .  For  conven- 
ienc»  (as  in  Chapter  I)  the  contact  wilj.  ta  taken  from  x  =  0 
to  X  =  i.  The  boundary  conditions  (3.2a,b)  become 

=  0  ,  -00  <  X  <  00  (3.3h) 

(3.3b) 

^r,(x,0)  =  6  +  OX  +  w*(x-b)  j  0  <  X  <  £ 


The  unknown  pressure  is  q(x)  =  -g  (x, 0)  for  0  <  x  < 
Other  conditions  to  be  satisfied  are 


q(x)  >  0  ,  0  <  X  <  f  ,-  q(^)  =  q(0)  =  0  (3.3c) 


u„  and  du  /^x 
z  z 


continuous  and  finite 


at  x=0,  x=£ 

(3.3d) 


and  from  (1,45), 


l-v„  , 

- |x|  as  |x|  (3.3e) 


The  displacement  must  be  measured  relative  to  some  point 
near  the  load  (see  Section  1.1).  This  point  will  usually  be 
taken  as  x=0,  z=0. 
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The  load  components  (per  unit  axial  length)  are 


o  o 


Further  discussion  of  moving  contacts  will  be  restricted  to 
two-dimensional  problems. 


3.2  Outline  of  the  T'wq -Dimensional  Solution 

For  a  typical  material  element  at  the  surface,  neither 
stress  nor  displacement  is  knov.'n  for  all  time,  so  the  bound¬ 
ary  conditions  cannot  be  transformed  to  give  an  associated 
elastic  problem.  To  obtain  a  solution,  the  unknown  pressure 
distribution  is  treated  as  if  it  were  known.  The  analysis 
is  carried  out  to  satisfy  the  prescribed  conditions,  which 
then  determine  the  pressure.  With  given  pressure  the  prob¬ 
lem  is  one  of  a  moving  load,  and  the  expressions  for  steady 
state  displacement  given  in  Chapter  I  will  apply.  The 
following  dimensionless  quantities  (most  of  them  introduced 
in  Chapter  I)  will  be  used: 

?=W’  =  W  (3.15a) 


v(?)  s  u^(x)/VT  ;  w(e)  =  w*(x)/VT 


(3.5b) 

(3.5c) 
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N  = 


i  ^  =  -  /  Q(l}a| 


(3.5.d) 


^■''o  IF*  f 

=  —  ?  w  -  -  j 


hv 


6% 


(3.5e) 


Given  properties  of  the  viscoelastic  base  include  7(C) 
(related  to  the  creep  functions,  equation  I.27),  the  time 
parameter  and  velocity  in  the  combination  VT,  and  the  initial 
elastic  response  given  is  the  shape  of  the 

rigid  body  surface  w*(x).  In  an  actual  problem,  the  verti¬ 
cal  load  N*  and  the  tilt  a  might  be  specified,  with  the  ex¬ 
tent  and  placement  of  the  contact  to  be  determined.  However, 
in  the  analysis  it  is  advantageous  to  proceed  otherwise. 
Because  b  enters  the  problem  in  a  complicated  way  in  the 
function  w*,  it  is  much  easier  to  solve  for  a,  which  occurs 
only  as  a  linear  factor  (equation  3.3b).  Thus,  the  extent 
and  placement  of  the  contact  region  are  specified  by  regard¬ 
ing  £  and  b  as  known.  The  tilt  a  and  pressure  Q(|)  are  then 
unknown  quantities.  After  solving  for  these,  the  necessary 
load  components  J  and  P  can  be  determlgned. 

Conditions  to  be  satisfied  are 


Q(e)  >  0  ,  0  <  I  <  >.;  Q(0)  =  Q(?v)  =  0  (3.6a) 

v(4)  and  Bv/3|  continuous  at  ^  =  0  (3.6b) 

and  I  =  “X 
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The  displacement  v(4)  is  measured  from  ^  =  0  ,  so  that 
here  v(|)  represents  what  was  written  v(4)“v(0)  in  Chapter 
I,  Then  v(0)  =  0  ,  so  that  (3.3h)  becomes 

v(|)  =  ai  +  w(4-p)  -  w(-p)  ,  0  <  ^  <  A  (3.7) 

The  relation  between  pressure  and  displacement  is  given 
by  (l.42a,b)  with  =  0*  Thus,  (3.7)  becomes 
A  00 

J  ||f|T|  +  J  7(C)  log|'^-$^!g  r|dc]  d^’  = 

o  o 

a|  +  w(C-p)  -  w(-p)  ,  0<e<x  (3.8) 

This  is  a  Fredholm  integral  equation  of  the  first  kind  for 
the  unknown  pressure  Q(^).  The  solution  gives  Q(^)  in  terms 
of  a,  and  a  is  then  determined  from  the  conditions  of  smooth 
contact  (3.6a).  The  load  components  are  then  given  by  (3.5d,e). 
Once  «(5)  is  known,  the  surface  displacement  outside  the 
contact  region  is  found  from  (1.42). 

To  summarize  the  general  analytic  procedure  for  a  two- 
dimensional  moving  contact  problem;  given  7(C)  ^  ^ii) 

a)  Choose  values  of  A  ,  p 

b)  Solve  the  integral  equation  (3.8)  for  Q(^) 

c)  Determine  a  from  (3.6a) 

d)  Find  the  required  N  and  F  from  (3.5d,e) 
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To  get  numerical  values  of  the  actual  physical  quantities 
[i,  h,  q(x),  N*,  F*],  VT  and  (l-v^)/^^^  must  be  given  also. 
The  integral  equation  (3#8)  can  be  written 


Q(^’)[K(-4‘)  -  K(e-|’)]d|' 
o 


=  ae  +  w(4-p)  ,  0  <  ^  <  A 


where 


00 


K{5)  H  iog|?|  +  [  7(0  iog|l+0  ac 


(3.9) 


The  kernel  K( ' )  is  the  displacement  at  ^  due  to  a 
unit  concentrated  load  at  If  real  material  properties 

are  used,  the  creep  functions  will  be  available  as  curves  dr 
tabulated  data.  Then,  7(C)  is  evaluated  numerically  as  outr 
lined  in  Section  1.3^  and  K(0  would  also  have  to  be  evalii.^ 
ated  numerically.  If  the  material  is  represented  by  a 
mechanical  model  (general  Voigt  solid), 

7(C)  =  S  from  Section  1.3.  Then,  (3.9)  gives 

f  f* 

K(e)  =  (1  +  2  gi)  log  |e|  -  S  E.(-b.O  (3.10) 

The  general  behavior  of  K(^)  is  similar  to  that  shown  in 
Fig.  1.13.  The  function  is  logarithmically  infinite  as 
I  ^  I  -*•  0  and  as  |  ^  |  -♦  00  . 

With  an  analytic  form  for  the  Kernel,  there  is  some  hope 
of  obtaining  an  analytic  solution  to  the  integral  equation 
(3.8)  for  a  given  w(^).  The  stress-strain  laws  for  materials 
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represented  by  mechiaxilcal  models  can  be  given  in  differential 
form  as  well  as  the  integral  form  (3.1b,c).  This  differen¬ 
tial  form  makes  it  possible  to  reduce  certain  viscoelastic 
problems  to  "equivalent  elastic  problems"  (abbreviated  "EEP") 
(not  to  be  confused  with  the  "associated  elastic  problems 
of  Chapter  l) . 

If  the  EEP  can  be  solved,  the  viscoelastic  problem  is 
reduced  to  the  Integration  of  ordinary  differential  equations. 
In  moving  contact  problems,  this  procedure  (called  the  "EEP 
method")  can  be  carried  out  only  for  very  restricted  material 
properties.  In  these  cases,  the  result  is  equivalent  to 
analytic  solution  of  (3.8).  The  results  are  useful  in  spite 
of  the  restrictions,  and  indicate  what  can  be  expected  in 
more  general  cases.  The  next  sections  will  consider  the  EEP 
method  for  steady  state  moving  load  and  moving  contact 
problems. 

3.3  The  Equivalent  Elastic  Problem  for  a  Moving  Load 

(steady  State) 

When  the  stress-strain  laws  can  be  given  in  differential 
form,  such  as  equations  (l.lla,b),  an  EEP  can  be  formed  from 
the  viscoelastic  problem.  The  method  will  be  described  here 
for  the  steady  state  moving  load  problems  of  Chapter  I. 

Since  time  appears  only  in  the  term  x*-Vt  =  x  ,  ^/dt  be¬ 
comes  -V  d/dx  ,  and  the  differential  operators  ( 1.10b) 
become 
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^  \c  ^  Ic 

=  s  Pk(-v)  ^  ^  =  A 'ikt hj 


k=o 


^x 


k=0 


3x 
(3.11) 


*1 

The  stress-strain  laws  (l.lla,b)  become-^ 


PxLS^j(x,y,z)]  =  2\i^  Q^[e^j(x,y,z)]  (3.12a; 


Px[9(x,y,z)]  =  3Kq  Q_J.[e(x,y,z)] 


(3.12b; 


Prom  the  actual  stresses  and  strains  are  obtained 
’’derived  stresses  and  strains,"  denoted  by  and  defined 

ty 

®1J  =  '•  ®  (3.13a: 

'  ®x(®lj^  ’  ®  ”  '^x(®^  (3.-l3b: 


Then,  the  stress-strain  laws  (3.12)  become 


S 


ij 


®  =  3K„  e 


(3.14) 


which  are  the  same  as  the  elastic  stress-^straln  laws.  The 
total  derived  stresses  and  strains  are  defined  by 


€ 


ij  " 


+  J  «1J  ®  =  Px(“lj)  "  5  *lj(VI’x)(®) 


For  convenience,  the  initial  moduli  are  Introduced  as 
Indicated.  Then,  the  P  and  Q  operators  are  all  dlm-;nsionless 
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to  the  equilibrium  equations  gives 


Applying  the  operator  P 


Thus,  a .  .  ,  =  0  in  general  only  if  P„  =  P '  .  Similarly, 

i  J  >  J  A  A 

the  derived  strains  satisfy  the  compatibility  equations  in 
general  only  if 

If  P^  =  P‘  and  Q  =  Q ’  ,  the  viscoelastic  material 
has  identical  behavior  in  shear  and  dilatation,  which  corre- 
spcnds  to  a  constant  Poisson’s  ratio  v.  This  behavior  is  a 
considerable  restriction  on  even  the  idealization  of  a 
mechanical  model,  and  is  not  usually  found  in  real  materials 
such  as  polymers.  However,  by  making  this  restriction  re¬ 
sults  can  be  obtained  approximating  more  realistic  behavior, 
and  indicating  general  features  of  Interest. 

It  will  be  assumed  in  what  follows  that  v  is  constant. 
Then,  ar^j  =  ^ij  ^  derived 

stresses  satisfy  the  equilibrium  equation  ^  =  0  ,  and 

the  derived  strains  are  compatible.  Derived  displacements 
are  defined  by  u.j  =  then  =  l/2(u^  i^  * 

Thus,  the  derived  stresses  and  strains  satisfy  the  same 
relations  as  an  elastic  problem.  The  viscoelastic  problem 
is  reduced  to  its  "equivalent  elastic  problem"  in  the  de¬ 
rived  quantities,  with  boundary  conditions  derived  by  apply¬ 
ing  P^  to  prescribed  tractions  and  Q  to  displacements.  The 

A  X 
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EEP  is  a  problem  in  the  theory  of  elasticity,  and  if  it  can 


be  solved  the  derived  quantities  a^y  etc.,  will  be 
known  explicitly. 

Since  P  and  Q  are  differential  operators,  each  visco- 
elastic  quantity  is  related  by  a  differential  equation  (in 
X,  with  constant  coefficients)  to  the  corresponding  derived 
quantity.  For  example,  ^  =  u^(x). 

Solution  of  these  differential  equations  admits  certain 
arbitrary  functions  of  y  and  z,  which  must  be  determined  to 
complete  the  viscoelastic  solution.  This  requires  consider¬ 
ation  of  stated  or  implied  conditions  of  the  viscoelastic  “ 
solution,  such  as  continuity  or  behavior  at  infinity,  and 
the  boundary  conditions.  Every  viscoelastic  problem  will 
yield  only  one  EEP.  But  a  given  EEP  may  lead  to  many  visco¬ 
elastic  problems,  and  the  additional  conditions  are  needed 
to  distinguish  the  proper  solution. 

It  is  necessary  to  be  very  careful  in  formulating  the 
EEP  from  the  given  viscoelastic  boundary  conditions.  Dis¬ 


continuities  in  the  x-dlrecticn  in  prescribed  boundary 
functions  or  their  derivatives  will  become  singularities  in 
the  EEP  boundary  conditions  after  applying  the  differential 
operators  P  ,  Q  .  These  singularities  must  be  included  in 
obtaining  the  elastic  solution  of  the  EEP.  For  example,  if 
a  prescribed  boundary  traction  has  a  finite  jump,  its  first 
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derivative  becomea  a  concentrated  load  in  the  EEP;,  its  second 
derivative  a  concentrated  coiple,  and  so  on. 

The  EEP  method  of  solving  viscoelastic  pioblems  is  use¬ 
ful  principally  when  it  leads  to  a  closed  form  analytic  solu¬ 
tion.  This  effectively  limits  consideration  to  low  order 
differential  operators,  i.e.,  a  mechanical  model  of  few 
elements.  Otherwise,  the  elastic  solution  of  the  EEP  would 
be  complicated  by  the  presence  of  high  order  singularities, 
and  solution  of  the  differential  equations  would  be  diffi¬ 
cult.  Thus,  in  addition  to  the  requirement  of  constant  v, 
the  material  behavior  is  further  restricted  to  a  model  with 
a  small  number  of  discrete  retardation  times.  Even  with 
these  severe  restrictions,  an  analytic  solution  is  useful 
in  studying  significant  features  of  the  problem  and  identi¬ 
fying  characteristic  quantities.  This  is  particularly  true 
for  moving  contact  problems,  which  will  be  discussed  in  the 
next  seccicxT. 

To  demonstrate  the  method,  a  simple  moving  load  example 
will  be  considered.  The  two-dimer sional  problem,  of  a  constant 
load  q^  on  a  length  i  is  treated  for  the  standard  linear 
solid  (Pig.  1.4).  In  dimensionless  f-  rm.  with 

ft  -  2L.  ^  -  i_ 

^  VT  ^  "  VT 
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the  viscoelastic  operators 
are 

^  ^  ■  It 
^  -  It 

The  load  can  be  written 
q(C)  =  <io[H(e)  -  H(5-X)] 

PIG.  3-2 

The  boundary  condition  is  the i 


q(?)  =  =  qo(l+l')[H(e)-H(|-X)]  -  -  «(5-?')] 


The  discontinuity  in  the  load  gives  rise  to  the  concentrated 
forces  (Pig.  3  '^)  •  For  this  load,  the  elastic  solution  is, 
with  v(^)  =  u..',x)/^?T  , 


^  v(e)  =  (i+f)[£i  -  i  I0g|?|  -  (>>-4)  iog|x-4l]  +  :i.og|4|-iogi4-?'| 

O 

where  A  is  arbitrary.  The  viscoelastic  solution  is  deter¬ 
mined  from 

V  =  Q^(v)  =  V  -  -  e^  (e"^v) 

Integrating  this  with  the  v  given  above  yields 

^  v(4)  =  Ae^  +  A  +  (l+f)[|4-A|log|4-A|  -  4  log|4|] 

+  f[log|4-A|-e*'-‘’'  El{A-4)  -  log|4|  +  e^El(-4)] 
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where  A  is  the  arbitrary  constant  of  integration  As 
v(|)  ~  -  QqA  (1+f)  logjlj  must  be  satisfied.  This  will  be 
true  only  if  A  =  0.  The  final  result  agrees  with  (1.62)  of 
example  c)j  Section  1.6,  with  A  such  that  v(0)  =  0. 


In  this  section  the  EEP  method  is  applied  to  the  two- 
dimensional  steady  state  moving  contact  problem  of  Section  3*2. 
As  before,  the  method  is  restricted  to  materials  with  differ¬ 
ential  streSw's-strain  laws  and  constant  v.  The  differential 
operators  (3.11)  expressed  in  dimensionless  form  are 


P.  =  F.  +  2  pJ-1)  2-^;  Q.  =  1  +  2  qJ-1) 

^  J  k=i  ^  k=i  ^  a?*' 

(3.15) 

Within  the  contact  the  displacement  (3.7)  is  prescribed, 
while  the  pressure  Q(^)  is  unknown. 

Differentiation  does  not  alter  the  contact  region,  so 
the  EEP  is  again  a  contact  problem.  Applying  the  operator 
to  the  prescribed  displacement  (3-7)  gives  the  EEP  boundary 
condition 


v(0  =  Q|[v(4)]  =  v(^)  -  a  +  2  qi^(-l) 


k 


Hid 


.  0  <  4  <  A 

(3.16) 


Solution  of  this  elastic  contact  problem  gives  the  pressure 
Q(|).  Then, 
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(3.17) 


is  the  differential  equation  for  the  actual  viscoelastic 
pressure  e(e). 

Part  of  the  viscoelastic  displacement  is  purely  elastic 
(see  equation  1.58a),  while  the  rest  is  an  integral  of  the 
elastic  displacement.  So  the  viscoelastic  displacement  and 
derivatives  will  have  the  same  (and  no  worse)  discontinuities 
as  the  elastic.  By  limiting  the  problem  to  smooth  contact, 
the  viscoelastic  displacement  and  slope  will  be  continuous. 

p  2 

However,  the  second  derivative  of  the  displacement,  , 

has  a  discontinuity  of  unkn.own  magnitude  at  each  end  of  the 
contact  (^=0  and  i='h)»  and  all  higher  derivatives  have  singu¬ 
larities  there.  This  means  a  of  second  order  (or  higher) 
will  give  singularities  in  displacement  of  unknown  magnitude 
at^  =  0  and  (*=A  in  the  EEP.  Problems  with  such  dis¬ 
placement  singularities  could  be  formulated  mathematically 
in  elasticity  theory,  but  they  would  have  nc  direct  physical 
interpretation.  It  is  unlikely  that  elastic  solutions  to 
such  problems  are  available  or  could  be  readily  obtained. 

Thus,  to  be  able  to  apply  the  EEP  method  with  any  hope 
of  significant  results,  the  material  must  be  further  re¬ 
stricted  to  a  model  with  differential  operators  no  higher 
than  first  order.  Then,  because  of  the  conditions  of  smooth 
contact,  the  displacement  in  the  EEP  is  continuous,  and  the 
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possibility  of  a  concentrated  load  at  the  ends  of  the  contact 
region  is  eliminated.  If  smooth  contact  were  not  required, 
the  EEP  would  have  infinite  jumps  in  displacement  and  con¬ 
centrated  loads  of  Infinite  magnitude  at  each  end  of  the 
contact  region. 

Consideration  is  thus  limited  to  a  Kelvin  solid  (no 
initial  elasticity),  for  which  =  1  ,  or  a  standard  linear 
solid,  for  which  P^  =  1+f  -  ^/S,;  ,  where  1+f  =  initial 
elastic  modulus/final  elastic  modulus.  For  both  solids, 

Q|  =  1  -  .  Even  though  limited  to  these  two  simple 

models,  the  EEP  method  gives  results  that  are  useful.  The 
standard  linear  solid  exhibits  some  significant  characteristics 
of  more  general  viscoelastic  solids,  so  its  behavior  indicates 
what  might  be  expected  in  more  general  moving  contact  problems. 

From  (3.16),  the  displacement  in  the  EEP  is 

v(C)  =  +  w(^-p)  -  w(-p)  -  a  -  ,  0  <  e  <  A 

(3.18) 

and  V  is  continuous  across  the  ends  of  he  contact  region. 

This  is  a  contact  problem  in  the  plane  theory  of  elasticity, 
for  which  general  methods  of  solution  are  well  known.  The 
solution  gives  the  pressure  Q(0  and  the  complete  displace¬ 
ment  v(^)  except  for  an  arbitrary  constant.  Thus,  Q(0 
includes  a  "flat  punch"  term  of  arbitrary  magnitude 
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(3.19) 


ftfp  =  %/V^U^  >  0  <  C  <  I' 

which  gives  constant  displacement  in  the  contact  region. 

For  the  Kelvin  solid,  Q(4)  =  Q(0  since  =  1.  For 
the  standard  linear  solid,  from  (3.17) 

a{4)  =  (1  +  f)  Q  -  II 

and  thus  ^ 

ft(4)  =  H-  f  ft(C)dC  (3.20) 

Where  A  and  d  are  arbitrary  constants.  The  constants,  Q^, 
and  the  unknown  tilt  a  ,  are  determined  by  the  smooth  con¬ 
tact  conditions  (3.6a,b)  and  the  asymptotic  form  of  the  dis¬ 
placement  (from  equation  3.3e) 

v(^) - (1+f)  N  logj^l  as  1^1  -♦oo  (3.21) 

Once  q’(^)  and  a  have  been  found,  the  solution  is  essentially 
complete,  and  other  quantities  of  ■'.nterest  can  be  determined. 
Application  of  this  method  to  a  particular  problem,  the  roll¬ 
ing  cylinder,,  is  given  in  the  next  chapter. 
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CHAPTER  IV 


THE  ROLLING  CYLINDER  PROBLEM 

One  of  the  principal  problems  of  moving  contact  in 
viscoelasticity  is  the  rolling  of  a  circular  cylinder  (or 
sphere)  on  a  flat  surface.  Osborne  Reynolds  made  an  exten¬ 
sive  study  of  the  problem  in  I876,  and  since  then  several 
other  attempts  to  explain  and  predict  rolling  resistance 
have  been  made^.  The  most  significant  recent  developments 
have  resulted  from  the  work  of  Tabor  and  his  colleagues 
([^  [16],  C30].  [31]^  [32])  since  about  1950.  The  ex¬ 

perimental  work  showed  that  rolling  resistance  is  nearly 
Independent  of  surface  lubrication^  and  arises  principally 
from  energy  dissipation  within  the  rolled  material.  This 
has  led  to  two  general  theories  attempting  to  explain  roll¬ 
ing  resistance.  One,  developed  chiefly  by  the  Tabor  group, 
makes  use  of  measurements  of  hysteresis  loss  in  simple 
loading -unloading  cycles,  together  with  the  solution  for 
rolling  on  an  elastic  base.  A  similar  procedure  has  been 
suggested  by  Drutowski  [ 7  ] .  A  second  possibility  is  to 
consider  materials  which  can,  at  least  approximately,  be 
regarded  as  linearly  viscoelastic.  As  already  explained  in 
Section  1.5j  viscoelastic  materials  dissipate  energy, which 

I  ^ 

Discussion  and  additional  references  in  Bowden  and 
Tabor  [4],  Kelly  [18],  Drutowski  [6]. 
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gives  rise  to  resistance  to  moving  loads.  Flom  ([10],  [11  ]> 
[12])  has  sought  to  connect  rolling  friction  with  simple 
energy  loss  measurements  on  viscoelastic  materials.  His 
tests  of  hard  spheres  rolling  on  several  viscoelastic  poly¬ 
mers  show  that  the  variation  of  coefficient  of  friction  with 
velocity  of  rolling  is  qualitatively  similar  to  the  descrip¬ 
tion  of  Section  1.5  (Fig.  9 )♦  There  are  few  other  ex¬ 
perimental  results  available  for  spheres  or  cylinders  roll¬ 
ing  on  materials  of  viscoelastic  character.  Some  tests  of 
Tabor  [l6]  on  rubber  indicate  that  for  cylinders  longer  than 
several  times  the  diameter,  the  behavior  is  essentially  two- 
dimensional. 

This  chapter  considers  in  detail  the  two-dimensional 
problem  of  a  rigid  circular  cylinder  rolling  on  a  semi¬ 
infinite  viscoelastic  solid.  The  contact  is  supposed  per¬ 
fectly  lubricated,  so  there  are  no  tangential  forces  at  the 
contact  surface  and  there  is  no  distinction  between  rolling 
and  sliding.  This  is  treated  as  an  example  of  the  general 
moving  contact  problem  of  Chapter  III.  Several  published 
treatments  of  the  problem  are  examined.  Each  involves  spe¬ 
cific  limitations  which  prevent  more  general  application. 

A  different,  less  limited,  method  is  proposed  and  illustrated 
in  the  next  chapter. 


Ill 


4.1  Formulation  of  the  Rolling  Cylinder  Problem 


The  rigid  cylinder,  of  radius  R,  is  considered  fixed, 
and  the  viscoelastic  base  moves  under  it  with  velocity  V. 

The  contact  region  extends 
from  X  =  0  to  X  =  i,  and  the 
downstream  end  x  =  0  is  a 
distance  b  from  the  center  of 
the  cylinder  (Fig.  4.1).  De¬ 
formations  are  assumed  small 
(i.e.,  i  «  R),  so  that  within 
the  contact  region  the  surface 

PIG.  4.1 

of  the  cylinder  can  be  approxi¬ 
mated  by  a  parabola^: 

w(e-p) 

In  this  problem  w(^-p)  contains  a  term  linear  in  ^  with 
unknown  magnitude  p,  so  there  is  no  need  to  Introduce  the 
additional  tilt  a.  But  now  B  is  to  be  found,  and  not, 
as  in  Section  3-2,  assigned  initially.  For  a  given  i,  vari¬ 
ation  of  b  in  this  special  case  amounts  to  a  tilt  of  the 
conta  ct  surf a  ce . 

The  displacement  within  the  contact  region  is,  from  (3-7) 
v{?)  =  I  (PI  -  il^)  ,  0  <  I  <  A 

^  Using  quantities  defined  by  (3«5a-e)  and  p  =  R/VT. 
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It  is  convenient  to  introduce  new  dimensionless  quantities 


qU)  =  p  Q(?)  = 

0  < 

0 

—  q(x) 

TT  ' 

(4.2a) 

n 

l-v^  0 

»v*  r' 

(4.2b) 

N'  =  2p  N  =  ^ 

**0 

^  =  2  j  Q'(^)de 

The  integral  equation  (3-8)  then  becomes 


f  Q’(^’)[K(-^')  -  KU-V)]6V  =  - 1  r  >  0  <  e  <  X 

°  (4.3) 


Proceeding  as  in  Chapter  III,  the  creep  function  7(C) 
is  given,  and  a  value  of  the  contact  length  A  is  chosen. 

The  solution  of  (4.3),  along  with  the  smooth  contact  condi¬ 
tions  (3.6a,b),  gives  the  pressure  Q'(^)  and  the  parameter 

o 

p.  The  total  load  N’  is  then  found  from  (4.2b).  Rolling 
resistance  is  conveniently  represented  by  the  coefficient 
of  friction  x*  Using  (3.5d,e),  (4.1),  and  (4.2a,b), 


The  actual  physical  and  geometrical  quantities  can  be 
conveniently  represented  by  dimensionless  combinations  in¬ 
volving  only  the  radius  R,  actual  total  load  N*,  and  initial 
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elastic  constants  These  combinations,  given  on 
the  left-hand  sides  of  (4.5a-d),  are  formed  as  follows  (the 
quantities  on  the  right-hand  sides  having  been  found  in  the 
solution  for  a  given  Tv)  : 


velocity 

VT 

^o 

) 

(4.5a) 

pressure 

^o 

q{x) 

F  ^’(e) 

^o 

{4.5b) 

Coefficient  of  friction 

^0  ^ 

VT 

a“ 

0 

(4.5o) 

placement 

where 

f  1-v^ 

L_ 

a 

0 

P 

1 

i2 

VT 

—  B 
% 

(4.5d) 

^o  = 

1.; 

-  R  N' 

TT 

*] 

(4.6) 

is  one-half  the  contact 

length 

for  an 

elastic  solid. 

It  is 

then  possible  to  exhibit  the  variation  with  velocity  of  q(x), 
X,  or  b,  for  a  given  total  load  and  radius.  For  example,  a 
plot  of  Rx/Sq  ve  VT/a^  shows  the  variation  of  coef¬ 
ficient  of  friction  (rolling  resistance)  with  velocity. 

In  Section  1.4  it  was  shown  that  the  limiting  cases 
V  0  (?v  -►  oo)  and  V  -►  »  (Tv  -*•  0)  are  elastic.  This  can  be 
used  to  find  the  corresponding  limiting  values  of  certain 
quantities  in  the  rolling  cylinder  problem.  For  V  -*■  w  , 
the  material  exhibits  its  initial  elastic  behavior.  The 
solution  for  given  displacement  (4.1)  is 
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b  =  I  ;  q(x)  =  j:~-  ^  ./x{£~: 


.  N*  =  ■—  ^ 

"  "  B  1-v,  R 


For  V  0,  the  material  exhibits  its  final  elastic  behavior. 

The  solution  is  the  same  as  for  V ,  except  that  I^o/(1“Vq) 

is  replaced  by  m.^/(  1-v^)  =  l/(  1+f )  •  [x^/(  1-v^)  ,  where 
00 

f'  =  /  7(C) <3^  .  Some  limiting  values  are  given  in  Table  4.1. 


2  ^  ^o  R  N* 
TT  i  i 


PX 

O  O 


Table  4.1 


Limiting  Value 


?v  -*■  0 
V 


£  VT  , 


2  v/I+f 


Q’  ^  ''o  1  R  ,  . 

A  =  HT-  ?  j  1 


=  .15915. 


1 

5?  ^+f 


n«  VT  ^O  f  . 

Q  T-  =  ^  q  (niax) 


i  =  .31831. 


1  1 
^x/I+f 
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4.2  Analytic  Solution  by  the  ESP  Method 


The  rolling  cylinder  problem  can  be  solved  by  use  of 
the  EEP  method  outlined  in  Section  3.3*  The  viscoelastic 
material  must  be  a  standard  linear  (or  Kelvin)  solid  with 
constant  v.  This  is  a  severe  restriction  on  the  material 
behavior,  but  it  allows  a  complete  analytic  solution.  The 
same  solution  has  been -previously  obtained  by  Hunter  [I?]  . 

The  EEP  method  makes  more  evident  the  underlying  assumptions 
and  Conditions  of  the  problem,  and  is  clearer  and  more  direct 
than  the  somewhat  different  procedure  of  Hunter.  The  final 
results  of  the  EEP  method  are  the  same  as  those  of  Hunter. 

Details  of  the  procedure  for  a  standard  linear  solid 
are  now  given.  The  differential  operators  are 

‘’e  =  ^  ^  ■  Ir  ’  ^  ■  5f 

Applying  to  the  displacement  (4.1)  gives 

p  v(e)  =  D  +  (1+p)  0<e<A  (4.7) 

where  D  is  an  arbitrary  constant.  The  elastic  contact  prob¬ 
lem  with  prescribed  displacement  (4.7)  has  a  pressure^ 

A+B(e-^)  _ 

Q(e)  =  —  .  0  <  e  <  X  QJ 

^  Muskhelishvili  [25],  Section  ll6a. 


Il6 


Frcm  (3-20)  and 


where  A  is  arbitrary  and  B  =  i  +  g  ^  . 

(4.2a),  the  viscoelastic  solution  is 

Q’(£)  =  “  f  \ - - ^- +  VCCX-O  ]<3C 

^  (^.9) 

The  upper  limit  on  the  integral  is  here  chosen  so  the  con¬ 
dition  Q*(^)  =0  is  satisfied.  The  condition  Q‘(0)  =-.  0 
must  also  he  satisfied.  Setting  ^  0  in  (4.9),  the  inte¬ 

gral  can  be  evaluated^.  It  must  vanish  identically,  giving 
the  relation 

=  2  I^(h)  "  I+f  ^  (4.-10) 

where  It ,  I^  are  modified  Bessel  functions  of  the  first 
1  o 

kind,  and 

h  =  (1  f)  I  {‘i.lX) 

Another  relation  is  needed  between  A  and  B.  To  obtain 
this,  the  displacement  must  be  considered.  Outside  the  con¬ 
tact  region,  the  elastic  displacement  due  to  the  pressure 
(4.8)  is 

p  v(e)  -  ,D  -p  e®  -(a  +  ^)oosh"^  l|i  .  i| 

-  Bi|i  -  i  15  -  |l]  v^frF>T,  5<0  or  5  >  Tv 

_  ('>•12) 

^  Erdelyi  [8],  equation  7.12(10). 
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where  C  is  an  arbitrary  constant  (which  may  have  different 
values  inside  and  outside  the  contact  region)  and  the  upper 
Ij.mit  is  chosen  for  convenience.  As  [^|  ,  this  expres¬ 

sion  gives 

pv(^)  ~  pv(0  +  c^e^  =  C^e^  _  (A  +  ^)  log  |^|  +  constant  + 

This  will  have  the  proper  asymptotic  form  (3.21)  only  if 
=  0.  Then,  substituting  for  pv  from  (4.7)  or  (4.12), 
and  choosing  D  =  -  p  gives 

pv(e)  =  PI  -  I  {2  _  ^  I  <  0  or  e  >  A  (4.133) 

where 

00  2 

I(^)  =  f  [(A  +  3-)  cosh"^  ||^  -  1| 

i 

+  (B  i|i  -  I  |C  -  ||)  x/CTFXT  ]e-^dC  (4.13b) 
(for  i  <  0,  the  upper  limit  is  0  rather  than  00) 

Inside  the  contact  region,  the  displacement  is 

pv(e)  =  -  I  +  [Cg  -  1(A) ]e^ 
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ChccsirLG!  o:lveB  t‘h<=*  cOT’r*pp.-h  <?n1pr*Ampnr 

To  meet  the  condition  of  continuous  displacement  at  4  = 
requires j  on  comparing  (^.l)  and  (4.13a),  that  l(?v)  =0. 
Evaluation  of  the  integral  l(?v)  leads  to  the  relation 


A  =  — 
^  2 


X  K^(X/2) 


(1“B) 


(4.1 


where  K^,  are  modified  Bessel  functions  of  the  second 
kind.  Prom  the  two  relations  (4.10)  anc  (4.l4)  A  and  B 
can  be  determined,  and  then  p  =  B  -  1  +  ^  . 

Continuity  of  v(|)  at  i  =  7^  has  been  guaranteed,  and 
from  (4.13a,b)  this  is  true  at  ^  =  C  also.  Since 
V  =  V  -  and  both  v  and  v  are  continuous,,  the  slope 

is  also  continuous.  From  (3.5d),  and  since 

« =  f|(a)  =  (i+f)  « -  II 

the  total  load  is  given  by 

A  ?v 

(i+f)  N  =  (i+f)  J  Q(?)ae  =  J ft(e)  a?  +  (}(7v)  -  qCo) 

o  o 

Thus,  using  (4.8),  with  Q(7v)  =  Q(0)  =  0  , 


(l+f)  pN  =  i  J  [ 


?  rA+B(e  -  h 


d|  =  A  4-  ^ 


Erdelyi  [8],  equation  7.3(15) 


1] 


The  asymptotic  form  of  the  displacement  (4.13a,b)  as  j^j  -*■  « 

v(|)  ~  +  ^)  log  |5|  =  -  (1+f)  N  log  |?1 

(4.15) 

whli'h  satisfies  the  condition  (3.21).  All  the  conditions 
of  the  problem  are  now  satisfied. 

This  completes  the  solution  for  a  rolling  cylinder  on 
a  standard  linear  solid  base.  Numerical  results  are  obtained 
most  directly  by  specifying  \  (with  a  given  value  of  f), 
then  proceeding  with  evaluation,  r^'lng  the  following  formulas: 


1) 

h  =  1  (1+f) 

(4.11) 

2) 

Ko(V2)I,(h) 

[0  <  m  <  1] 

(4.16) 

3) 

,  4  f  '^l'’'/^) 

*  “  2  Itf  K  (X72) 

m 

1+m 

(4.17a) 

4) 

5) 

6) 

7) 


B  =  1  - 


m 


1+f  1+m 


ft  _  A  -JL  JL. 

P  "  2  "  1+f  1+m 


[p<|] 


«'(0  =  j  e 


l^(l+f)TiV2 


J 


■^i-t^ 


where  t  =  2^/X-l 
2 


N'  =  2p  N  =  ^ 


A  +  -g- 


(4.17b) 

(4.18) 

+  ]<2t 

(4.19) 

(4.21) 
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Evaluation  of  the  coefficient  of  friction  from  (4.4)  leads 


to 


Bfl  + 


81 


-1 


(4.21) 


The  surface  displacement  for  the  rolling  cintact  is  given  by 

pv(|)  =  @5  -  I  e®  -  I  J  3(t)  *  dt  ,  5  <  0 

^  (4.22a) 

=  0  <\  .(4.22b) 

=  PI  -  I  -  I  J  g(t)  ^  dt  ,  i  >-K 

^  (4.22c.) 

where 

2 

g(t)  =  (a  +  g-^cosh  ^jt  j  +  ^  (B  -  ^  t)  Vt ' £ 

( 4 . 23 ) 


To  obtain  the  pressure  distribution  over  the  contact  region 
a  numerical  evaluation  in  step  5)  is  necessary.  Hunter  [I7] 
gives  a  plot  of  Q’  (|)  for  the  one  case  f  =  1,  Tv  =  1.6.  A 
sketch  of  the  general  case  is  given  in  Fig.  4.2.,  along  with 
an  elastic  contact  for  comparison.  The  effect  of  viscoelas¬ 
ticity  is  seen  in  “^he  non-symr.  of  contact  region  and 
pressure  distribution,  and  the  resultant  horizontal  force. 

The  variation  of  rolling  resistance  with  velocity  is 
best  presented  by  plotting  Rx/^q  vs  VT/a^  ,  as  suggested 
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in  t/iis  pirsccding  ssction.  r*Oi'  s  fixeu  soiU'-'j.ons  lOi* 
different  velocities  are  obtained  by  varying  "h.  The  result 
is  a  curve  of  the  form  of  Fig.  1.9,  with  a  single  maximum. 
Typical  curves  are  given  by  Hunter*  for  f=landf=9* 

In  an  actual  rolling  cylinder  problem^  the  physical 
quantities  f,V,TjN*,Rj  might  be  specified.  To  determine  the 

fi 

contact  region  and  pressure,  a  plot  of  A  vs  VT/a^  (for  a 

given  f)  is  useful.  The  given  data  determines  VT/a  ,  and 

<*«/ 

the  appropriate  A  can  be  found  from  the  plot.  Then  the 
other  quantities  can  be  determined  by  the  procedure  just 
given. 


Elastic  Viscoelastic 


Note  that  Hunter’s  a  is  i/l+f  times  the 
here.  ° 


a^  used 
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same  way  as  for  the  standard  linear  solid.  There  is  no 
initial  elasticitj^  so  the  final  elastic  response  is  used. 
For  this  material,  =  1  and  =  1  -  ^/S|  .  Only  a 
need  be  specified  to  obtain  a  solution.  The  results  are 
summarized  below: 


=  +  K^(X/2) 


(4.24) 


P  =  B  +  $  -  1 


(4.25) 


p  q(e)  =  b/^  + 


(4.26) 


'f  1  N*  -h 


TT  VT  2 


=  2  (®  f) 


(4.27) 


A  6b+A  _  33+3  “A  „  /jt  oQ^ 

P-  “  1.  4b¥a  "  3  =  a  4p+4-A  "  ^  ^ 


The  Kelvin  solid  is  of  interest  as  the  limiting  case  of  a 
standard  linear  solid  with  the  initial  elastic  response  much 
smaller  than  the  final  response  (i.e.,  very  large  f ) .  In 
this  limiting  case,  the  pressure  is  infinite  at  the  upstream 
end,  I  =  A  (see  Fig.  4.2). 


4-3  DjLrsct-  Anslytlcsl  Solution  (Morlsnd) 

♦ 

Starting  from  the  fundamental  field  equations  (3.1a) 
and  the  general  stress-strain  law  in  integral  form  (1.9a,b), 
viscoelastic  problems  can  be  treated  directly  as  boundary 
value  (space)  and  initial  value  (time)  problems.  Norland 
[24]  has  investigated  the  rolling  cylinder  problem  in  this 
way.  Conditions  to  be  satisfied  are  given  by  (3.3a-e).  The 
general  procedure  is  outlined  be\cv,  and  some  d eta i 3, s., of 
Norland’s  solution  are  given. 

At  first,  the  analysis  applies  to  the  general  two- 
dimensional  problem  of  a  normal  load  moving  on  the  surface 
of  a  semi-infinite  base  (as  in  Chapter  I).  Norland  takes 
a  Fourier  transform  with  respect  to  time.  The  transformed 
problem  is  then  essentially  an  elastic  problem,  and  the 
solution  is  expressed  as  Fourier  integrals  (in  x*) .  At  this 
stage  the  condition  of  steady  state  is  introduced,  and  after 
considerable  manipulation  the  time  variable  is  replaced  by 
the  steady  state  variable  x  =  x*  -  Vt.  The  result  is  then^ 
(in  the  notation  of  equation  1.29) 

00 

Q(4)  ^  f  i^^‘i(s)  cos  s^  +  f2(s)  sin  s^]  ds  (4.29a) 

<J 

•  o 


The  v;hole  analysis  to  this  point  could  be  done  more 
directly  and  simply  by. considering  a  moving  load,  steady 
state  problem  from  the  beginning,  including  the  appropriate 
stress-strain  laws  (3.1b,c).  Then,  a  Fourier  transform  with 
respect  to  ^(=x/VT)  gives  an  elastic  problem  which  leads 
directly  to  equations (4. 29a,b) . 
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00 

v(^)  =  J  0(s)  |[f3^(s)  +  ^(s)f2(s)](  cos  -  1) 
o 

+  [fgCs)  -  ^(s)f^(s)]  sin  sc|  J 


(4.29b) 


where  f^Cs)  =  fi(-s)  and  ^2^®)  ^  funci:ions 

0(s)  and  ^(s)  are  related  to  the  creep  functions  and  are 
defined  by 


0(s)[l  -  i^(s)]  =  1  +  7(is)  (4.30) 


where  7(s)  is  the  Laplace  transform  of  7(C)  and 
i  =-v/^  .  For  a  given  moving  load,  f^  and  f^  are  known,  and. 
(4.29b)  gives  the  resulting  deflection. 

The  moving  contact  problem  is  considered,  as  before, 
by  assuming  some  length  of  contact,  and  finding  the  neces¬ 
sary  total  load.  Then,  and  f2(s)  are  unknown  functions. 

Once  they  have  been  found,  the  complete  solution  can  be 
obtained.  It  is  convenient 
here  to  take  ^  =  0  at  the 
center  of  the  contact  region. 

This  region  is  then  j^j  <  7^/2 
(see  Fig.  4.3).  The  boundary 
conditions  (3-3)  become 


:.25 


FIG.  4.3 


(4.313) 


Q(?)  =  0  ,  III  >  I 
v(l)  =  i  [(P  -  5)  I  -  I  1^3  ,  |e!<|  (4.31b) 

It  is  necessary  to  add  an  additional  flat  punch  pressure, 

Qo[(V2)^  “  to  Q(4)  given  by  (4,29a).  The  re¬ 

sulting  displacement  is  found  by  the  method  of  Chapter  I, 
and  is  added  to  (4.29b).  For  this,  Morland  assumes 
7(0  =  2  f^e"^^^ — i.e.,  a  mechanical  model — for  the  visco¬ 
elastic  behavior.  The  additional  displacement  is  then 

V„(l)  =  ’MJo  f  ^  >  1^1  s  I  ('‘•32) 

with  this  assumption  for  y{C)> 

if*  b  f* 

♦(s)  =  1  +  2  ,  i,(s).^(s)  =  S  2 -gi-s  (4.33) 

i  b^+s*^  i  b^+s 

Satisfying  the  boundary  conditions  (4.31a,b)  leads  to  the 
equations 

00 

J  [^^^(s)  cos  s%  +  f2(s)  sin  s|]  ds  =  0  ,  |C|  >  |  (4.34a) 
o 

J  «[(f;^+Kg)(cos  s|-l)  +  (f2-V.f3^)slr.  s|]|&  +  irQg  2  ^  KQ(bj^  |)(l-e‘’l5) 

O  i  i  ■ 

=  |((P  -  |)  I  -  I  1^]  ,  |l|>|  (4.34b) 
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These  are  dual  integral  equations  for  the  unknown  functions 
f^(s),  f^Cs).  Norland  procedes  by  splitting  each  of  (4.34a,b) 
into  odd  and  even  parts.  Then,  a  solution  is  assumed  in  the 
form  of  series  of  Bessel  functions 


fl(s)  =  ttQ^ 


00 

2 

m=l 


2m  "2m^2 


(|s) 


00 

o^^T^2m-l‘^2m- 
m=l 


(4  s) 


1^2 


which  satisfy  (4.34a)  identically.  Introducing  these  into 
the  odd  and  even  parts  of  (4.34b)  leads  eventually  to  two 
Infinite  sets  of  equations,  which  are  combined  into  one 
infinite  set  of  linear  equations  for  the  unknown  constants 
a  .  The  coefficients  of  in  the  set  are  integrals  of 
<t>  and  Tp  with  Bessel  functions,  a  typical  integral  being 

00 

/ 

o 

The  solution  of  this  infinite  set  of  equations,  with  the 
additional  conditions  ;i(+  A/2)  =  0  ,  gives  the  a^,  Q^, 
and  B.  This  formally  solves  the  contact  problem,  giving 
f^  and  f2  as  an  infinite  series.  The  contact  pressure  is 
then  given  by 

=  I  II  ^  \  ®2m  +  Van-1  ] 

m=i-  m=l 
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where  |  =  X/2  sin  9.  From  this  it  follows  that 

N  =  ,  PX  =  I  [1  -  +  I  ] 

Thus,  for  a  given  >v,  the  coefficients  of  a^^^  must  be  determined 
and  an  infinite  set  of  linear  equations  solved.  Then  the 
pressure,  total  load,  friction  coefficient  and  other  quantities 
of  interest  can  be  found.  The  whole  process  must  be  repeated 
for  each  value  of 

To  carry  out  the  calculations,  Norland  makes  some  simpli¬ 
fying  approximations.  The  material  behavior  is  approximated 
by  assuming  0  =  constant  and  f  =  constant  in  the  integral 
pare  of  (4.34b).  The  coefficient  integrals  can  then  be 
evaluated  in  closed  form.  For  the  "flat  punch"  displacement 
(4.32),  a  different  approximation  is  used  by  taking  a  single 
term  in  the  sum  (i.e.,  standard  linear  solid  with  constant  v) . 

The  solution  is  further  approximated  by  using  a  ^x^  block 
Instead  of  the  infir-ite  set  of  equations.  Even  with  these  ajproxl- 
rnations,  much  numerical  work  is  needed  to  get  results  for 
a  single  value  of  Tv.  Norland  gives  one  numerical  example 
and  the  resulting  curve  of  pressure  distribution  over  the 
contact  region.  The  results  are  given  below  for  this  exam¬ 
ple,  in  which  Tv  =  1,  f  =  ,5*  Also  given  are  the  corres¬ 
ponding  results  by  the  method  of  Section  4.2  (Hunter)  fer 
comparison. 
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Quantity 


Method 


Hunter 

Norland 

VT 

a 

2.40 

2.21 

^0 

.446 

.490 

PX 

.067 

.072 

4.4  Other  Treatments  of  Rolling  Contact  Problems;  Comments 

Although  each  of  the  methods  of  Sections  4.2  and  4.3 
can  be  applied  only  with  certain  restrictions  and  approxi¬ 
mations,  each  starts  from  an  ’’exact"  formulation  of  the  roll-: 
Ing  contact  problem  In  the  sense  that  the  viscoelastic  base 
Is  treated  as  a  continuum,  and  the  mutual  effect  of  adjacent 
elements  Is  taken  Into  account.  However,  only  two-dimensional 
problems  can  be  satisfactorily  treated,  and  analytic  solutions 
are  possible  only  for  limited  kinds  of  viscoelastic  behavior. 

Rolling  contact  problems  have  also  been  discussed  by 
adopting  a  representation  of  material  behavior  which  is  more 
widely  applicable,  but  much  more  approximate.  The  visco¬ 
elastic  base  is  assumed  to  be  made  up  of  Independent  verti¬ 
cal  columns,  like  a  Winkler  foundation.  Each  colifinn  is  a 
one -dimensional  viscoelastic  rod  with  regard  to  vertical 
stress  and  strain,  and  the  deformation  of  any  one  column 
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has  no  effect  on  and  is  not  affected  by  any  other  columns. 
Thus,  some  of  the  features  of  viscoelasticity  are  present, 
including  delayed  recovery  and  energy  dissipation,  but  the 
overall  behavior  of  the  base  as  a  continuous  material  is 
only  approximately  represented. 

With  this  simplified  representation,  both  three  and 
two-dimensional  rolling  contact  problems  can  be  treated. 

The  base  has  a  finite  thickness  h,  and  the  vertical  dis¬ 
placement  of  any  column  is  thus  he  ,  where  e  is  the 
strain.  The  stress-strain  law  can  be  written 

u 

u 

a(t)  -  ^  J  G(t-T)  e(T)  dT  (^.35) 

o 

where  G(t)  is  the  relaxation  function.  If  the  surface  dis¬ 
placement  is  prescribed,  (4.35)  gives  the  necessary  contact 
pressure.  By  prescribing  a  spherical  or  cylindrical  dis¬ 
placement,  the  rolling  contact  problem  is  essentially  solved. 
Due  to  the  delayed  recovery  behln'^  the  roller,  the  contact 
region  and  pressure  are  not  syrari'-tric.  This  gives  a  resist¬ 
ing  force  and  a  corresponding  coefficient  of  friction,  and 
thus  provides  some  qualitative  description  of  rolling 
resistance. 

Applications  of  this  simplified  treatment  of  rolling 
contact  problems  Include:  Flom  and  Bueche  [13],  sphere  on 
a  Kelvin  solid;  May,  Morris  and  Atack  [23]^  cylinder  on 
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Maxwell  fluid,  standard  linear  and  more  general  solids; 

Norman  [26],  cylinder  on  a  Kelvin  solid.  Results  are  quali¬ 
tatively  similar  to  those  of  the  more  “exact”  methods  of 
Hunter  and  Morland.  In  particular,  the  coefficient  of 
friction  varies  with  velocity  in  the  manner  previously  dis¬ 
cussed,  increasing  with  velocity  to  a  single  maximum,  then 
decreasing  to  zero.  When  a  mechanical  model  is  used  for  the 
material  behavior,  G(t)  is  a  sum  of  negative  exponentials. 
Then  the  results  can  be  expressed  in  analytic  form,  in 
terms  of  elementary  functions. 

In  summary,  this  simplified  treatment  gives  some  use¬ 
ful  results  of  a  qualitative  nature  for  the  rolling  contact 
problem,  and  involves  only  straightforward  (although  perhaps 
tedious)  integrations.  However,  it  is  a  somewhat  crude 
representation  of  the  actual  behavior  of  the  base.  For  the 
more  accurate  representation  of  the  base  as  a  continuum, 
analytic  solutions  (such  as  Hunter,  Section  4.2,  and  Norland,, 
Section  ^.3)  are  available  only  in  limited  cases  of  two- 
dimensional  problems.  The  direct  method  (Morland)  applies 
in  principle  for  general  viscoelastic  behavior,  but  in 
practice  approximations  are  needed.  Even  then,  an  infinite 
set  of  equations  must  be  solved,  which  requires  further 
approximation  and  much  numerical  work  for  each  particular 
example.  The  EEP  method  gives  a  complete  closed -form  ana¬ 
lytic  solution,  with  no  approximation  in  the  analysis,  and 
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numerical  work  is  needed  only  to  get  the  detailed  pressure 
distribution.  However,  this  method  can  be  used  only  for 
simple  models  (standard  linear  or  Kelvin  solids,  with  con¬ 
stant  v)  that  in  some  respects  represent  the  behavior  of 
viscoelastic  materials.  Thus  there  are  significant  limita¬ 
tions  to  each  of  these  methods  of  obtaining  an  analytic  solu¬ 
tion  of  the  rolling  contact  problem. 
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CHAPTER  y 


A  NUMERICAL  METHOD  FOR  TWO -DII4ENSIONAL 
MOVING  CONTACT  PROBLEMS 

There  is  no  complete  analytic  solution  of  the  moving 
contact  problem  available  for  a  general  viscoelastic  material. 
The  EEP  method  discussed  in  Chapters  III  and  IV  can  be  used 
only  for  a  material  of  very  restricted  type,  represented  by 
a  simple  model  with  constant  v.  For  a  rolling  cylinder  this 
method  allows  a  complete  analytic  solution,  but  for  other 
shapes  of  contact  surface  the  analysis  is  much  more  compli¬ 
cated  and  does  not  in  general  yield  closed-form  expressions. 
The  direct  analytical  method  discussed  in  Section  4.3  can  in 
principle  be  applied  to  any  contact  surface  on  any  viscoelas¬ 
tic  material.  But  involved  analysis  is  required,  and  then 
much  approximation  and  numerical  work  is  needed  to  obtain,  a 
definite  result. 

For  the  general  moving  contact  problem,  with  arbitrary 
contact  surface  on  any  viscoelastic  solid,  it  is  evident 
that  the  solution  murt  be  obtained  numerically.  This  can 
best  be  achieved  by  using  numerical  procedures  from  the  be¬ 
ginning,  thus  avoiding  lengthy  •'=»nd  elaborate  analytical 
procedures.  Creep  functions  available  only  as  numerical 
data  can  be  used  directly,  eliminating  the  laborious  task 
of  trying  to  represent  the  data  analytically.  This  chapter 
presents  a  numerical  method  for  solving  the  two-dimensional 
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steady  state  moving  contact  problem  (as  formulated  in  Chapter 
III).  The  procedure^  restricted  only  by  the  requirement  of 
smooth  contact,  is  simple  and  direct,  and  readily  adapted 
for  digital  computer.  The  method  is  illustrated  with  the 
specific  problems  of  a  rolling  cylinder  and  a  "nearly  flat 
punch " . 

5.1  Numerical  Evaluation  of  Surface  Displacement 

The  steady  state  surface  displacement  due  to  a  two- 
dimensional  pressure  distribution  Q(^)  (0  <  ^  <  A)  moving 
on  a  semi-infinite  viscoelastic  base,  given  by  (l.42a,b), 
is  ^ 

v(5)  -  vcig)  =  J a(e')tK(5.-e')-K(e-e')id?'  {s.i) 

o 

where 

00 

K{?)  =  log  i?|  +  J  7(0  log  (?+0  dC  {5.2') 

o 

and  is  an  arbitrary  point  on  the  surface  (Fig.  5.1).  The 
dimensionless  quantities  given  by  (3.5a-e)  are  used.  It  is 
supposed  that  the  pressure  is  everywhere  finite,  and  that 
Q(0)  =  Q(^)  =  0.  This  is  the  case,  for  example,  when  Q(|) 
arises  from  a  "smooth"  moving  contact  (definecii  in  Section  3.1). 

To  evaluate  numerically  the  integral  in  (5.1),  a  smooth 
pressure  distribution  Q(4)  on  0  <  ^  <  A  is  replaced  by  n 
trapezoidal  load  elements  (Pig.  5.1).  Within  a  typical 
segment  ^  =  1,2,  ...,n)  the  pressure  is  repre¬ 

sented  by 
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PIG.  5.1 


Q(^)  =  (i)  - 

where  Qj^  =  ^  ^o  ^  ^  ^  ^n  ^  ''*  displacement  (5-i) 

then  becomes 

n  5*^ 

V(«)  -  vClg)  =  2  /  Q‘‘(«')[K(e3-e')-K(|-4')]at' 

Ic—  X 

«k-l 

Substituting  for  OiH^)  from  (5.3)  ant'  for  K(4)  from  (5.2), 
and  with  =  Q  =  0  ,  this  can  be  reduced  to 

v(e)  -  v(?3)  =  I  V 


1 


a.  (!i) 

'  K' 


9 


k+l 


^  \v  ft\ 


^  ^  /  i»»  \  ->  ic  i  w  I 

0,,  e,,. , 


'k 


1 

Wl 


(5.4b) 


Kj,(?)  ^cD(4-ej,)  = 


(5.4c) 


ffl(l)  =  e®log|4|  +  y'7{C)(5+C)®log|i+C|dC  (5.4d) 

o 


The  Integral  in  (5.1)  is  thus  replaced  by  a  finite  sum.  For 
a  given  choice  of  the  division  points  .  the  functions 
Hj^d)  depend  only  on  the  creep  function  y{t^) ,  and  may  be 
regarded  as  evaluated  once  for  all,  for  some  given  viscoelas¬ 
tic  material.  Then  (5.4a)  gives  the  displacement  for  any 
pressure  distribution  in  terras  of  its  values. 

The  trapezoidal  pressure  is  continuous,  but  has  dis¬ 
continuous  slope  dQ/d^  at  each  The  resulting  displace¬ 

ment  is  consequently  continuous  with  finite  slope  dv/d| 
everywhere.  But,  as  may  be  verified  from  (5.4a),  d  v/d4 
will  be  infinite  at  each  The  pressure  could  be  better 

approximated  by  polynom.ials  of  higher  degree,  which  would 
give  smoother  pressure  and  displacement.  But  the  added 
smoothness  does  not  justify  the  increased  complexity  in¬ 
volved.  It  is  better  to  use  the  simple  representation 
given  by  (5.3)^  taking  more  divisions  if  increased  accuracy 
is  desired. 
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WhcGH  th-G  pi-G55ui"c  Q(^)  is  known^  tfiG  toT/ai  ioao  is 


given  by 


N 


J  ^U)  d? 


Substituting  from  (5.3)  i'or  Q(4)>  this  becomes 


The  displacement  at  any  point  on  the  surface  can  be 
determined  from  (5.^a>b)  for  any  viscoelastic  solid.  The 
elastic  solid  is  covered  by  7(C)  =  0.  If  7(C)  is  in 
numerical  form,  determined  from  measured  creep  functions 
(as  described  in  Section  1.3)i  then  the  integral  in  (5.^d) 
would  have  to  be  evaluated  numerically.  If  7(C)  has  an 
appropriate  analytical  form,  a)(^)  can  be  expressed  analyti¬ 
cally.  For  example,  if  the  viscoelastic  material  can  be 

m 

represented  by  a  general  Voigt  solid,  so  that  7(C)  =  2  f. e 

i=l  ^ 

then 

m  f .  o  m  f. 

a)(?)  =  (1  +  S  gi)  r  login  +  ?(l+2  login)  S 

i=l  i  i=l  bjj^ 

(5.6) 


-biC 


+  [log|e|  -  e^i^  S^(-bj^e)]  2 


m  2f , 


i=l  b. 


where  E^(-x)  is  the  Exponential  Integral. 
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5»2  Numerical  Procedure  for  Moving  Contact  Problems 


The  moving  contact  problem  was  described  in  Chapter  III, 
and  the  general  procedure  for  obtaining  a  solution  was  out¬ 
lined  in  Section  3.2.  Values  of  p  and  A  (Fig.  5.1)  are 
chosen,  and  the  unknown  pressure  is  found  from  the  integral 
equation  (3.8).  Since  $  and  X  are  not  independent  in  a 
particular  smooth  contact,  it  is  necessary  to  include  an 
unknown  angle  of  tilt  a  along  with  the  prescribed  contact 
surface  w(^-p).  Then,  reasonable  values  of  p  and  X  can  be 
chosen  independently,  and  the  solution  will  give  the  neces¬ 
sary  tilt  for  this  to  be  a  possible  smooth  contact. 

For  a  numerical  solution,  the  unknown  smooth  pressure 
is  approximated  by  n  trapezoidal  elements.  The  number 
and  spacing  of  the  division  points  will  be  governed 

by  the  accuracy  desired  and  the  character  of  the  specific  % 

problem.  Following  the  procedure  of  Section  ^.1,  the  integral 
equation  (3.8)  becomes 

ic— 1 

0  <  5  <  (5.7) 

where  n.(^)  is  given  by  (5-^h-d)  and  f  is  arbitrary. 

There  are  n  unknown  quantities,  the  n-1  pressure  ordi¬ 
nates  =  Q(^j^)  ,  k  =  1,2,  ...,n-l,  and  the  tilt  a. 
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To  solve  for  these  unknowns^  the  method  of  collocation 

is  used.  If  (5.7)  is  satisfied  at  each  of  n  points  ^  =  1^, 

0  <  <  A  (j  =  1,2, ...,n),  there  will  be  n  equations  for 

the  n  unknowns.  Iich  of  the  unknowns  occurs  as  a  linear 
1 

factor  ,  so  the  problem  is  reduced  to  solving  a  set  of  simul¬ 
taneous  linear  algebraic  equations.  For  given  divisions  of 
the  contact  region  and  given  points  of  collocation  (^*) 

there  will  be  a  unique  solution  for  the  and  a.  With  this 
solution  the  left  side  of  (5*7)  will  exactly  equal  the  right 
side  only  at  the  points  ,  but  agreement  will  be  close 
throughout  the  contact  region  if  these  points  are  suitably 
chosen. 

It  is  convenient  to  choose  one  of  the  division  points 
as  the  reference  point  4^,  and  the  remaining  as  the 
collocation  points  .  The  collocation  procedure  then  in¬ 
sures  the  prescribed  displacement  at  each  end  of  the  contact 
region  and  at  the  n-1  intermediate  points  at  which  the 
pressure  is  determined.  Let 

^ jk  ~  ^k^ ^ J ^  ~  ^k^ ^a ^  *  k=  1, 2, ...,n-l  (5. 8a ) 

hn  =  (5.&;>) 

«j  =  w(5j  -  @)  -  -  e)  (5.80) 

This  is  the  principal  reason  why  the  small  tilt  a  is  in¬ 
troduced  as  an  unknown -- because  it  is  a  linear  term.  If  the 
tilt  were  prescribed  or  implicit  in  the  function  w(4-p),  6 
could  not  be  prescribed,  but  would  have  to  be  found  in  the  so¬ 
lution.  Since  p  is  not  a  linear  factor  in  general,  the  whole 
procedure  would  be  greatly  complicated. 
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>(here  J  =  0,1,..., n  except  j  =  a.  Then,  (5-7)  becomes 


1  n-1 
k=l 


\  *Jk  +  “  *jn  =  ' 


0, 1, . . . , n  except  j  =  a 

(5.9) 


This  is  the  set  of  n  equations  to  be  solved  for  the  n 
unknowns  Q,^{^  =  1,2,  ...,n-l)  and  a,  which  gives  an  approxi¬ 
mate  solution  of  the  integral  equation  (3.8).  The  smooth 
contact  conditions  (3.6a,b)  are  automatically  satisfied  in 
this  formulation  of  the  problem. 

Once  the  are  known,  the  total  load  N  is  given  by 
(5.5).  The  coefficient  of  friction  is 


1 

N 


^k 


*^k-l  (5.10) 


The  complete  numerical  procedure  for  solution  of  a  two- 

dimensional  moving  contact  problem  is  summarized  below; 

Given:  'y(C)"3  creep  function  of  the  base  material 
(see  Section  1.3) 

w( C) -prescribed  shape  of  contact  surface 


1) 

2) 

3) 


Choose  values  of  A,  p. 


Divide  the  contact  region  into  n  segments  at  the 
points  =  0,  •  •  •  •>  *  •  •  •»  Choose  one 

as  reference  point 

Determine  and  w^  from  (5.8a-c),  (5.^b,  c),  using 
(5-^8)  for  a  general  viscoelastic  material,  or  (5-6) 
for  a  Voigt  solid. 
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4)  Solve  the  set  of  n  simultaneous  linear  equations 
(5*9)  for  a  and  (k  =  1 ,2, . . . ,n-l) .  The  pressure 
distribution  is  then  given  by  (5. 3) --straight  lines 
connecting  the  ordinates  at  the  points  . 

5)  Find  the  total  load  N  from  (5*5)>  the  coefficient 
of  friction  from  (5* 10). 

6)  Displacements  at  any  points  on  the  surface  are  found 
from  (5.4a). 

For  each  pair  of  values  of  p  and  "K  this  procedure  gives 
one  solution,  with  particular  values  of  N  and  a  .  If 
either  the  actual  load  N  or  the  velocity  V  is  then  speci¬ 
fied,  the  other  is  determined  by  this  solution.  This  means 
that  if  certain  values  of  V,  N*,  and  a  are  wanted,  trial 
solutions  must  be  made  with  various  values  of  p  and  X  until 
the  desired  results  are  obtained. 

The  actual  numerical  calculations  are  straightforward > 
and  can  be  programmed  for  a  digital  computer.  Evaluation 
of  the  coefficients  will  depend  on  the  form  in  which 

7(C)  is  expressed.  For  a  Voigt  solid,  the  general  analytic 
form  (5-6)  can  be  used.  If  the  creep  functions  are  given 
numerically,  7(C)  can  be  found  as  described  in  Section  -.3. 
Then,  the  integral  in  (5.4d)  is  evaluated  numerically.  Com-^ 
puter  library  routines  for  solving  simultaneous  linear 
equations  are  usually  available. 

To  summarize,  this  numerical  procedure  yields  a  solu¬ 
tion  to  the  moving  contact  problem  by  representing  uhe  unknown 
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pressure  distribution  as  n  trapezoidal  elements.  The  pre¬ 
scribed  displacement  is  then  satisfied  at  certain  points  in 
the  contact  region,  leading  to  n  simultaneous  equations  to 
be  solved  for  the  unknown  pressure  ordinates.  The  principal 
advantages  of  the  method  are  its  simplicity  and  flexibility, 
and  that  it  can  be  used  for  any  viscoelastic  (or  elastic) 
solid  and  for  any  shape  of  contact  surface  subject  to  the 
conditions  of  smooth  contact.  In  the  next  sections,  numeri¬ 
cal  results  are  obtained  in  some  specific  applications. 

5.3  Numerical  Solution  of  the  Rolling  Cylinder  Problem 

The  numerical  method  of  Section  5.2  will  be  applied  to 
the  rolling  cylinder  problem  as  formulated  in  Section  4.1, 
where  it  was  indicated  that  P  is  to  be  found  and  it  is  not 
necessary  to  Introduce  a  .  The  physical  quantities  are  shovm 
in  Pig.  I  (at  the  end  of  the  chapter),  and  the  dimensionless 
quantities  are  defined  there.  It  is  convenient  to  choose 
^  =  ^  =  0  .  Then,  using  (4.1)  and  with  Q,  ’  =  q’(^,  )  , 

(5.7)  becomes 

s  %'  [nj^CD  -  !!^.(o)]  =  2pe  -  I®  ,  0  <  i  <-K  (5.11) 

Collocating  at  the  points  ^  ^  "  l,2,..-,n  ,  the  set  of 
equations  corresponding  to  (5*9)  is 

n-1  2 

s  \  +  2p(-ej)  =  -  .  J  =  (5-12) 
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where  is  given  by  (5.8a)  with  =  0.  Once  the  solu¬ 
tion  has  been  obtained ^  the  total  load  is  (from  5.5) 


^  ~  (^k+l"^k-l^ 


(5.13) 


k=l 


and  the  coefficient  of  friction  is  (from  4.4,  4.2a,b,  and 
5.3) 

^k+1  ^k  ^k-1 


n-1 


~  N*  \  ^^k+l"^k-l^ 


-  P  (5.1^) 


From  the  elastic  solution  and  the  existing  solutions 
for  standard  linear  solids  (Section  4.2,  Fig.  4.2),  it  is 
expected  that  the  pressure  distribution  will  be  smooth  for 
a  general  viscoelastic  solid.  This  suggests  division  of  the 
contact  region  into  equal  segments.  Then,  =  k0  ,  where 


^  “  ^k  “  ^k-1  “  n  ^  k  -  1,2, ...,n 
The  set  of  equations  (5.12)  becomes 


(5.-15) 


n-1 


2  Aj^  -  2P0J  =  -jV  ,  j  =  1,2,  ...n 


(5.16) 


where  now 


*Jk  “ 


(5.17a 


“k(5j)  =  ft2“(?r«k)  -  “(?j-?k-l)  -  “(?J-^k+l)]  t5.17b 
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and  cod)  is  given  by  (5.4d)  or  (5*6).  The  total  load  {5-13) 


becomes 


N 


26 


n-1 

2 

k=l 


and  the  coefficient  of  friction  (5.14)  becomes 


PX 


-  P 


(5.18) 


(5.19) 


Division  into  equal  segments  is  used  for  all  rolling  cylin¬ 
der  examples  given  in  this  chapter. 

A  particular  problem  for  a  given  material  [given  7(C)] 
is  fixed  by  choosing  the  number  of  divisions  n  and  the 
length  6  of  each  segment  (or  X  =  n0).  The  coefficients 

determined  from  (5.17a^t>)j  the  set  of  n  simul¬ 
taneous  equations  (5.16)  is  solved.  This  determines 
(k  =  l,2,...,n-l)  and  p.  From  (5.I8),  (5.19)  N*  and 
PX  are  found.  The  results  are  most  conveniently  expressed 
in  terms  of  the  quantities  defined  by  (4.5a-d).  Each  value 
of  A  (or  0)  corresponds  to  a  certain  V,  if  N*  is  specified, 
or  to  a  certain  N  ,  if  V  is  specified.  To  illustrate  the 
procedure,  some  specific  examples  are  given  in  this  and  the 
next  sections.  The  numerical  work  was  performed  on  a  digi¬ 
tal  computer^,  using  an  existing  library  routine  for  solving 
the  simultaneous  equations. 

The  IBK  7090  at  the  Stanford  University  Computation 
Center  (using  eight  significant  figures). 
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A  good  check  on  the  numerical  procedure  is  provided  by 
the  elastic  base,  for  which  7(C)  =  0*  The  known  analytic 


solution  is 

a'd)  =  ivtUv'T  (5.20) 

N  =  i  ^  ~  2 

and  the  pressure  distribution  is  symmetric  about  =  A/2. 
Numerical  results  for  n  =  10  and  n  =  40  are  given  in 
Table  5-1.  In  each  case  p  =  A/2  exactly  (to  eight  signifi¬ 
cant  figures),  and  symmetrically  located  values  of  are 
identical  to  at  least  six  significant  figures.  Values  of 
for  n  =  10  are  about  2^,  those  for  n  =  40  about  .55^, 
lower  than  the  exact  solution.  Values  of  are  less  than 
the  exact  value  of  k.8%  for  n  =  10,  1.2^  for  n  =  40. 

When  the  material  is  a  standard  linear  solid  with  con¬ 
stant  V  [i.e.,  7(C)  =  fe  results  of  the  numerical  method; 
can  be  compared  with  the  complete  analytic  solution  given  in 
Section  4.2.  A  few  such  examples  are  given  in  this  section 
to  indicate  the  accuracy  of  the  method.  Increasing  n 
beyond  a  certain  value  will  not  necessarily  mean  more  accur^- 
ate  results.  When  the  segments  are  very  short  (0  small, 
n  large),  increased  precision  (more  significant  figures)  is 
necessary,  requiring  careful  attention  to  the  details  of 
calculation,  and  perhaps  more  sophisticated  numerical 
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techniques.  Trial  solutions  for  several  values  of  n  will 
enable  selection  of  a  suitable  value  of  n  for  a  particular 
class  of  problems. 

Two  rolling  cylinder  problems  are  solved  with  n  =  10^ 
20,  and  4o.  In  each,  A  =  1  and  7(C)  =  ,  with  f  =  1 

in  one  example,  f  =  10  in  the  other.  Results  are  given  in 
Table  5-2,  and  pressure  distributions  are  shown^  in  Pig.  II. 
Numerical  values  are  consistently  improved  with  larger  n, 
the  error  being  roughly  halved  when  n  is  doubled.  Except 
near  the  ends,  the  pressure  distributions  for  all  three 
values  of  n  are  nearly  the  same.  On  the  basis  of  such 
examples,  n  =  20  was  chosen  for  use  in  all  subsequent  roll¬ 
ing  cylinder  examples.  Comparison  with  the  analytic  solu¬ 
tion  indicates  errors  of  2  to  5^.  The  computer  time  needed 
for  a  typical  solution  with  n  =  20  was  .1  to  .5  minutes. 

For  a  given  base  material  [given  7(C) ]>  the  complete 
range  of  rolling  cylinder  problems  is  covered  by  varying  A 
from  0  to  00.  To  illustrate  typical  results,  numerical  solu¬ 
tions  for  some  representative  values  of  A  are  given  in 
Table  5.3j  along  with  the  analytic  solution  for  comparison. 
In  these  examples,  7(C)  =  Numerical  values  com¬ 

pare  with  the  exact  solution  as  follov/s:  3  -  very  close 
agreement;  N'  -  numerical  solution  consistently  2  to  3^  low; 
PX  -  error  of  5%  or  less.  The  limiting  cases  A  =  0  and 

T 

Figures  labeled  with  Roman  numerals  are  at  the  end  of 
the  chapter. 


X  =  m  tb?*  ssni6  th6  6l8stic  Solution  (Tsble  4,1, 

Section  4.1),  These  limits  are  closely  approached  for 
7\  =  .01  and  X  =  100  respectively.  The  pressure  distri¬ 
butions  are  shown  in  B’ig.  IV.  For  A  =  .01  and  X  =  100  the 
pressure  is  very  nearly  symmetric,  while  for  X  =  1  it  is 
noticeably  non-symmetric .  For  large  and  small  X,  p/X 
approaches  .5,  while  the  lowest  value  of  p/X  occui’s  for 
X  =  1.  The  variation  of  coefficient  of  friction  with  veloc¬ 
ity  is  shown  in  Pig.  VI,  where  Rx/a^  is  plot tec  against 
VT/a^  (log  scale).  The  friction  curve  has  a  pronounced  peak, 
with  the  maximum  occurring  for  VT/a^  =  1.42  (x  =  1.6). 

The  results  of  this  set  of  examples  indicate  some  fea¬ 
tures  which  are  characteristic  of  any  viscoelastic  rolling 
cylinder  problem.  The  elastic  solution  is  symmetric  and 
gives  no  resisting  force.  A  viscoelastic  material,  due  to 
its  delayed  recovery  and  energy  dissipation,  gives  rise  to 
asymmetry  in  the  solution  and  a  resisting  force  (indicated 
by  x)*  The  asymmetry  is  evident  in  two  respects  (see  Fig. 
4.2):  1)  the  contact  region  is  displaced  upstream,  so  that 

its  center  (4  =  ^/2)  does  not  lie  directly  below  the  center 
of  the  cylinder  (4  =  P)j  this  means  p/X  is  less  than  1/2 
(the  elastic  value),  the  difference  indicating  the  amount 
of  a symmetry j  2)  the  pressure  distribution  is  not  symmetric 
in  the  contact  region j  the  maximum  pressure  occurs  upstream 
from  the  center  of  contact.  These  features  are  clearly  shown 
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in  the  examples  just  given^  particularly  for  -  1.  As  the 
velocity  goes  from  0  to  the  resistance  (for  a  given  total 
load)  increases  from  zero,  reaches  a  maximum  for  A  near  1, 
then  decreases  to  zero.  This  confirms  the  predictions  made 
in  Section  1.5  about  the  variation  of  x  with  velocity. 

The  examples  give  consistently  low  values  of  n  (about 
2  to  3^  for  n  =  20) .  This  error  is  due  principally  to  an 
inaccuracy  inherent  in  the  approximate  representation  of  the 
pressure  distribution.  The  pressure  elements  on  each  end 
are  triangular  (Fig.  5*2),  with  finite  slope  at  ^  -  0  and 
^  =  7^.  The  exact  pressure  distribution,  however,  will  have 
a  vertical  slope  at  the  ends"’, 
numerical  solution  tends  to  gi\e 
Q^'ardQ^  1  larger,  and 

the  other  slightly  smaller, 
than  their  exact  values.  This 
tendency  can  be  seen  in  Table 
5.1.  Even  with  this  correc¬ 
tion,  the  total  load  (area 
under  the  pressure  curve) 

found  in  the  numerical  solution  is  smaller  than  the  exact 
value.  The  di'^ference  is  indicated  by  the  siiaded  regions 
i.n  Pig.  5.2.  With  smaller  segments  (larger  n),  a  steeper 

1 

This  is  true  for  the  elastic  solution  (5-20),  and  is 
also  true  for  viscoelastic  solids  with  initial  elastic 
response . 


To  compensate  for  this,  the 


PIG.  5.2 
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s±op0  st  t-h?  0nd  is  possiblSj  thus  moi?^  closslj'^  I’spr’sssnt tng 
the  exact  pressure  distribution.  This  is  evident  in  Fig. 
Ill,  where  the  same  problem  is  solved  for  different  values 
of  n.  The  numerical  procedure  could  be  improved  somewhat 
by  using  unequal  divisions  of  the  contact  region.  Smalle’^ 
segments  could  be  used  near  the  ends,  since  that  is  where 
the  pressure  changes  rapidly. 

Comparison  of  the  numerical  and  analytic  solutions  for 
these  examples  shows  that  the  num.erical  method  with  n  =  20 
gives  results  within  a  few  per  cent  of  the  exact  values. 

This  suggests  that  the  same  procedure  could  give  equally 
satisfactory  results  in  problems  for  which  no  analytic  solu¬ 
tion  can  be  I'ound.  An  analytic  solution  is  available  only 
for  a  standard  linear  solid  with  constant  v.  Two  principal 
restrictions  of  this  model  arq  1)  similar  behavior  in  shear 
and  dilatation  [B(^)  =  J(C)#  thus  v  =  constant],  and  2)  a 
single  retardation  time.  Although  some  general  features  of 
rolling  cylinder  problems  are  exhibited  with  this  model, 
other  important  aspects  are  not  adequately  represented.  In 
the  next  section  some  examples  will  be  given  for  which  the 
material  behax^ior  is  not  so  restricted. 

5.4  Examples  of  Rolling  Cylinder  Solutions 

More  general  and  realistic  viscoelastic  behavior  is 
illustrated  by  simple  Voigt  solids,  each  having  several 
retardation  times  and  not  restricted  to  constant  v.  Then, 
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m  ,  ^ 

yiO  =  2  ,  with  m  >  2,  and  a^H)  is  again  given  by 

i=l  ^ 

(5.6).  The  numerical  procedure  of  Section  5.3  is  used,  with 
the  contact  region  divided  into  20  equal  segments. 

The  dilatation  behavior  of  certain  viscoelastic  materials 
can  sometimes  be  considered  as  elastic.  In  this  case,  B(0=  1 
and  V  varies  in  time,  v;ith  >  ''o*  simplest  such 

material  behaves  in  shear  as  a  standard  linear  solid,  with 
j(^)  =  1  +  f^(l-e“^)j  this  will  be  used  for  illustration. 

In  addition  to  f^,  the  initial  Poisson's  ratio  must  be 
specified.  For  this  set  of  examples,  f^  =  1  and  =  1/3^ 
and  the  final  Poisson's  ratio  is  =  .412.  Then,  from 
(1.28b) 

7(C)  =  .750  e"^  +  .027778  e"  (5.22) 

Results  for  representative  values  of  A  are  given  in  Table 
5.4.  The  limiting  cases  A  -  0  and  A  »  are  found  from 
Table  4.1  (f  =  .765  in  these  examples^  The  pressure  dis¬ 
tribution  for  A  =  1  is  plotted  in  Fig.  V.  The  variation  of 
friction  coefficient  with  velocity  is  shown  in  Fig.  VI. 

Some  corresponding  results  for  the  material  having  the 
same  behavior  in  shear,  but  constant  v,  are  shown  in  Table 
5.3,  and  Figs.  V,  VI.  Comparison  shov;s  the  results  for  the 
two  materials  (differing  only  in  dilatational  behavior)  to 
be  nearly  the  same.  This  indicates  that  in  rolling  cylinder 
problems  the  exact  nature  of  the  behavior  in  dilatation  is 
not  a  critical  feature.  With  elastic  dilatation,  the  effects 


h 

t 


of  viscoelafjticity  (minimum  p/A,  maximum  Rx/a^)  are  again 
most  pronounced  for  A  in  the  neighborhood  of  1,  but  the  mag¬ 
nitude  of  these  effects  is  somewhat  less  than  for  the  material 
with  constant  v. 

In  all  the  examples  considered  thus  far,  the  viscoelas¬ 
tic  material  has  been  a  model  v/ith  a  single  retardation  time 
in  shear.  For  such  a  material,  most  of  the  transition  from 
initial  to  final  elasticity  takes  place  in  a  short  time  span 
(on  a  log  scale),  about  one  decade  on  each  side  of  the  re¬ 
tardation  time^.  This  means  most  of  the  viscoelastic  be¬ 
havior  is  concentrated  in  this  short  time  span.  In  the  roll¬ 
ing  cylinder  problem,  A  is  the  dimensionless  time  needed  for 
a  given  element  of  material  to  traverse  the  contact  region. 
When  A  is  near  1,  the  element  undergoes  most  of  its  transi¬ 
tion  from  initial  to  final  behavior  during  the  time  of  con¬ 
tact,  and  the  viscoelastic  effects  are  very  pronounced. 

This  is  somewhat  analogous  to  resonance  in  damped  vibrations. 
When  A  is  less  than  about  .1,  the  contact  time  is  too  short 
for  little  more  than  the  initial  elastic  behavior  to  be 
effective.  Likewise,  when  A  is  larger  than  about  10,  the 
contact  time  is  long  enough  that  the  final  elastic  behavior 
predominates.  Thus,  the  asymmetry  of  the  solution  is  most 
pronounced  for  A  in  the  neighborhood  of  1,  and  the  friction- 
velocity  curve  has  a  definite  peak  in  this  neighborhood. 

^  This  retardation  r-ime  is  used  to  form  the  dimension¬ 
less  time  variable  t  =  t/T,  so  it  has  the  value  r  -  1, 
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These  aspects  of  the  results  indicate  that  the  model 
with  a  single  constant  retardation  time  is  very  artificial. 
Rather  than  discrete  retardation  times,  actual  materials 
appear  to  have  a  continuous  distribution  over  a  large  time 
span.  The  behavior  of  an  actual  material  over  the  full 
range  of  contact  times  &/Y  can  be  adequately  represented  by 
a  single  model  of  discrete  retardation  times  only  if  there 
is  a  broad  spectrum  of  such  times.  The  viscoelastic  effects 
will  be  significant  over  a  much  wider  range  Ox’  A,  but  the 
magnitudes  of  these  effects  will  be  considerably  decreased. 

This  will  be  illustrated  by  a  simple  model  with  constant 
V  and  five  redardation  times,  t  =  .01,  .1,  1,  iO,  and  100. 

The  retardation  times  have  equal  weight  in  the  creep  be¬ 
havior,  and  the  final  elasticity  is  the  same  as  for  the 
set  of  examples  of  the  last  section  (f-1).  Thus, 

J(t)  =  2  - 

7(C)  =  .2(  e‘^+  10e‘^°^+ 

(5.23) 

The  results  for  A  =  1  are  given  in  Table  5.5^  3nd  the  pres¬ 
sure  distribution  shown  in  Fig.  V.  For  comparison,  these 
include  also  results  previously  given  for  two  materials  with 
a  single  retardation  time  in  shear.  A  few  points  on  the 
friction-velocity  curve  are  shown  in  Fig.  VI.  Although 
only  five  discrete  times  were  used,  the  results  are  in 
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striking  contrast  to  the  single  time  model.  For  X  =  Ij 
the  5-tim3  model  gives  a  much  more  symmetric  solution,  with 
pA  much  closer  to  .5.  The  friction  curve  is  broad  and 
smooth,  with  no  pronounced  peak.  The  maximum  again  occurs 
near  ?\  =  1,  but  the  curve  is  relatively  flat  over  the  five 
decades  covering  the  retardation  times.  The  maximum  value 
of  the  coefficient  of  friction  is  considerably  reduced, 
being  only  about  1/5  of  the  peak  for  the  single  time  model 
with  constant  v. 

These  results  and  other  examples  lead  to  the  general 
conclusion  that  a  single-retardation-time  model  is  not  ade¬ 
quate  for  describing  quantitatively  the  behavior  of  actual 
materials  when  subject  to  rolling  contact  over  a  wide  range 
of  contact  times.  Actual  materials  with  a  broad  spectrum  of 
retardation  times  will  not  produce  the  pronounced  viscoelas¬ 
tic  effects  of  a  single-time  model.  Instead,  these  examples 
indicate  that  the  effects  will  be  much  reduced  in  magnitude 
and  spread  over  a  broader  range  of  contact  times.  This  is 
evident  also  in  the  example  with  elastic  dilatation.  The 
difference  in  shear  and  dilatation  behavior  introduces  a 
secondary  retardation  cime  of  relatively  small  weight  (see 
5.22).  But  even  this  produces  a  noticeable  decrease  in  the 
friction  and  a  more  symmetric  solution  (see  Table  5.5). 

For  a  given  final  elasticity  f  and  contact  time  a  model 
with  several  retardation  times  will  have  a  more  symmetric 
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solution  and  lower  friction  than  the  corresponding  single - 
time  model.  There  are  definite  computational  advantages  in 
using  a  model  to  represent  an  actual  material^  but  the  model 
must  have  at  least  several  retardation  times  for  the  results 
to  be  at  all  realistic. 

The  discussion  so  far  has  been  concerned  principally 
with  those  aspects  of  viscoelastic  behavior  which  are  influ¬ 
enced  by  the  retardation  times.  But  also  important  is  the 
magnitude  of  the  change  from  initial  to  final  behavior. 

This  is  most  conveniently  indicated  by  the  ratio  of  final 
to  initial  deformation  in  a  creep  test.  In  particular,  for 
moving  load  or  contact  problems  this  ratio  is  given  by  1  +  f, 

CX) 

where  ^  ~  7(C)^C  •  Small  values  of  f  indicate  nearly  elas¬ 

tic  behavior.  Large  values  of  f  mean  that  the  final  response 
is  much  greater  than  the  initial  response.  The  limiting 
case  as  f  -*■  »  is  a  material  with  no  initial  elasticity. 

This  means  no  free  spring  in  a  mechanical  model  (Fig.  1.3). 
The  simplest  such  model  is  a  Kelvin  solid,  consisting  of  a 
spring  and  dashpot  in  parallel  and  having  one  retardation 
time. 

To  Illustrate  the  effects  of  the  value  of  f,  a  standard 
linear  solid  with  constant  v  is  used,  thus  7(C)  =  fe~^. 

Pig.  Ill  shows  the  pressure  distribution  for  f  =  10,  A  =  1. 
There  is  a  slg'^ificant  peak  in  the  pressure  very  close  to 
the  upstream  enf^  jf  the  contact  region.  The  results  for 
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various  values  of  f,  with  the  same  contact  time  ?\  =  1,  are 
given  in  Table  5.6,  and  the  corresponding  pressure  distri¬ 
butions  in  Pig.  VII.  The  case  f  =  0  is  the  elastic  solution 
which  is  symmetric.  As  f  is  increased,  the  viscoelastic 
effects  become  more  and  more  pronounced;  p  decreases,  indi¬ 
cating  the  contact  region  is  displaced  more  and  more  upstream; 
the  maximum  pressure  increases  and  the  peak  shifts  closer  to 
the  upstream  end;  the  friction  increases  (although  the  veloc¬ 
ity  remains  nearly  constant) .  For  the  limiting  case  as 
f  -♦  00  ,  the  Kelvin  solid,  there  is  an  exact  analytic  solu¬ 

tion  given  in  Section  4.2  (4.24-28).  In  this  limiting  case, 
the  pressure  peak  is  infinite  and  is  at  the  front  end  of  the 
contact  region.  Since  the  Kelvin  solid  has  no  initial  elas¬ 
ticity,  the  results  in  Table  5.6  and  Fig.  VII  are  based  on 
the  final  elastic  constants  (l-v^)/li^  to  allow  comparison. 

These  results  indicate  that  for  large  f  the  viscoelastic 
effects  (asymmetry,  friction  coefficient)  will  be  much  in¬ 
creased  over  the  effects  for  small  f.  In  particular,  there 
will  be  a  significant  peak  in  the  pressure  distribution  near 
the  upstream  end.  The  presence  of  this  peak  is  due  to  the 
small  initial  response  relative  to  the  final  response  when 
f  is  large.  The  deformation  near  the  front  of  the  contact 
is  produced  principally  by  the  initial  elastic  response, 
while  farther  downstream  the  longer  time  of  contact  pro¬ 
duces  a  much  greater  response.  Thus,  a  very  high  pressure 
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comparable  in  magnitude  to  the  deformation  in  the  rest  of 
the  contact  region.  As  f  gets  larger^  a  greater  pressure  is 
needed,  so  the  peak  increases.  In  the  limit  f  ,  there 
is  no  initial  response,  so  an  infinite  pressure  is  needed. 

The  example  used  here  is  a  situation  most  favorable  to  pro¬ 
ducing  a  peak.  The  time  of  contact  (X=l)  coincides  with  the 
single  retardation  time  of  the  material,  so  the  viscoelastic 
effects  are  exaggerated  as  explained  earlier.  With  this 
same  material,  if  'K  »  1,  the  contact  time  is  long  enough 
so  the  final  elastic  response  is  effective  over  nearly  all 
the  contact  region.  In  this  case  the  results  are  nearly  the 
same  as  the  elastic  solution  for  all  values  of  f.  For 
A  «  1  ,  the  short  contact  time  means  only  the  initial  elas¬ 
tic  response  is  effective  over  most  of  the  region.  The  solu¬ 
tions  are  thus  close  to  elastic  again  when  f  is  not  too 
large.  However,  for  very  large  f  the  initial  response  is 
very  small  and  not  significant,  so  the  solutio.n  approaches 
that  of  the  Kelvin  solid. 

The  significance  of  several  rather  than  just  one  retar¬ 
dation  time  in  the  material  behavior  has  already  been  dis¬ 
cussed.  This  is  particularly  evident  when  f  is  large.  As 
an  example,  the  results  for  the  5-retardation-time  model  used 
previously,  this  time  with  f  =  100  and  A  =  1,  are  given  in 
Table  5*6.  The  solution  is  considerably  more  s^nranetric  and 
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tne  friction  coefficient  much  less  than  for  the  corresponding 
single-time  model.  The  pressure  distribution  is  very  nearly 
the  same  shape  as  the  f  =  1  curve  in  Fig.  VII  (with  slightly 
greater  ordinates)  and  there  is  no  pronounced  peak. 

Actual  materials^  such  as  polymers,  may  have  values  of 
f  as  high  as  10^.  However,  when  f  is  very  large  there  are 
difficulties  in  the  numerical  calculations  due  to  the  sharp 
peak.  Many  divisions  of  short  length  are  needed  in  the 
vicinity  of  the  peak.  This  is  evident  in  Fig.  Ill,  where 
results  for  different  values  of  n  can  be  compared.  But,  as 
mentioned  in  Section  5.3>  other  difficulties  arise  when  very 
short  segments  are  used.  Attempts  were  made  to  estimate  the 
peak  by  further  theoretical  considerations,  and  so  remove 
it  from  the  numerical  calculations,  but  without  success. 

3.5  Numerical  Solution  of  the  Nearly  Flat  Punch  Problem 

As  an  illustration  of  another  application  of  the  numeri¬ 
cal  method  for  moving  contacts,  the  problem  of  moving  "nearly 
flat"  punch  is  considered  in  this  section.  A  "flat  punch" 
is  8  rigid  (two-dimensional) 
body  with  a  plane  surface  at 
some  angle  a  pressed  into  the 
viscoelastic  base  (Fig.  5-3) 
by  a  force  N*.  If  the  punch 
is  perfectly  flat,  the -contact 
will  be  a  straight  line  with  ‘ 


FIG.  5.3 
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sure  will  be  infinite  and  the  slope  of  the  deformed  surface 
will  be  discontinuous  at  one  or  both  ends  of  the  contact. 

Real  corners^  however,  will  be  slightly  rounded,  not  perfectly 
sharp.  The  pressure  at  the  corners  may  be  high,  but  finite, 
and  the  contact  will  be  smooth.  The  important  question  is 
how  the  peak  pressures  depend  on  the  shape  of  the  rounded 
corners.  Such  a  "nearly  flat  punch"  may  score  or  gouge 
the  base  surface  if  peak  pressures  are  too  high.  Avoidance 
of  such  damage  is  a  practical  problem  of  frequent  occurrence. 

Some  aspects  of  the  moving  contact  of  a  nearly  flat 
punch  will  be  investigated  The  corners  will  be  assumed  to 
be  circular  arcs  (radius  R),  with  the  length  of  the  rounded 
portion  very  small  relative  to  the  total  length  of  contact 
(Fig.  5*^).  The  problem  is  then  treated  as  one  of  smooth 
contact.  The  length  d  of  the  flat  portion  is  fixed.  Then, 
the  length  of  the  rounded  corner  in  contact  on  the  left  end 
is  b,  and  for  a  total  contact  length  £  the  length  of  the 
rounded  corner  on  the  right  end  is  ^-d-b  =  r^b,  with 
b  «  i  .  In  general  7^  1  because  the  contact  is  not 
symmetric , 

The  elastic  solution  for  smooth  contact  of  a  nearly 
flat  punch  (reported  elsewhere)  shows  that,  while  the 
pressure  is  everywhere  finite  and  is  zero  at  the  ends,  very 
sharp  peaks  occur  just  inside  the  ends  of  contact.  As  the 
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length  of  the  corners  becomes 


smaller  (i.e.,  b/d  -♦  O),  the 


peak  gets  sharper,  and  in  the 


limit  the  solution  becomes  that 


for  a  perfectly  flat  punch. 


This  suggests  that  the  visco¬ 


elastic  solution  may  also  have 


sharp  peaks  in  the  pressure 


FIG.  5.4 


distribution  near  the  ends. 


The  viscoelastic  oroblem  is  treated  as  outlined  in 


Section  5.2.  Lengths  are  expressed  as  dimensionless  ratios 


to  VT: 


-  iL. 
P  -  VT 


Ti  =  (a-6-P)/P 


(5.24) 


The  other  dimensionless  quantities  given  by  (3.5a-c^)  ^re 
used.  The  shape  of  the  contact  surface  (with  the  circular 
arcs  approximated  by  parabolas)  is  given  by 


w(i-p)  =  -  ^  ,  0  <  e  <  p 


,  P  <  5  <  P  +  6 


-  Ip  (5-e-<5)®  ,  P  +  6  <  e  < 


The  n  collocation  points  are  c  sen,  and,  taking  ^  =  P 


(Fig.  5.1)^  w(^Q-p)  =  0.  Then,  from  (5-9)^ 


“Jk  "  ",in  = 


where 


0, 1, . . . , n,  except 

ra 

II 

(5.25) 

n-1 

k=l 

(p  =  2pa  t  and 

wj  s  2pv.j  =  -  ,  0  <  <  p 

=  0  >p<lj<^  +  p  (5.26) 

=  -(?j-p-a)^  ,  p  +  6  <  |_j  <  ^ 


The  k.y^  are  given  by  (5.8a,b)  and  (5.^0-d)  with  =  p. 

For  a  particular  problem^  values  of  p  and  "K  are  chosen,  thus 
determining  t)  (6  is  a  given  quantity).  The  solution  then 
gives  the  pressures  Qj^  and  the  tilt  0.  To  obtain  a  partic¬ 
ular  value  of  0,  for  example,  several  trials  with  different 
A  and/or  p  may  be  necessary 

For  illustration,  the  results  of  one  set  of  examples 
are  given.  A  standard  linear  solid  with  constant  v  is  used, 
so  7(C)  =  e~^.  In  these  examples,  5  =  .99  ^  =  .005j 

the  only  variation  is  in  t]  (thus  A  varies  slightly).  The 
division  of  the  contact  region  and  the  choice  of  n  is  based 
on  the  expectation  of  sharp  peaks  near  the  ends  and  rela¬ 
tively  little  variation  near  r.he  center.  After  some  pre¬ 
liminary  trials,  n  =  22  was  chosen.  A  typical  division  is 
shown  schematically  for  t]  =  I.67: 
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(not  to  scale) 


The  results  for  several  values  of  t)  are  given  in  Table 
5.7i  and  some  pressure  distributions  are  shown  in  Pig.  VIII. 

The  sharp  peaks  in  the  pressure  are  present  as  anticipated, 
and  the  expected  asymmetry  is  shown  by  the  greater  height 
of  the  upstream  peak.  Another  expectation  is  verified  by 
the  regular  change  in  tilt  as  rhe  upstream  corner  is  lenghened, 
other  factors  being  constant.  The  tilt  is  nearly  zero  when 
T\  =  1.67>  and  the  body  rotates  clockwise  as  ii  is  increased. 

The  pressure  distribution  is  very  little  changed  near  the 
downstream  end,  but  there  is  significant  change  near  the 
upstream  end  as  tj  is  varied.  The  peak  pressure  and  total 
load  N’  increase  with  t],  and  the  line  of  action  of  the  verti¬ 
cal  resultant  moves  upstream.  For  the  zero  tilt  example 
(ti  =  1.67),  the  line  of  action  of  N'  is  at  ^  =  .60. 

These  examples  indicate  the  general  suitability  of  the 
numerical  method  for  viscoelastic  moving  contact  problems 
of  this  type.  Even  with  a  relatively  small  n,  the  important 
features  are  evident  in  the  solution.  The  method  can  of 
course  be  applied  for  mere  general  viscoelastic  materials 
if  desired.  The  shapes  of  the  rounded  corners  can  also  be 
easily  changed.  For  example,  the  stress  peaks  could  perhaps 
be  made  more  nearly  equal  by  using  different  radii  at  the 
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two  ends.  It  is  interesting  to  note  that  (within  the 
limitations  of  small  displacem.ents  and  smooth  contact)  the 
radius  of  the  rounded  corner  enters  only  as  a  scale  factor 
in  the  pr-  3sure  and  load  expressions.  Thus,  if  b,  d,  and  T]b 
(Pig.  5.4)  are  kept  the  same,  the  pressure  distribution 
expressed  by 


1-y 


o  1 

TT 


q(x) 


remains  unchanged  if  the  radius  R  is  changed.  However, 
the  total  load  changes  in  proportion  to  R  ,  so  the  actual 
magnitude  of  the  pressure  q(x)  changes  similarly. 
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(ao/2N*)q  (Oo/2N*)q 


FIG.  II  ROLLING  CYLINDER  PRESSURE 
DISTRIBUTIONS  FOR 
VARIOUS  n  ,  f  =  1. 


FIG.  Ill  ROLLING  CYLINDER  PRESSURE 
DISTRIBUTIONS  FOP- 
VARIOUS  n  ,  f  -  10. 
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(ao/2N*)q  fao/2N’^)q 


C/X 


FIG.  IV  ROLLING  CYLINDER  PRESSURE 

DISTRIBUTIONS  FOR  VARIOUS  X  . 


FIG.  V  ROLLING  CYLINDER  PRESSURE 
DISTRIBUTIONS  FOR 
VARIOUS  MATERIALS. 
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Rx/0 


10**  10**  10*'  i  10  100  1000 


VT/ao(log  scole) 

FIG.  VI  ROLLING  RESISTANCE  VS.  VELOCITY 
FOR  ROLLING  CYLINDER. 


FIG.  VII  ROLLING  CYLINDER  PRESSURE 
DISTRIBUTIONS  FOR 
VARIOUS  VALUES  OF  f  . 
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FIG.  VIII  PRESSURE  DISTRIBUTIONS  FOR 
NEARLY  FLAT  PUNCH. 
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TABLE  3.1 


ROLLING  CYLINDER  ON  ELASTIC  BASE, 
NUMERICAL  AND  ANALYTIC  SOLUTIONS 


Values 

of  for 

A=10 

Numerical 

Ana lyti 

40^A 

o 

1 — ( 

II 

c: 

n=40 

1 

.508 

.497 

2 

.675 

.693 

3 

.825 

.839 

4 

.977 

.942 

.955 

6 

1.125 

1.13^ 

8 

1.234 

1.263 

1.271 

10 

1.368 

1.379 

12 

1.429 

1.449 

1.460 

14 

1.509 

1.518 

16 

1.529 

1.550 

1.560 

18 

1.575 

1 . 584 

20 

1.561 

1.583 

1.590 

N* 

23.80 

24.71 

25.00 

pa 

.5000 

.5000 

.5000 
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QiABLE  3-2 


ROLLING  CYLINDER  RESULTS  FOR  VARIOUS  n 
7(C)  =  fe"^;  -K=l 


Numerical  Ana  lytic 

n=10  n=20  n=40 


PA 

.4025 

.4014 

.4008 

.4003 

N' 

.2048 

.2097 

.2121 

.2144 

px 

Q'max. 

.1213 

.1350 

.1231 

.1359 

.1239 

. 1364  . 

.1249 

i  * 

PA 

.1981 

.1965 

.1958 

.1952 

N' 

.0697 

.0712 

.0718 

.0724 

f. 

px 

Q ’max. 

.4257 

.0835 

.4291 

.0719 

.4308 
.  0665 

.4322 

X  • 

Corresponding  pressure  distributions 
Figures  I  and  II. 

are  shown 

TABLE  5.3 


ROLLING  CYLINDER  RESULTS  FOR  VARIOUS  >v 


7(0=e"^ 


A 

pa 

N‘A^ 

PX 

VT/a^ 

I^/ao  ( 

a^/2N*)q  max 

0 

.5000 

.2500 

0 

2/>^ 

0 

.3183 

.01 

.4932 

.2500 

. 0000670 

200.0 

.01340 

.  4934 

.2440 

. 0000692 

202.4 

.01402 

.3186 

.1 

.4612 

.2491 

.00380 

20.0 

.0760 

.  4620 

.2430 

.00410 

20.28 

.0832 

.3180 

1 

.4003 

.2144 

.1249 

2.16 

.270 

.4014 

.2097 

.1231 

2.184 

.  2688 

.2967 

10 

.4768 

.1377 

.4292 

.270 

.116 

.4768 

.13^3 

.4276 

.2729 

.  1167 

.2273 

100 

.  4975 

,1263 

.4925 

.0281 

.0138 

.4975 

.1227 

.4828 

.0285 

.01378 

.2254 

00 

.5000 

.1250 

.5000 

2.828A 

0 

.2250 

A 

-  analytic  solution 

.,  Section 

4.2 

N  -  numerical  solution.  Section  5-3 
*  -  limiting  case,  exact  solution 

Pressure  distributions  are  shown  in  Figure  III. 
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TABLE  5-^ 


NUMERICAL  RESULTS  FOR  ROLLING  CYLINDERS 


Standard  linear  solid  in  sliear,  elastic  dilatation 

7(0  =.750e"^  +  .02778e~^*^^^^ 

=.333  ,  =.412  ,  f  =.765 


-K 

-JZl 

N'A^ 

^0 

PX 

BxAo 

0 

.  500 

.250 

00 

0 

0 

.318  (exact) 

.01 

.495 

.244 

202 

. 000054 

.0109 

.319 

.1 

.471 

.243 

20.3 

.00317 

.0643 

.318 

.5 

.430 

.233 

4.14 

.  0403 

.167 

.312 

.75 

.423 

.225 

2.81 

.0692 

.195 

.306 

1.0 

.422 

.217 

2.15 

.0976 

.210 

.301 

1 . 6 

.427 

.200 

1.40 

.157 

.219 

.289 

2.0 

.  433 

.192 

1.14 

.189 

.215 

.282 

5.0 

.464 

.162 

.496 

.308 

.153 

.250 

10 

.480 

.150 

.258 

.364 

.0940 

.242 

100 

.498 

.139 

.0268 

.413 

.0111 

.240 

1000 

.  500 

.138 

.  0027 

.415 

.0011 

.240 

00 

.  500 

.142 

0 

.432 

0 

.240  (exact) 
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TABLE  5-5 


ROLLING  CYLINDER  RESULTS  FOR  VARIOUS  MATERIALS 


A  =  1 


g)  Standard  linear  solid, 
constant  v 

b)  Standard  linear  solid 

in  shear,  elastic 
dilatation  (v^  =-333) 

c)  Five  retardation  times 

in  shear,  constant  v 


J(0  = 

7(0  = 

J(C)  = 
7(0  = 


1  +  (1  -  e'^) 

e 

1  +  (1  -  e"^);  B(C:  -  1 
.750e~^  +  .02778e'~‘®^^^ 

see  equation  (5* 23) 


pa 

N'/A^ 

.  9X 

a) 

.4014 

.2097 

.1231 

b) 

.4215 

.2171 

.0976 

c) 

.4694 

.1738 

'402 

2.184 

.269 

.297 

2.146 

.210 

.301 

2.399 

.096 

.268 

72 


TABLE 


b 


ROLLING  CYLINDER  RESULTS  FOR  VARIOUS  f 


A  =  1  in  all  cases 

Results  are  based  on  the  final  elastic  response:  a^,.=a_yG.+f 


f 

&/-K 

il+flN' 

-PX.. 

VT/a^ 

Rx/Sf 

Single 

retardation  time 

0 

.  500 

.25c 

0 

2.000 

0 

.318  (exact) 

1 

.401 

.419 

.123 

1.546 

.190 

.420 

10 

.197 

.790 

.429 

1.129 

.434 

.892 

100 

.149 

.872 

.  524 

1.06^ 

.559 

3.08 

00 

.143 

.893 

.538 

1.058 

.569 

“  (exact) 

Five 

retardation  times 

1 

.469 

.348 

.040 

1.696 

.  068 

.379 

100 

.401 

.576 

.133 

1.319 

.175 

.515 

TABLE 

3^1 

NUMERICAL  RESULTS  FOR 

NEARLY  FLAT  PUNCH 

p  =.005,  <5  =.99,  7(C)  =  e"‘=  ^ 

Q'  X  10^^ 


(p  xlO^ 

T 

\ 

4 

N’xlO 

left  max 

right  max 

min 

-2. 52 

1.0 

1.000 

2.52 

11.9 

12.0 

.71 

-  .32 

1.6 

1.003 

3.56 

12.0 

19.1 

.91 

.05 

1.67 

1.00335 

3.74 

12.0 

20.1 

.93 

.53 

1.8 

1.004 

3.97 

11.9 

pi  c; 
da. .  p 

.96 

1.33 

2.0 

1.005 

4.36 

12.0 

23.6 

.99 

6.13 

3.0 

1.010 

6.66 

11.9 

35.7 

1.23 
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